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1. SCOPE 


This report documents the software which has been deve'ioped in response to 
Action Document 63-21-37-3343-01, "Sampling Scheme for Early Warning Alarm.' 
The total software system referred to as the Patch/Skip Sampling System is 
described in detail. 
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2. SYSTEM OVERVIEW 


The programs contained in this document are intended to provide a way to 
extract a sample from a full-frame scene and summarize it in a useful way. 

The sample in each case was chosen to fill a 512-by-512 pixel (sample-by-line) 
image since this is the largest image that can be displayed on the Integrated 
Multivariant Data Analysis and Classification System (IMDACS) system. T!iis 
sample size provides one megabyte of data for manipulation and storage a. id 
contains about 3 percent of the full-frame data. 

Two extraction strategies (methods) are represented in the programs. In the 
skip method, every sixth pixel from every sixth line is taken to show a 
sampled image which preserves the gross structure of the full-frame image. In 
the patch method, 32- sample by 32- line pixel squares (patches) are taken from 
the full-frame image on a 16-by-16 square grid. These 256 squares of data can 
provide a summary of the full frame while retaining the full resolution field 
structure for analyst interpretation. 

In both cases, the 512-by-512 pixel image is further summarized by the patch 
image processor. This processor computes means for 256 32-by-32 pixel squares 
which constitute the 512-by-512 pixel image. For each square, the screening 
algorithm is run and the channel means are computed along with a screening 
summary of the rejected pixels. The mean vegetation index, truncated at 0, is 
computed for values GREEN, AVI, PVI, TVI7, and LAI. The values of GIN5, 

GIN15, and GIN30 are also computed. These may be computed with or without 
X-STAR haze correction (ref. 1). 

Thus, 256 measurements for each of 8 vegetation indexes are available over a 
100-mile (161-kilometer) square, computed in 2 waysi 

With these data from several segments, the vegetation indexes and the two 
modes of sampling the full-frame image are compared. 

Flow diagrams of Subsystems 1, 2, and 3, are presented in figures 2-1, 2-2, 
and 2-3. 
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Figure 2-1.- Flow diagram of Subsystem 

















SKIP 



Figure 2-2.- Flow diagram of Subsystem 


















Figure 2-3.- Flow diagram of Subsystem 3 










3. SUBSYSTEM 1 - READTP 


The READTP program is the main routine of the patch image processor. It is 
the driver routine for reading Landsat 1-, Landsat 2-. or Universal-formatted 
tape. In addition, it is the driver routine for creating the patch image in 
the IMDACS form. 


CALLING SEQUENCE : PROGRAM READTP 

CALLING ARGUMENTS : Not applicable. 

FILES ; 

File Usage Record format 

RPARAM.OAT User interface (17X,I5./,17X,I5,/12X,I3./,10X,6A1) 

??????. IMH IMDACS header file Refer to IMDACS documentation, volumes 1 and 2 
??????. IMO IMDACS data file Refer to IMDACS documentation, volumes 1 and 2 


COMMON BLOCKS ; The following abbreviations are used in tables throughout 
this document. 

A » alphanumeric 

I = integer 

I/O = input output 

L = logical 

R » real 


Common blocks for the READTP routine are /COMH/, /LABEL/, and /HIST/. 
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/COMH/COMMON BLOCK 


Name 

Relative 
word number 

Usage 

Jm. 

ST 

1 

I/O 

I 

SE 

2 

I/O 

1 

LS 

3 

I/O 

I 

LE 

4 

I/O 

I 

NOCHAN 

5 

I/O 

I 

NOSAMP 

6 

I/O 

I 

SCNID 

7-12 

Header 

A 

SUNEL 

13 

Header 

I 

SUNAZ 

14 

Header 

I 

SNSHD 

15 

Heujier 

I 

NREV 

16 

Header 

I 

EXYR 

17 

Header 

I 

EXDAY 

18 

Header 

I 

DAY 

19 

Header 

I 

HR 

20 

Header 

I 

MIN 

21 

Header 

I 

SEC 

22 

Header 

I 

SEQNO 

23-24 

Header 

I 

FLAT 

25-26 

Header 

R 

FLONG 

27-28 

Header 

R 

NLAT 

29-30 

Header 

R 

NLONG 

31-32 

Header 

R 


Definition 

Starting pixel number on computer com 
patible tape (CCT) . 

Ending pixel number on CCT. 

Line start number on full frame. 

Lins end number on full frame. 

Number of channels on tape. 

Number of pixels per scan line per 
channel . 

Scene identification. 

Sun elevation angle. 

Sun azimuth angle. 

Sensor heading. 

Orbit revolution number. 

Exposure year. 

Exposure day (Julian). 

Scene identification day. 

Scene identification hour. 

Scene identification minute. 

Scene identification seco-nd. 

CCT sequence number (N of M), 

Center of the full frame (latitude). 
Center of the full frame (longitude). 
NADIR (latitude). 

NADIH (longitude). 
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/LABEL/COMMON BLOCK 


Relative 

Name word number Usage Type Definition 
NAME 1-3 Error A Name of routine. 

MISNO 4 Header I Mission number. 

E 5 Header I Mission number. 

/HIST/COMMON BLOCK 

Relative 


Name 

word number 

Usage 

Jm. 

Defi nition 

H4 

512 

Array 

R 

Histogram for band 4 

H5 

1025 

Array 

R 

Hi stogram for band 5 

H6 

1537 

Array 

R 

Histogram for band 6 

H7 

2049 

Array 

R 

Histogram for band 7 


SUBROUTINES CALLED: 


Reference 


Name 

Definition 

(section) 

HDREAD 

Header record read. 

3.1 

INTOCT 

Integer to OCTAL. 

3.2 

HDRFLN 

Header file. 

3.3 

HDT 

High density full-frame format. 

3.4 

LANDS AT 

Landsat (1 and 2) format. 

3.5 

HDRHIS 

Header histogram. 

3.6 

CALLED 

BY : Not applicable. 
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LOCAL VARIABLES; 


Name Definition 

BFSZ Buffer size in bytes. 

BUF Buffer for image data. 

OELTAL Change in center point line. 

DELTAP Change in center point pixel. 

IE Error flag. 

NOML Nominal center point line. 

NOKP Nominal center point pixel. 

NPOINT Array containing nominal center point in degrees. 

STARTL Line in which extracting should start. 

STARTP First pixel to be extracted during patch extraction from full frame. 
STOPP Last pixel to be extracted during patch extraction from full frame. 
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3.1 UNIT 1 - HDPEAD 


The HDREAD subroutine reads the header data records from a foreign imagery 
tape in one of the •following formats: Universal, Landsat 1 or 2, or 

LARSYS II or III. 

CALLING SEQUENCE : HDRE AD(LUN, 1 FRM, HEADER ,BFSZ. EOF , PC , IE .BUFFER , INS , HDRADR , I SS ) 

CALLING ARGUMENTS: 


Name 

Usage 

Jm. 

Definition 

LUN 

I/O 

I 

Logical unit number assigned to the tape. 

IFRM 

Process 

I 

Format type of the image tape. 

HEADER 

Array 

I 

Buffer. 

BFSZ 

Process 

I 

Buffer size in bytes. 

EOF 

Process 

I 

Flag which, if=l, indicates end of file. 

PC 

Process 

I 

Parity counter. 

IE 

Process 

I 

Error code. 

BUFFER 

Array 

I 

Buffer for data. 

INS 

Process 

I 

Number of scan lines. 

HDRADR 

Process 

I 

Address of the array for WTQIO. 

ISS 

Process 

I 

Starting scan l^ne. 

FILES: 

File 

TAPE 

Usage 

Input Refer to 

Record format 

Landsat or Universal formats. 

COMMON 

BLOCKS : ‘ 

The common block for subroutine HDREAD is /HCOM/. 
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/HCOM/COMMON BLOCK 


Relative 


Name 

word number 

Usage 

liiEi 

Definition 

SS 

1 

Process 

I 

Sample start. 

SE 

2 

Process 

I 

Sample end. 

LS 

3 

Process 

I 

Line start. 

LE 

4 

Procesf 

I 

Line end. 

NRPDS 

5 

Process 

I 

Number of records per data set 
(data set is a scan line). 

NDSPR 

6 

Process 

I 

Number of data sets per. record. 

NCPR 

7 

Process 

I 

Number of channels per record. 

NRPC 

8 

Process 

I s 

Number of records per channel . 

ANCL 

9 

Process 

I 

Length of ancillary block. 

NC 

10 

Process 

I 

Number of channels in the image. 

NS 

11 

Process 

I 

Number of samples per channel 
per scan line. 

NBIT 

12 

Process 

I 

Number of bits in a pixel. 

DO I 

13 

Process 

I 

Data order indicator. 

NCAR 

14 

Process 

I 

Same as NCPR. 

SVD 

15 

Process 

I 

Start of the video within a data set. 

RSIZ 

16 

Process 

I 

Tape record size. 

PSKIP 

17 

Process 

I 

Pixel offset for reading data from tape. 

HSIZ 

18 

Process 

I 

Number of words in the header OAT. 

CALP 

19 

Process 

I 

Used in calculating offset bytes in data 

CERR 

20 

Process 

I 

Catastrophic error. Calling program 
should terminate if error i 0. 

NSPR 

21 

Process 

I 

Number of pixels per record. 
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SUBROUTINE CALLED: 


Reference 


Name 

Definition 

(section) 

HPRC3 

Header process. 

3.1.2 

HDCODE 

Header decode. 

3.1.1 


CALLED BY : The HDREAD subroutine Is called by the driver routine READTP 

(see section 3). 

LOCAL VARIABLES : Local variables for subroutine HDREAD are listed In common 

block /HCOM/ of the HDREAD subroutine In this section. 
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3.1.1 UNIT 1, LEVEL 2 - HOCODE 

The HDCOOE subroutine decodes selected portions of tape header records in one 
of the following formats: Universal, Landsat 1 and 2, or LARSYS II and III. 


CALLING SEQUENCE: 

SUBROUTINE HDCODE(IF,H,IE,LVL) 

CALLING ARGUMENTS: 


Name Usage 

Jm. 

Definition 

IK Process 

I 

Variable indicating format of header; if = 1, 
indicates Universal, if = 2, indicates Landsat, 
if = 3, indicates LARSYS. 

H Array 

A 

Header record exactly as read from tape. 

IE Process 

I 

Variable indicating error code; if = 0, no error; 
if = 1, not 1, 2, or 3. 

LVL Process 

I 

Variable which indicates Landsat level indicator; 
if = 1, level 1, if = 2, level 3. 

FILES: 

File 

Usage 

Record Format 

Tape header 

Input 

Landsat format; refer to the Landsat documentation. 

COMMON BLOCKS: 

Common blocks for the HDCODE subroutine are /COMH/, /LATON/ 


and /LABEL/. For descriptions of common blocks /COMH/ and /LABEL/, refer to 
program READTP (section 3). 


/LATLON/COMMON BLOCK 

Relati ve 


Name 

word number 

Usage 

Im. 

Definition 

LATD 

1 

Process 

A 

Latitude direction. 

LONGD 

2 

Process 

A 

Longitude direction 
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SUBROUTINES CALLED: 


Name 

Definition 

Reference 

(section) 

FRONT 

Delete blanks from string. 

3. 1.1.1 

ETAATE 

EBCDIC to ASCII. ASCII to EBCDIC. 

3. 1.1. 2 

B2I 

Binary to integer. 

3. 1.1. 3 

A2I 

Alpha to integer. 

3. 1.1. 4 

LATLON 

Converts from degrees and minutes to degrees and 
tenths of a degree. 

3. 1.1. 5 

JULIAN 

Converts Gregorian calendar to Julian calendar. 

.3. 1.1. 6 


CALLED BY ; The HDCODE routine is called by subroutine HDREAD (see 
section 3.1). 

LOCAL VARIABLES : Local variables for subroutine HDCODE are listed in common 

block /COMH/. Refer to the driver routine READTP, section 3. 
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3. 1.1.1 Unit 1. Level 3 - FRONT 


The FRONT subroutine left justifies an ASCII character string, removes any 
leading blanks, and blank fills the array after the last character is shifted. 

C ALLING SEQUENCE : SUBROUTINE FRONT(I.N) 

Nam Usage Type Definition 

I .rray A String which will be jxamined for blanks. 

N Process I The length of the string in characters. 

FILE : Not applicable. 

COMMON BLOCKS : Not applicable. 

SUBROUTINES CALLED : Not applicable. 

CALLED BY : Subroutine FRONT is called by the subroutine HDCODE (see 

section 3.1.1). 

LOCAL VARIABLES : 

Name Definition 

K Counter for the number of characters examined. 

FILES : Not applicable. 

COMMON BLOCKS: Mot applicable. 
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3. 1.1. 2 Unit 1. Level 3 - ETAATE 

The ETAATE subroutine converts an EBCDIC character to an ASCII character and 
an ASCII character to an EBCDIC character one byte at a time. 

CALLING SEQUENCE : SUBROUTINE ETAB(H.ARG.NCHAR) 

Name Usage Type Definition 

H Array A String input of EBCDIC character; may be one character. 

AR6 Array A String output of ASCII character; may be one character. 

NCHAR Process I Number of characters to be converted., 

SUBROUTINES CALLED : Not applicable. 

CALLED BY ; The ETAATE subroutine is called by 'subroutine HDCODE (see 
section 3.1.1). 

LOCAL VARIABLES: Not applicable. 
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2. 1.1.3 Unit 1. Level 3 - B2I 


The B2I subroutine converts binary data to integers. 

CALLING SEQUENCE : Function B2I(H,Bl,B2) 

CALLING ARGUMENTS : 

Name Usage Type Definition 

H Array binary Array containing the string on which to operate. 

B1 Process I First character in string on which to operate. 

B2 Process I Last character in string on wh :h to -operate. 

FILES : Not applicable. 

COMMON BLOCKS : Not applicable. 

SUBROUTINES CALLED : Not applicable. 

CALLED BY : Ihe B2I subroutine is called by subroutine HDCODE (see 

section 3.1.1). 

LOCAL VARIABLES : 

Name Definition 

C The scratch pad word containing the byte on which to operate. 
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3. 1.1.4 Unit 1, Level 3 - A2I 

The A2I subroutine converts ASCII data to Integer data. 

CALLING SEQUE NCE: Function A2I(H,B1,B2) 

CALLING ARGUMENTS ; 

Name Usage Type De finition 

H Array Binary Array containing the string on which to operate. 

B1 Process I First character on which to operate in string. 

B2 Process I Last character on which to operate in string. 

FILE : Not applicable. 

COMMON BLOCKS ; Not applicable. 

SUBROUTINE CALLED ; Not applicable. 

CALLED BY: The A2I subroutine is called by subroutine HDCODE (see 
section 3.1.1). 

LOCAL VARIABLES ; 

Name Definition 

C The scratch pad word containing the byte on which to operate. 
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3. 1.1. 5 Unit 1. Level 3 - LATLON 

The LATLON subroutine converts the latitude and longitude from degrees and 
minutes to degrees and tenths of a degree. 

CALLING SEQUENCE ; REAL FUNCTION LATL0N(H,B1 ,B2) 

C ALLING ARGUMENTS : 

Name Usage Type Definition 

H Array Byte Array containing the latitude and longitude. 

B1 Process I First character in the array H desired. 

B2 Process I Last character in the array H desired. 

FILES ; Not applicable. 

COUMC r; B LXKS ; Not applicable. 

SUBROUTINES CALLED ; Not applicable. 

CALLED BY : The LATLON subroutine is called l*y the HOCOOE subroutine (see 

section 3.1.1). 

LXAL VARIABLES ; 

Name Definition 

LM Decimal that places calculation for numeric packing. 

LL Decimal that places calculation for numeric packing. 
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3. 1.1.6 Unit 1. Level 3 - JULIAN 


The JULIAN subroutine converts the Gregorian calendar date to the Julian 
calendar integer. 

CALLING SEQUENCE ; FUNCTION JULI AN(M,D,Y) 

CALLING ARGUMENTS ; 

Name Usage Type Definition 

M Process I Month value. 

0 Process I Day value. 

t 

Y Process I Year value. 

FILES ; Not applicable. 

COMMON BLOCKS ; Not applicable. 

SUBROUTINE CALLED ; Not applicable. 

CALLED BY ; The JULIAN subroutine is called by the HDCODE subroutine (see 
section 3.1.1). 

LXAL VARIABLES ; 

Name Definition 

Julian The Julian date in integer form. 
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3.1.2 UNIT 1, LEVEL 2 - HPROS 

The HPROS subroutine builds a common block of header data for subsequent 
program use. 

CALLING SEQUENCE : HPR0S(LUN1, HEADER, FMT, EOF, PRTY.IPRM) 

CALLING ARGUMENTS : 

Name Usage Type 
LUNl Input I 
HEADER Array A, I 

FMT Process I 

EOF Flag I 

PRTY Process I 
IPRM Process I 

FILES : 

File Usage Record Format 
Tape Input Input data tape. 

COMMON BLOCKS : The common block for the HPROS subroutine is /HCOM/. Refer to 

subroutine HDREAD (section 3.1) for a description of this common block. 

SUBROUTINES CALLED : Not applicable. 

CALLED BY : The HPROS subroutine is called by subroutine HDREAD (see 

section 3.1). 

LOCAL VARIABLES : Local variables for subroutine HPROS are listed in common 

block /HCOM/. Refer to the HDREAD subroutine, section 3.1. 


Definition 

Logical unit number assigned to the tape input. 

Buffer array used to store the header. Data from 
the tape header. 

Format number of the headers, 

1 = Universal, 

2 = Landsat, 

3 = LARSYS. 

End-of-file indicator; 1 = EOF. 

Parity count. 

Parameter value. 
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3.1.3 UNIT 1, LEVEL 2 - IBYTE 

The IBYTE subroutine converts integers to byte form, one byte at a time. 
CALLING SEQUENCE: FUNCTION IBYTE(START, STRING) 

CALLING ARGUMENTS : 

Name Usage Type Definition 

START Process I Starting position on string. 

STRING Array I Array containing the values which require change. 

FILES ; Not applicable. 

COMMON BLOCKS ; Not applicable. 

SUBROUTINE CALLED ; Not applicable. 

CALLED BY ; The IBYTE subroutine is called by subroutine HDREAD (see 
section 3.1). 

LOCAL VARIABLES ; Not applicable. 
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3.2 UNIT 2 - INTOCT 

The INTOCT subroutine converts the Julian date to OCTAL and converts the 
version number to OCTAL. 

CALLING SEQUENCE ; Subroutine INTOCT(INTEGR, OCTAL) 
calling ARGUMENTS: 


Name 

Usage 

Jige 

Definition 

INTEGR 

Process 

I Integer to be 

converted to OCTAL. 

OCTAL 

Process 

A Integer value 

following conversion to OCTAL. 

FILES: 

Not applicable. 



COMMON BLOCKS : Not applicable. 

SUBROUTINES CALLED: Not applicable. 

CALLED BY : The INTOCT subroutine is called by the driver routine REAOTP (see 

section 3). 

LOCAL VARIABLES : 

Name Definition 

VAL Temporary storage word for calculation. 

X Remainder of a division. 

Z Counter for the number of iterations. 
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3.3 UNIT 3 - HDRFLN 


The HDRFLN subroutine writes a header record in IMDACS consisting of a file 
header record and a Universal imagery header record (SPU format header). The 
file header and the format header record consist of 1536 words each. 

CALLING SEQUENCE: SUBROUTINE HDRFLN(LUN2,FLNIMH,UIFHR,SUNEL,FLAT,FL0NG) 

CALLING ARGUMENTS: 


Name 

Usage 

Im. 

Definition 

LUN2 

In/Out 

I 

Logical unit, number assigned to a file on the disk. 

FLNIMH 

Array 

A 

Array containing the file name. 

UIFHR 

Array 

A, I 

Array containing the header data. 

SUNEL 

Process 

I 

Sun elevation angle. 

FLAT 

Process 

R 

Format center latitude. 

FLONG 

Process 

R 

Format center longitude. 

FILES: 

File 

Usage 

Record Format 


??????. IMH Output Header file in IMDACS format. 
COMMON BLOCKS : Not applicable. 

SUBROUTINE CALLED : 

Reference 

Name Definition (section) 

UNVRSL Universal 3.3.1 


CALLED BY : The HDRFLN subroutine is called by the driver routine READTP (see 

section 3). 

LOCAL VARIABLES : Not applicable. 
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3.3.1 UNIT 3, LEVEL 2 - UNVRSL 

The UNVRSL subroutine builds an SPU universal header record for imagery files. 
CALLING SEQUENCE: SUBROUTINE UNVRSL(LUN2.SUNEL.FLAT.FL0NG) 


CALLING ARGUMENTS: 


Name Usage 
LUN2 I/O 
SUNEL Process 
FLAT Process 
FLONG Process 


Type Definition 

I Logical unit number assigned to a file on the disk, 
I Sun elevation angle. 

R Frame center latitude. 

R Frame center longitude. 


FILES : Not applicable. 

COMMON BLOCKS : Common block for the UNVRSL routine are /LABEL/ and /LATLON/. 
For a description of /LABEL/, refer to the driver routine REAOTP (section 3). 
For a description of /LATLON/, refer to subroutine HDCODE (section 3.1.1). 


SUBROUTINES CALLED : Not applicable. 

CALLED BY : The UNVRSL subroutine is called by subroutine HDRFLN (see 

section 3.3). 

LOCAL VARIABLES : Local variables for the UNVRSL routine are listed under 

variables in the IMDACS documentation for the header file. 
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3.4 UNIT 4 - HOT 


The HOT subroutine extracts patch images from the full-frame image, with 
processing taking place one strip of CCT at a time. It further builds a 
16-by-16 patch frame image consisting of 32-sample by 32-line pixel patches. 
Each pixel in a patch is composed of bands 4 through 7. This subroutine uses 
HOT full-frame images. 

CALLING SEQUENCE : Subroutine HDT(LUN1,LUN2,BUF,BFS2,FLNIMD, EOF, IE, BUFFER) 

CALLING ARGUMENTS: 


Name 

Usage Type 

Definition 

LUNl 

In/out 

I 

Logical Unit nun^ber assigned to the tape. 

LUN2 

In/out 

I 

Logical unit number assigned to the disk. 

BUF 

Array 

I 

Array containing the data. 

BFSZ 

Array 

I 

Size of the data array. 

FLNIMD 

Array 

I 

Array containing the name of the output file. 

EOF 

Array 

I 

End of file. 

IE 

Array 

I 

Error code. 

BUFFER 

Array 

I. A 

Used as a scratch pad for pixel and header data. 

FILES: 





File Usage Record Format 

??????. IMD Output file Consult with the IMDACS system documentation. 


COMMON BLOCKS : The common blocks for the HOT routine are /COMH/ and 

/LABEL/. For a description of each, refer to routine READTP (section 3). 
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SUBROUTINES CALLED: 


Reference 

Name Definition (section) 

LBORDR Line border 3.4.1 

PBORDR Pixel border 3.4.2 

LDCODE Line decode 3.4.3 

HIST Histogram 3.4.5 

CALLED BV : The HOT subroutine is called by the driver routine READTP (see 

section 3). 

LOCAL VARIABLES ; 

Name Definition 

BLOCKN Block number. 

ENDCCT Parameter for reading to the end of the file. 

LCOUNT Line counter for the patch image. 

LEND Ending line within the full frame. 

LSTART Starting line within the full frame. 

PATCH Current patch number. 

PIXELS Pixel location in PIXBUF. 

PTCH Flag indicating that the patch counter should be decremented. 

SAVEP Last patch processed on previous CCT. 

STARTP Starting pixel location in patch extraction. 

STOPP Last pixel location in patch extraction has been readjusted. 
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3.4.1 UNIT 4, LEVEL 2 - LBORDR 

The LBORDR subroutine determines when imagery data desired from full frame is 
not available; patch lines are outside of the full frame. 


CALLING SEQUENCE: 

SUBROUTINE LBORDR(STARTL, NEXTL, LS, LE , LEND, LSTART) 

CALLING 

ARGUMENTS 



Name 

Usage 

Type 

Definition 

STARTL 

Process 

I 

Calculated starting line from calling routine. 

NEXTL 

Flag 

L 

Flag which causes the line not to be processed. 

LS 

Process 

I 

Starting line of full frame. 

LE 

Process 

I 

Last line of the full frame. 

LEND 

Process 

I 

Last line within the full frame. 

LSTART 

Process 

I 

First line within the full frame. 

FILES: 

Not applicable, 

• 


COMMON BLOCKS : Not applicable. 

SUBROUTINES CALLED : Not applicable. 

CALLED BY ; The LBORDER subroutine is called by subroutines HDT and LANSAT 
(see sections 3.4 and 3.5, respectively). 

LOCAL VARIABLES : 

Name Definition 
LEND Line end. 

LSTART Line start. 

NEXTL Flag which, if = TRUE, instructs program to process next full 
frame scan line. 
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3.4.2 UNIT 4, LEVEL 2 - PBOROR 

The PBOROR subroutine determines when the imagery data desired from the full 
frame is not available; patch pixels are outside of the full frame. 

CALLING SEQUENCE : SUBROUTINE PBORDR(STARTP.STOPP, STOP, NEXTP.SE, SAT) 

CALLING ARGUMENTS : 

Name Usage Type 
STARTP Process I 

STOPP Process I 

STOP Process I 

NEXTP Flag L 

SE Process I 

ST Process I 

FILES ; Not applicabl 

COMMON BLOCKS : Not applicable. 

SUBROUTINES CALLED ; Not applicable. 

CALLED BY : The PBORDR subroutine is called by subroutines HDT and LANSAT (see 
section 3.4 and 3.5, respectively). 

LOCAL VARIABLES : 

Name Definition 

NEXTP Flag which, if = TRUE, instructs program to process next patch. 


Definition 

Starting pixel number provided by the calling 
routine; altered if necessary for each patch. 

Stop pixel for each patch; altered pdrtion of 
patch is on the next CCT. 

Stop pixel for each patch. 

Flag indicating if the patch should be processed. 
Last pixel in the scan line. 

First pixel in the scan line. 
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3,4.3 UNIT 4, LEVEL 2 - LDCODE 

The LDCODE subroutine decodes image data from a foreign imagery tape in one of 
the following tape formats: Universal, Landsat 1 and 2, or LARSYS II and III. 

CALLING SEQUENCE : SUBROUTINE LDCOOE(LUN, IF, BUFFER, BFSZ,DLIN.RCHAN,LSS, 

LSE,BUF,EOF,PC,IE,INIT) 

CALLING ARGUMENTS: 


Name 

Usage 

T^ 

Definition 

LUNl 

In/out 

I 

Logical unit number assigned to the tape. 

IF 

Process 

I 

Format type of the tape. 

BUFFER 

Array 

Byte 

Array containing the scan line data and scratch 
pad buffer. 

BFSZ 

Process 

I 

Size of array buffer. 

DLIN 

Process 

I 

Requested line number; starts with line 1. 

RCHAN 

Process 

I 

Requested channel number, starts with 1. 

LSS 

Process 

I 

Pixel with which to start data move; assumes 
first pixel on tape. 

LSE 

Process 

I 

Last pixel desired. 

BUF 

Array 

Byte 

Decoded pixels LSS through LSE for calling program 

EOF 

Process 

I 

Flag which, if = 1, indicates end of file; if = 0, 
indicates more data is to come. 

PC 

Process 

I 

Parity check counter. 

IE 

Process 

I 

Error code. 

INIT 

Process 

I 

Initialization flag which if = 0, indicates first 
pass; if = 1, not first pass. 

FILES: 

Name 

Usage 

Record Format 


Tape Input Data tape being read. 
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COMMON BLOCKS : Common blocks for the LDCODE subroutine are /HCOM/ and 

/LABEL/. For a description of each, refer to sections 3.1 and 3, 

>"especti vely . 

SUBROUTINES CALLED : Not applicable. 

CALLED BY : The LDCODE subroutine is called by subroutines HOT and LAN5AT (see 

sections 3.4 and 3.5, respectively. 

LOCAL VARIABLES ; Local variables for the LDCODE subroutine are listed in 
common block /HCOM/ in the subroutine HDREAD (section 3.1). 
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3.4.4 UNIT 4, LEVEL 2 - PARTP 

The PARTP subroutine adjusts the parameters when a patch is divided between 
two CCr strips. 

CALLING SEqUENCE ; Subroutine PARTP(PTCH.STARTP.STOPP,STOP.SEqNO.SE,ST,PIXELS, 


I, PATCH) 



CALLING 

ARGUMENTS: 


Name 

Usage 


Definition 

PTCH 

Flag 

L 

Flag which determines if the patch counter needs 
to be decremented. 

STARTP 

Process 

I 

Starting pixel number provided by the calling 
routine; altered if necessary for each patch. 

STOPP 

Process 

I 

Stop pixel for each patch; altered portion of patch 
is on the next CCT. 

STOP 

Process 

I 

Stop pixel for each patch. 

SEQNO 

Process 

I 

Sequence number for each CCT. 

SE 

Process 

I 

Last pixel in the scan line for the current CCT. 

ST 

Process 

I 

First pixel in the scan line for the current CCT. 

PIXELS 

Array 

I 

Array containing the desired patch pixels. 

I . 

Process 

I 

Band number'' ' through 4. 

PATCH 

Process 

I 

Patch number for resetting counter if necessary. 

FILES: 

Not applicabl 

e. 

COMMON BLOCKS: 

Not applicable. 

SUBROUTINES CALLED: 

Not applicable. 


CALLED BY ; The PARTP subroutine is called by subroutines HDT and LANSAT (see 
section 3.4 and 3.5, respectively). 
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LOCAL VARIABLES: 

Name Definition 

CCTNO The CCT strip number that Is last stored. 

PIECE Array containing the parameter values to be used for the remainder 
of the rightmost patch located on the next CCT strip. 
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3.4.5 UNIT 4, LEVEL 2 - HIST 

The HIST subroutine builds histograms for each channel. 


CALLING SEQUENCE ; Subroutine HIST(BUF,STARTP,STOPP.RCHAN.ST) 

CALLING ARGUMENTS : 

Name Usage Type Definition 

BUF Array Byte Array containing the scan line from the full frame. 

STARTP Process I Starting pixel number provided by the calling 

routine. 

« 

STOPP Process I Stop pixel for each patch. 

RCHAN Process I Current channel number. 

ST Process I Starting pixel number for the cu*^rent CCT. 

FILES : Not applicable. 

COMMON BLOCKS : Common blocks for the HIST subroutine are /HIST/ and 

/LABEL/. For a discription of each, refer to routine READTP (section 3). 

SUBROUTINES CALLED : Not applicable. 

CALLED BY ; The HIST subroutine is called by subroutines HDT and LANSAT (see 
sections 3.4 and 3.5). 

LOCAL VARIABLES : Not applicable. 
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3.5 UNIT 5 - LANSAT 


The LANSAT subroutine extracts patch images from the ful 1-frame Image, with 
processing taking place one strip of CCT at a time. It further builds a 


16-by-16 patch frame Image consisting of 32-sample by 32-1 Ine pixel patches. 
Each pixel In a patch Is comprised of bands 4 through 7. The subroutine uses 
Landsat 1 and 2 full-frame size Images. 

CALLING SEQUENCE: SUBROUTINE LANSAT(LUin,LUN2,BUF,BFSZ,FLNIMD, EOF. IE, BUFFER, 

IFRM,SS,INIT) 

CALLING 

ARGUMENTS: 


Nama 

Usage 

Type 

Definition 

LUNl 

I/O 

I 

Logical unit number assigned to the tape. 

LUN2 

I/O 

I 

Logical unit number assigned to the disk file. 

BUF 

Array 

I 

Buffer containing the pixel data per scan line. 

BFSZ 

Process 

I 

Buffer size to be used. 

FLNIMD 

Array 

I 

Array containing the file name of the output file. 

EOF 

Process 

I 

End of file. 

IE 

Process 

I 

Error code. 

BUFFER 

Array 

I 

Scratch pad array containing header or image data. 

IFRM 

Process 

I 

Format type of the tape. 

SS 

Process 

I 

Start pixel number. 

INIT 

Process 

I 

Initializing flag. 


For further Information on the LANSAT subroutine, consult the documentation 
written for subroutine HOT (see section 3.4). 
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3.6 UNIT 6 - HDRHIS 


The HDRHIS subroutine writes the histograms for each band (4-7) to the IMDACS 
header record file. 

CA LLIN6 SEQUENCE ! SUBROUTINE HDRHIS(LUN2,FLNIMH,UIFHR) 

CALLING ARGUMENTS : 

Name Usage Type Definition 

LUN2 I/O I Logical unit number assigned to a file on the disk. 

FLNIMH Array A Array containing the file name. 

UIFHR Array A, I Array containing the header data. 

FILES ; 

File Usage Record Format 

??????. IMH Header file Consult IMDACS documentation for information on 

format . 

COMMON BLOCKS : 

Rel ati ve 


Name 

word number 

Usage 


Definition 


H4 

512 

Array 

R 

Histogram for band 

4. 

H5 

1025 

Array 

R 

Histogram for band 

5. 

H6 

1537 

Array 

R 

Histogram for band 

6. 

H7 

2049 

Array 

R 

Histogram for band 

7. 


SUBROUTINE CALLED ; Not applicable. 

CALLED BY ; The HDRHIS subroutine ‘ called by the driver routine READTP (see 
section 3). 

LOCAL VARIABLES : Not applicable. 
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4. SUBSYSTEM 2 - SKIP 


The SKIP routine is the driver routin'!' for the skipped image generation 
program. 

CALLING SEQUENCE : Not applicable. 

CALLING ARGUMENTS : Not applicable. 

FILES ; 

Fi le Usage Record format 

Input Cards Card 1, columns 1-6; File name. 

Card 2, columns 1-2; Line start number. 

Card 3, columns 1-2; Pixel start number. 

COMMON BLOCKS ; Common blocks for the SKIP routine are COMH, LABEL, and 
HIST. Refer to driver routine READTP (section 3) for a description of each. 

SUBROUTINE CALLED; 


Reference 


Name 

Definition 

(section 

ASNLUN 

Assign logical 

unit number. 

System 

CL OSS 

Close. 


System 

6ETADR 

Get address. 


System 

HDREAD 

Header record 

read. 

3 

HDRFLN 

Header record 

file name. 

3.3 

HDRHIS 

Header record 

histogram. 

3.6 

OPENS 

Open System. 
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Name Definition 

QIO Queue I/O. 

SKPRD Skip read. 

CALLED BY: Not applicable. 


Reference 

(section) 

System 

2 
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4.1 UNIT 1, LEVEL 1 - SKPRD 

The SKPRD subroutine controls the reading and writing of the data to be 
selected from the full Landsat frame for the skipped image. 

CALLING SEQUENCE: CALL SKPRD(LUN1, LUN2,BUF,BFSZ, EOF, IE, BUFFER. IFRM.HDRADR, 


ICCT.LINST.IPIXST) 


CALLING 

ARGUMENTS: 


Name 

Usage 


Definition 

uUNl 

I/O 

I 

Logical unit for input tape. 

LUN2 

I/l 

I 

Logical unit for output data files. 

BUF 

Array 

Byte 

Buffer for pixel data. 

BFSZ 

Process 

I 

Size of buffer BUF. 

EOF 

Process 

I 

End of file indicator for tape read; if = 1, end 
file is found. 

IE 

Process 

I 

Error indicator for tape read; if = 0, no error. 

BUFFER 

Array 

Byte 

Buffer for reading data from tape to be used by 
LDCODE. 

IFRM 

Process 

I 

Format of input tape; IFRM = 2 for LANSAT. 

HDRADR 

Process 

I 

Address of buffer. 

ICCT 

Process 

I 

The CCT strip number that is being processed. 

LINST 

Process 

I 

Start line number for building image. 

IPIXST 

Process 

I 

Pixel start number for building image. 

FILES: 

Not applicable. 

9 

COMMON 

BLOCKS: 

The common block for the SKIP routine is /HIST/. Refer 


driver routine READTP (section 3) for a description of this common block. 
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SUBROUTINE CALLED : 

Reference 

Name Definition (section) 

LDCOOE Line decode 4.1.1 
SKPWRT Skip write 4.1.2 

CALLED BY ; The SKPRD subroutine is called by the driver routine SKIP (see 
section 4). 

LOCAL VARIABLES : 

Name Definition 

I LINE Line number to be read from tape. 

LINE Line number to write to disk file. 

LSS Pixel start for CCT strip; always equal 1. 

LSE Pixel stop for CCT strip; equals number of pixels per scan on 
this strip. 

IPIXL First pixel to remove from this strip. 

IPADO Skip factor to add for pixel skipping; skip factors are 6, 6, 6, 
and 7. 

LADD Skip factor to add for line skipping; and skip factors are 4 and 5. 

ISTRT Index into array IPADD; assures that the skip factor is not changed 
across CCT strips. 
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4.1.2 UNIT 1. LEVEL 2 - SKPWRT 

The SKPWRT routine writes the selected pixels for the skipped image to the 
output disk file. 

CALLING SEQUENCE: CALL SKPWRT{LUN2,BUF.INDX,ICHAN,ICCT,LINE) 


CALLING ARGUMENTS: 


Name Usage Type Definition 

LUN2 I/O I Logical unit number for output data file. 

BUE Array Byte Buffer containing pixel data to output. 

INDX Process I Count of number of entries to process from BUF 

ICHAN Process I Channel number for data in BUF. 

ICCT Process I The CCT strip number to process. 

LINE Process I Line number to write data in BUF. 


FILES : 

File Usage Record format 

FILENM.IMD Image data IMDACS data file format. 

COMMON BLOCKS : Not applicable. 

SUBROUTINES CALLED : 

Reference 

Name Definition (section) 

SKPHST Skip histogram 4. 1.2.1 

CALLED BY : The SKPWRT subroutine is called by routine SKPPD (see 

section 4.1). 
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LOCAL VARIABLES: 


Name 

ICNT 

OBUF 


Definition 


Count of number of pixels written to scan line for previous CCT 
strips; if ICCT » 1, then ICNT = 0. 

Output buffer to write to the file. 
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4.1.2. 1 UNIT 1. LEVEL 3 - SKPHST 


The SKPHST subroutine computes the histogram for the skip image program. 


CALLING SEQUENCE : CALL SKPHST(OBUF) 


CALLING ARGUMENTS: 


Name Usage Type 
OBUF Array Byte 


Definition 

Buffer containing pixel data to write skipped 
image data file. 


FILE: Not applicable. 


COMMON BLOCKS : The common block for subroutine SKPHST is HIST. Refer to the 

driver routine READTP (section 3) for a description of this common block. 


SUBROUTINES CALLED : Not applicable. 

CALLED BY : The SKPHST subroutine is called by subroutine SKPWRT (see 

section 4.1.2). 


LOCAL VARIABLES : Not applicable. 


5. SUBSYSTEM 3 - PATCHP 


The PATCHP program is the driver for the patch processor and produces analyst 
aid statistics for the 32-sample by 32-line pixel patches in a 512-sample by 
512-line pixel image file. 

CALLING SEQUENCE ; Not applicable. 

CALLING ARGUMENTS : Not applicable. 

FILES ; 

Fi le 
INPUT 

FILENM.IMH 
FILENM.IMD 
FILENM.DAT 

HEADER RECORD 


Col umns 

Format 

Contents 


1 

IX 

Blank. 


2-11 

lOAl 

Scene identification. 


12-13 

12 

Sun angle. 


14 

11 

Satellite number. 


15 

A1 

Latitude direction north 

or south. 

16-19 

14 

Latitude degrees. 


20-21 

12 

Latitude minutes. 


22 

A1 

Longitude direction east 

or west. 


Usage Record format 

CARDS Card 1, columns 1-30; File name for header record file = 
FILENM.IMH. 

Card 2, columns 1-30; File name for image data file = 
FILENM.IMD. 

Input IMDACS image header record file. 

Input IMDACS image data file. 

Output Analyst aid statistics file. 



Columns Format 


Contents 


23-26 

14 

Longitude degrees 

27-28 

12 

Longitude minutes 

29-30 

2X 

Blank. 

31-36 

A6 

File name. 


DATA RECORD - 256 RECORDS, ONE RECORD PER PATCH 
Columns Format Contents 


1 

IX 

B1 ank . 

2-4 

13 

Patch number. 

5-28 

4F6.2 

Channel means. 

29-32 

14 

Garbled pixel count 

33-36 

14 

Cloud pixel count. 

37-40 

14 

Water pixel count. 

41-46 

F6.3 

GAMA. 

47-49 

13 

Soil line number. 

50-53 

F4.1 

Green number. 

54-57 

F4.1 

AVI. 

58-61 

F4.1 

PVI. 

62-66 

F5.3 

TVI7. 

67-71 

F5.3 

LAI. 

72-74 

13 

GIN5. 

75-77 

13 

GIN15. 

78-80 

13 

GIN30. 
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COMMON BLOCK: The common block for program PATCHP is PATCH. 


PATCH 


Relative 


Name 

word number 

Usage 

Im. 

Definition 

CH 

1-4096 

Array 

I, A 

Array dimensioned (1024,4) containing 
1024 pixel values for each of 4 chan- 
nels belonging to a 32-by-32 patch. 

KH 

4096-9192 

Array 

I, A 

Array dimensioned (1024,4) containing 
1024 Kauth-transformed pixel values for 
each of 4 channels belonging to a 
32-by-32 patch. 

FLAG 

9193-9704 

Array 

Byte 

Array dimensioned (1024) to contain a 

flag value for each pixel after 

screening algorithm has been performed 

on a patch image; 

value of 1 = good pixel, 

value of 10-29 = water pixel, 

value of 30-39 = cloud pixel, 

value of 40-49 = garbled pixel. 

X 

9705-9720 

Array 

Real*4 

Four-member array containing mean pixel 
value for each good pixel on each of 
4 channels per patch. 

GAR 

9721-9722 

Process 

I*A 

Count of garbled pixels in patch. 

CLOUD 

9723-9724 

Process 

I*A 

Count of cloud pixels in patch. 

WATER 

9725-9726 

Process 

I*A 

Count of water pixels in patch. 

GAM 

9727-9730 

Process 

Real 

Integer form of GAMA. 

VEG 

9731-9778 

Array 

Real 

Array dimensioned (12) to save vege- 
tative indices. 

GIN 

9779-9784 

Array 

I*A 

Array dimensioned (3) to contain values 
for GIN5, GIN15, and GIN30. 

I AG 

9785-9786 

Process 

I*A 

Sun angle. 

SAT 

9787-9788 

Process 

I*A 

Satellite number. 

SLN 

9789-9790 

Pr'ocess 

I*A 

Soil line number. 

GAMA 

9791-9794 

Process 

Real 

X-STAR haze parameter. 
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SUBf^OUTINES CALLED: 


Reference 


Name 

Definition 

(section) 

CL0S$ 

Close 

System 

GAMMA 

X-STAJl depth correction 

5.3 

HAZCOR 


5.4 

HEADER 

Header process 

5.1 

OPENS 

Open 

System 

PATCRD 

Patch read 

5.2 

PATRPT 

Patch report 

5.6 

VIN 

Vegetation index 

5.5 

CALLED 

BY : Not applicable. 


LOCAL ' 

VARIABLES: 


Name 

Definition 


IPATC), 

Patch number being proce 

issed (1-266) 
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5.1 UNIT 1. LEVEL 1 - HEADER 

The HEADER subroutine reads the universal header record from the image header 
file. It decodes some information from it, and writes it for the header 
record of the output analyst aid statistics file. 

CALLING SEQUENCE ; CALL HEADER( IMHNAM,FILANAM,BUF,IAG,SAT,IE) 

CALLING ARGUMENTS; 


Name 

Usage 

Im. 

Definition 

IMHNAM 

Array 

Byte 

Name of IMDACS image header file. 

FILNAM 

Array 

Byte 

Name of analyst aid statistics output file. 

BUF 

Array 

Byte 

Array to read universal header record into 
record 2 of header file. 

I AG 

Process 

1*2 

Sun angle. 

SAT 

Process 

1*2 

Satellite number. 

IE 

Process 

1*2 

Error indicator; if error is found, set to 1. 


FILES : See program PATCHP (section 5) for files information. 

COMMON BLOCKS : Not applicable. 

SUBROUTINE CALLED : 

Name Definition Level 
CLOS$ Close System 

CALLED BY : The HEADER subroutine is called by the program PATCHP (see 
section 5). 

LOCAL VARIABLES : Not applicable. 
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5.2 UNIT 1, LEVEL 1 - PATCHRD 

The PATCRD subroutine reads four channels of imagery data for one patch into 
an integer *4 buffer. 

C ALLING SEQUENCE : CALL PATCRD(LUN, IPATCH,BUF ) 

CALLING ARGUMENTS: 


Name 

Usage 

Jm. 

Definition 


LUN 

I/O 

1*2 

Logical unit number for IMDACS 

image data file. 

IPATCH 

Process 

1*2 

Patch number for which data is 

desired. 

BUF 

Array 

1*2 

(1024,4) dimensioned array for 

imagery data. 


FILES : The PATCRD routine reads the IMDACS image dai:a file (for description, 

see PATCHP, section 5). 

COMMON BLOCKS : Not applicable. 

SUBROUTINES CALLED: Not applicable. 

CALLED BY : The PATCRD subroutine is called by the PATCHP program (see 

section 5). 

LOCAL VARIABLES : 

Name Definition 

ILIN Line number within image file where requested patch data starts. 

IPIX Pixel number within ILIN where requested patch data starts. 
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5.3 UNIT 1, LEVEL 1 - GAMMA 

The subroutine GAMMA computes the X-STAR optical depth parameter. 
CALLING SEQUENCE; CALL GAMMA 
CALLING ARGUMENTS ; Not Applicable. 

FILES : Not applicable. 

COMMON BLOCKS : The common block for the GAMMA subroutine is PATCH (see 

section 5). 

SUBROUTINE CALLED: 


Reference 


Name 

Definition 

(section 

CAL IB 

Cal ibrate mean 

5.3.1 

MEAN 

Calculate mean 

5.3,2 

SCRAL 

Screen 

5.3.3 

$SIN 

Sine function 

System 


CALLED BY : The GAMMA subroutits^i is called by the PATCHP program (see 

section 5). 

LOCAL VARIABLES. : 

Name Definition 

CR Sun-angle correction factor; corrects to 52 degrees. 

2 Sun-angle correction applied to each pixel value. 
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5.3.1 UNIT 1, LEVEL 2 - CALIB 

The CALIB subroutine applies calibration constants to the pixel data. The 
calibration constants are dependent on the satellite number in the following 
way, where LK = the satellite correction constant from which to choose: 

LK = -1 for postlaunch Landsat 2 

LK = 0 for no correction 

LK = 1 for Landsat 1 calibration 

LK = 2 for prelaunch Landsat 2 (no correction) 

LK = 3 for LEC Landsat 3 calibration 

LK = 4 for ERIM Landsat 3 calibration 

CALLING SEQUENCE: CALL CALIB(LK,N) 

CALLING ARGUMENTS: 


Name 

Usage 

Im. 



Definition 

LK 

Process 

1*2 

Indicates 

satellite correction constant to choose. 

N 

Process 

1*2 

Number 

equals 

of possible pixel values to be calibrated 
1024. 


FILES : Not applicable. 

COMMON BLOCKS : The common block for the CALIB subroutine is PATCH (see 

section 5). 

SUBROUTINE CALLED : Not applicable. 

CALLED BY : The CALIB subroutine is called by the subroutine GAMMA (see 

section 5.3). 
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LOCAL VARIABLES; 


Name 
CALF AC 


CALBAS 


Definition 

Calibration factor array (channel, LK), as followi 


1.040 

1.000 

1.090 

0.820 

1.275 

1.141 

1.098 

0.948 

1.161 

1.230 

1.246 

1.062 

1.137 

1.173 

1.247 

1.126 


Calibration bias array (channel, LK), as follows: 


- 5.790 

1.190 

- 2.910 

3.010 

- 1.445 

- 2.712 

- 2.950 

0.446 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 
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5.3.2 UNIT 1, LEVEL 2 - MEAN 

The MEAN subroutine computes the mean and standard deviation for a requested 
channel in the patch image. Only good pixels will be used in the 
computations. 

CALLING SEQUENCE ; CALL MEAN(ICHAN.N,NN,STDEV) 

CALLING ARGUME NTS; 


Name 

Usage 

Tj^ 

Definition 

I CHAN 

Process 

1*2 

Channel for which calculations are requested. 

N 

Process 

1*2 

Number of possible pixel values to use in computa- 
tion; equals 1024. 

NN 

Process 

1*2 

Count of number of pixel values actually used in 
computations. 

STDEV 

Process 

1*2 

Standard deviation for channel requested. 


FILES ; Not applicable. 

COMMON BLOCKS : The common block for subroutine MEAN is PATCH (see section 5). 

SUBROUTINE CALLED ; 

Reference 

Name Definition (section) 

$S0RT Square root System 

CALLED BY : The MEAN subroutine is called by GAMMA subroutine GAMMA (see 

section 5.3). 
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LOCAL VARIABLES ; 

Name Definition 

S Sum of good pixel values. 

SS Sum of square of good pixel values. 
VARI Variance on requested channel. 
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5.3.3 UNIT 1, LEVEL 2 - SCRALL 

The subroutine SCRALL uses the Kauth transformation and screens the pixel data 
and calculates GINS, GIN15, and GIN30. 

CALLING SEQUENCE ; CALL SCR ALL (N) 

CALLING ARGUMENTS ; 

Name Usage Type Definition 

N Process I'^2 Number of pixels to screen equals 1024. 

FILES ; Not applicable. 

COMMON BLOCKS : The common block for subroutine SCRALL is PATCH (see 

section 5). 

SUBROUTINE CALLED : 

Name Definition Level 
$SIN Sine computation System 

CALLED BY ; Subroutine SCRALL is called by GAMMA (see section 5.3 for 
description) . 

LOCAL VARIABLES : 

Name Definition 

RR Kauth transformation. 

2 Screening algorithm test vectors. 

CC Screening algorithm thresholds. 

iI,JJ Screening algorithm pointers. 

LL Screening algorithm output values. 
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5.4 UNIT 1, LEVEL 1 - HAZCOR 

The HAZCOR subroutine applies the XSTAR haze correction that is given to GAMMA 
and the Sun angle. 

CALLING SEQUENCE ; CALL HAZCOR 

CALLING ARGUMENT: Not applicable. 

FILES : Not applicable. 

COMMON BLOCKS : The common block for the HAZCOR subroutine is PATCH {see 

section 5 for a description). 

SUBROUTINES CALLED ; 

Reference 

Name Definition (section) 

KAUTH Kauth transform 5.4.1 

$EXP Exponent System 

$SIN Sine function System 

CALLED BY : The HAZCOR subroutine is called by the driver routine PATCHP (see 

section 5). 

LOCAL VARIABLES : 

Name Definition 

COR Sun-angle correction factor. 

A Gain for correction. 

B Bias for correction. 
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5.4.1 UNIT 1, LEVEL 2 - KAUTH 

The KAUTH subroutine computes the Kauth Thomas data transformation for the 
patch image. 

CALLING SEQUENCE : CALL KAUTH(N) 

CALLING ARGUMENTS : 

Name Usage Type Definition 

N Process 1*2 Number of pixels to transform equals 1024. 

FILES : Not applicable. 

COMMON BLOCKS : The common block for the KAUTH subroutine is PATCH 

(see section 5 for description). 

SUBROUTINES CALLED : Not applicable. 

CALLED BY ; The KAUTH subroutine is called by routine HAZCOR (see 
section 5.4). 

LOCAL VARIABLES ; 

Name Definition 


Kauth 

matrix. 

as follows 


.332 

.603 

.676 

.263 

-.283 

-.660 

.577 

.388 

-.900 

.428 

.076 

-.040 

-.016 

.131 

-.452 

.882 
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5.5 UNIT 1. LEVEL 1 - VIN 


The subroutine VIN calculates a requested vegetative index number. The 
following list describes the vegetative indexes that can be requested and 
their request identification (ID) number. 


ID Vegetative 
number index 


Computation 


1 

2 

3 

4 

5 

6 

7 

8 


9 

10 

11 

12 


Greenness 

Green number 

3 AVI 

bpvi 

^TVI6 

^TVI7 


Leaf Area Index 
Simple Form 


-.283*CH1-.660*CH2+.577*CH3+.388*CH4 

(-.283*CHl-.660*CH2+.577*CH3+.388*CH4)+5 

2*CH4-CH2 

ABS(.38425*(2.4’^CH4-CH2)) . 

SQRT(ABS((CH3-CH2)/(CH3+CH2)+.5)) 
SQRT(ABS((2*CH4-CH2)/(2*CH4+CH2)+.5)) 
1.653*CHl/CH2-1.698*CHl/CH3+.093 


Leaf Area Index - 2.677-3.694*CHl/CH2-2.309*CHl/CH3 
Complex form +2.875*CH1/CH4+.043*CH2/CH3 

-1.346*CH2/CH4+3.017*((CH1/CH2) 

-(CH1/(2*CH4)))*CH1/CH2 


Berkeley 2*CH4/CH2 


®DVI CH2>2.4*CH4 

^RVI CH2/CH4 


Unused 


^Ashburn Vegetative Index 
^Perpendicular Vegetative Index 
‘'Transformed Vegetative Index 6 
‘^Transformed Vegetative Index 7 
®Di fferencfi Vegetative Index 

f 

Red Vegetative Index 
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CALLING SEQUENCE: CALL VIN(ID,N) 


CALLING ARGUMENTS : 

Name Usage Type Definition 

ID Process 1*2 One of the 12 identification numbers listed above. 

N Process 1*2 Number of possible pixel values to use in computa- 

tion equals 1024. 

FILES ; Not applicable. 

COfiMON BLOCKS : The common block for the VIN subroutine is PATCH (see 

5 for a description). 

SUBROUTINES CALLED : 

Name Definition ' evel 
$SQRT Square root System 

CALLED BY ; The VIN subroutine is called by PATCHP (see section 5 for a 
description). 

LOCAL VARIABLES : 

Name Definition 

ICNT Count of pixels used in calculating requested vegetative index. 


section 
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5.6 UNIT 1, LEVEL 1 - PATRPT 

The PATRPT subroutine formats the analyst aid stocistics file and writes a 
printer report of the contents. 

CALLING SEQUENCE ; CALL PATRPT(FILNAM) 


CALLING 

ARGUMENTS t 


Name 

Usage Type 

Definition 


FILNAM Array Byte File name of the analyst aid statistics file to 

be output js a report. 

FILES ; This PATRPT subroutine reads the analyst aid statistics file 
documented in section 5. 

COMMON BLOCK ; Not applicable. 

SUBROUTINES CALLED ; 

Name Meaning Level 
CL0S$ Close System 

OPENS Open System 

CALLED BY ; The PATRPT subroutine is called by the driver routine PATCHP 
(see section 5). 

LOCAL VARIABLES : Not applicable. 
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6. POSPP 


The post patch processor program POSTPP computes the spatial response function 
by regression onto a quadratic surface. 

CALLING SEQUENCE ; Not applica/e. 

CALLING ARGUMENTS : Not applicable. 


FILES : 

Fi le Usage Record format 

FILENM.DAT Analyst aid Statistics file; see PATCHP (section 5) for format. 


COMMON BLOCKS : Not applicable. 

SUBROUTINES CALLED : 

Name Meaning Level 

0PEN$ Open System 


CALLED 

BY : Not applicaole. 

LOCAL 

VARIABLES: 

Name 

Definition 

I SUM 

Inaccurate pixel count. 

RN** 

Normalizing constants. 

P(K) 

Orthogonal polynomials. 

C(K) 

Raw orthogonal coefficients. 

VINH 

Estimated variable value. 

R(I.J) 

Residual computation. 

SSR 

Sum of square residual. 
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Name 


Definition 


SUM Sum of residuals, 

RMSE Root mean square error. 

CN** Coefficients for orthogonal polynomials. 

SAS** Coefficients for simple polynomials. 

VAR Variable of interest. 
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APPENDIX A 


LISTING OF PATCH PROCESSOR 


riQiRlNAL PAGE IS 


FORTRAN 

IV-PLU3 

002-51E 14; 35; 06 30-OCT-80 PAGE 1 

PATCHP. 

FTN 

/TR: CLOCKS/WR 

0001 

C 

PROGRAM PATCHP 


c 

c 

THIS PR“6oRAM will DE the driver for the PAls. PROCESSOR 

0002 


COMMON/PATCH/CH, KH» FLAG. X. OAR. CLOUD, WATER. CAM, VEO, GIN. lAO, 
1 SAT. SLN, GAMA * 

0003 

0004 


INTEGER*2 lAG, SAT, CH( 1024, 4) , KH< 1024. 4) 
INTE0ER»2 GAR. CLOUD. WATER. CIN<3>. SLN 

0003 

0006 


DYTE FLAG <1024) 

REAL GAM, X(4), VEG< 12), GAMA 

0007 

0 

DYTE FILNAMOO), iHHNAM(30), IMDNAM<30) 


c 

c 

READ NAME OF HEADER FILE AND IMAGE DATA FILE 

oooa 

0009 

1000 

READ< 5, iOOO) IMHNAM 
FORMAT ( 30A 1 > 

0010 

001 1 


DO 10 1 = 1.30 

IF< IMHNAM< I ). NE. ' ')GOTO 10 

0012 

0013 


IMHNAM< I )=0 
GO TO 15 

0014 

0015 

10 

15 

CONTINUE 

READ(5, 1000)IMDNAM 

0016 

0017 


DO 20 1 = 1, 30 ■ ■ ’ 

IF< IMDNAMl I >. NE. ' ')G0 TO 20 

0018 

0019 


lMDNAM(t)=0 
GO TO 25 

0020 

20 

c 

CONTINUE 


c 

c 

CONSTRUCT NAME OF OUTPUT FILE USE NAME OF IMO FILE DUT 
SUBSTITUTE DAT FOR IMD 

0021 

c 

25 

DO 30 1=1. 30 

0022 

0023 


FILNAM( I )=IMDNAM< I > 

IF< IMDNAM( I >. EQ. ')G0 TO 40 

0024 

0025 

30 

40 

CONTINUE 
FILNAM( I+1)='D' 

0026 

0027 


FILNAM< I+2)='A' 
FILNAMl J+3)='T' 

0028 

0029 


FILNAMl I+4)=0 
IE=0 

0030 

0031 


CALL HEADER( IMHNAM, FILNAM, CH< 1, 1). lAG, SAT, IE) 
IF( IE, NE 'GO TO 990 


C 

C 

OPEN IMAGE DATA FILE 

0032 

C 

OPEN(UNIT=l, NAME=rMDNAM, TYPE='OLn', SHARED, READONLY. 


C 

1 ACCE3S= 'DIRECT ERR=900) 


c 

c 

LOOP FOR PROCESSING PATCHES 

o 

o 

CJ 

c 

DO 100 1PATCH=1,256 
WRIT£<5, 1235) I PATCH 

0034 

0035 

i235 

FORMATd.X, ’PROCESSING PATCH ',15) * 

CALL PATCRD( 1, IPATCH, CH< 1,1)) 

0036 

0037 


CALL GAMMA 
CALL HA2C0R 
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OPICalNAL HAGL iS 

OF POOR QUALITY 


i FORTRAN 

IV-PLUS 

V02-51E 14: 35; 06 30-0CT-80 PAGE 2 

PATCHP 

FTN 

/TR; DLQCKS/WR 

0038 

0039 


CALL VIN<2. 1024> - • 

CALL VINO. 1024) 

0040 

0041 


CALL VIN<4. 1024) " ' ' 

CALL VIN<6, 1024) 

0042 

0043 


CALL VIN<7, 1024) 

WRITE<2. 1 100) IPATCH, <X< I ), 1-1,4). OAR. CLOUD, WATER, GAMA, SLN, VE0<2), 

0044 

1100 

1 VEOO), VEG<4). VEG<6), VE0<7), (0IN< I ), I-1-.3) - 

FORMAT( IX. 13. 4F6, 2. 314, F6. 3. 13, 3F4. 1, 2F5. 3, 313) 

0045 

100 

C 

CONTINUE 


C 

C 

FINISHED 

0046 

0047 


CLOSE < UNIT-1) 
CL0SE(UNIT-2) 


C 

C 

WRITE PATCH REPORT 

0040 

C 

CALL PATRPT(FILNAM) 

0049 

C 

bo TO 990 


C 

C 

ERROR 

0050 

0051 

900 

1200 

WRITE <6, 1200) 

FORMATdX, 'ERROR IN OPENING IMAGE DATA FILE') 

0052 

0053 

990 

CONTINUE 

END 

. 
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ORIGINAL PAGE IS 
OF POOR QUALITY 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


1 FORTRAN 

IV- 

■PLUS V02-51E 14:35:22 30-OCT-80 PAGE 1 

HEADER 

FTN 

/TR- BL8CKS/WR 

0001 

C 

SUBROUTINE HEADER ( IMHNAM. F ILNAM, BUF, lAO. SAT. IE) 


c 

c 

THIS SUBROUTINE WIlL READ THE UNIVERSAL HEADER RECORD 
FROM THE IMH FILE, IT WIL... DECODE SOME INFO AND WRITE IT TO 


c 

c 

AN OUTPUT . DAT FILE 


c 

c 

CALL ARGUMENTS: 

IMHNAM I NAME OF IMH FILE 


c 

c 

FILNAM I NAME OF OUTPUT DATA FILE 

BUF I DUFFER TO USE TO READ HEADER DATA INTO 


c 

c 

I AG 0 SUN ANGLE OF ELEVATION 

SAT 0 SA TELL I TE NUMBER 


c 

c 

IE 0 ERROR INDICATOR SET TO 1 FOR ERROR 

0002 

0003 


BYTE IMHNAMOO), FILNAMOO). BUF<3072). BW0RD<2> 
INTEGER SAT 

0004 

0005 


EQUIVALENCE < IWORD, DWORD) 

DPEN<UNIT=1. NAME«IMHNAN. ACCESS* 'DIRECT READONLY. 


c 

1 SHARED. TYPE* 'OLD'. ERR = 100) 


c 

c 

READ HEADER AND DECODE DATA 

0006 

0007 


IW0RD*0 

READ! 1 '2> (BUF; I >. 1 = 1, 3060) 


c 

c 

READ SUN ANGLE OF ELEVATION BYTES 1833-1834 

0008 

c 

BWDRD< 1 )=DUF< 1934) 

0009 

0010 


BW0RD(2)=DUI- ( 1833) 
IAG=IWORD 


c 

c 

READ SATELLITE NUMBER BYTES <b5-66 

001 1 

c 

BWORD< 1 )=BUF(66) 

0012 

0013 


BW0RD(2)=DUF<65) 
SAT= IWORD 


c 

c 

READ LATITUDE AND LONGITUDE BYTES 1985-1995 

0014 

c 

BWORDM )=DUF< 1987) 

0015 

0016 


DW0RD<2)=DUF< 1986) 
LATD=IWORD 

0017 

0018 


DWORD ( 1 )=DUF( 1989) 
DW0RD<2)=DUF< 1983) 

0019 

0020 


LATM=IWORD 
DWORD! 1 )=DUF( 1992) 

0021 

0022 


BW0RD(2)=DUF( 1991 ) 
LOND^IWORD 

0023 

0024 


DWORD! 1 )=DUF! 1994) 
DWORD ! 2) =nUF! 1993) 

0025 

0026 


LONM=IWORD 
CLOSE !UNIT=1) 


c 

c 

WRITE DATA TO OUTPUT FILE 

0027 

c 

OPEN (UNIT=2, NAM£=F1L.NAM, TYPE* 'NEW ERR=200. 
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FORTRAN IV-PLU5 V02-5 1E 14 35:22 30-0CT-e0 PAGE 2 

HEADER. FTN /TR : QLOCKS/UR 

1 FORh= 'FORMATTED ' ) 

0029 URITE(2. 1000) <DUF( I ), I»33. 42)< lAO. SAT, DUF ( 1985). LATD< LATH. 

1 DUF( 1790) . LOND, LONM. <FILNAM( I ), 1-1.6) 

0029 1000 FORWAT( IX, lOAl, 12, I 1, Al, 14. 12. Al. 14. 12. 2X. 6A1. 44X) 

C 

C CHECK FOR GOOD DATA FROM HEADER 


0030 

C 

IF< lAO. LT. 15. OR. lAO, OT. 75) lE-l 

0031 

0032 


IF<SAT. LT. 1. OR. SAT CT 3) IE-1 
1F< IE. EO. 1 ) WRITE (6. 1300) I AG, SAT 

0033 

1300 

FORMAtUX, 'ERROR IN SUN ANGLE ='.I5,2X, 'OR IN SATELLITE NUMBER 
1 15) 

0034 

0035 


IF(SAT. EQ. 2)SAT=-1 
RETURN 


c 

c 

ERROR MESSAGES 

0036 

c 

100 

WRITE(6, 1100) 

0037 

0038 

1 100 

FORMAT < IX, 'ERROR IN OPENING HEADER FILE') 
IE-1 

0039 

0040 

200 

RETURN 

WRITE(6. 1200) 

0041 

0042 

1200 

FORMAT! IX. 'ERROR IN DPENINQ OUTPUT DATA Fii_E') 
IE-1 

0043 

0044 


RETURN 

END 





UKliaii VI* i>. I''m’uc. e1> 

OF POOR QUALITY 



A-6 
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FORTRAN IV-PLUS VQ2-51E 

HA-iCQR. FTM /TR: BLOCKS/WR' 


14; 3S; 38 


30-0CT-B0 


HAZCOR FTN 


SUBROUTINE HAZCOR 
IMPLICIT INTEGER <A-Z) 

COMMON/PATCH/CH» KH. FLAG, X, OAR» CLOUD. UATER. GAM. VEG, G1 
I SAT. SLN. GAMA 

INTEGER*2 lAG, SAT. CH< 1024, 4), KH< 1024, 4) 

INTEGER*2 CAR, CLOUD, WATER, CIN<3,‘. SLN 
DTtE FLAG<r02*4) 

REAL GAM, X(4), UEG< 12), GAMA 
R~EAL A ( 4 ) ,’ B ( 4 ) , C~OR , ANGR AD 

'type 100 

FORMAT! ^ HAZCOR DEBUG ON M 

GAMA* ( GAM-300 ) / 1 000 

I F<GAMA. GT. 0. 7 .OR. GA M A. LT. -0. 7 ) GO TO 333 

IF<IAG LT. 15 .OR. IAG'GT. 75) GO TO 333 

ANGRAD-1AG*3, 1416/180 
COR* 77713/SIN<ANGRAD) 

A( 1 )=£XP<-1. 26a 0»QAMA> 

A<2)=EXP<-1. 0910*GAMA) 

A<3 )= EXP 1-0. S 358»QAMA) 

A<4’)=EXP(-0. 59Br»GAMA) 

B( 1 )*61. 9«( 1-A< 1 ) )ttCOR 

B(2)=66. 2*( 1-A<2) >*COR 

D(3)=e3. 2*( 1-A<3> )*COR ^ 

B<4)»33, 9*( 1-A<4) )*CQR 
A< 1 )=A( 1 )*COR 
A(2)*A<2)*C0R 
A<3)=A<3>*CQR 
■ A ( 4 ) =A ( 4 ) *CQR 
TYPE 102, A, B 

FORMAT!' HAZCOR A AND B',/ 8F6. 3) 


DO 2 K=l, 1024 

CH!K, 1 )=CH!K, 1 )#A! 1 >+B! 1 ) 
CH!K, 2)=CH!K, 2>*A!2)+B<2) 
CH!K-, 3)=CH!K, 3 ) ^»A ! 3 ) +B ! 3) 
CmK, 4)=CH!K, 4)^*A!4)-»-B!4) 
CONTINUE 
CALL K.AUTH!1024) 

SAT*“2 

IAG=51 

QAM=300 

RETURN 

TYPE" 200. GAM, GAMA, I A0‘, SAT 
FORMAT!' HAZCOR PROBLEMS', 
' GAM*', 16,/, 

' GAMA* ' , F8. 4, /, 


TO CONTINUE 


ACCEPT 201, I IT 


FQ,RMAT!A1 ) 
WRITE!6, 2100) 


S TO STOP') 
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FORTRAN IV-PLUS V02-51E 
HAZCOR FTN ' /TR. "blocks/ WR 


14 35 36 


3O-OCT-0O 


FORMATUX, 'HAZCOR PROBLEMS ENCOUNTERED') 

IF( 1 I T . Ea 'S' ) STOP 

CO TO 'l 
END 
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■ TT17: C3g0. 403A D/-SP=SYO: t320. 403PATCRD 



mi?- 


OF POCH QUALi J 




FORTRAN 

1 0-PLUS 

U02-51E 14.35.54 20-OCT-30 PAGE i 

PATCRD. FTN 

/TR: DLOCKS/WR 

0001 

c 

SUQRQUtiNE (*AT6rtD(LUN. IPA'ttR, ftu*^k 


c 

c 

TH 1 i SUBROUTINE WILL READ ONE tHANNEU 61^ tSNE PATCH IMAGE 


c 

c 

CALL ■Art'<ruHENTS: 

LUN I LOGICAL UNIT NUMBER TO READ PATCH INFO 


c 

c 

IPATCH I PATCH NUMBER REGUE'^ED 

BUF 0 BUFFER FOR DATA FROM PATCH 

0002 

c 

BYTE DDUF<21 ia>, QW0RD<2> 

0003 

0004 


INTEGER IDUF( 1059), BUF( 1024. 4) 
EQUIVALENCE l IQUF, DBUF) , < IWORD, BWORD) 

0005 

c 

IW0RD=»0 


c 

c 

CALCULATE RECCRD NUMBER FOR FIRST LINE OF PATCH DATA 
ILIN=LINE NUMBER OF IMAGE WHERE PATCH DATA BEGINS 


c 

c 

IPIX=PIXEL NUMBER FOR START OF PATCH DATA IN ILIN 

0006 

0007 


ILIN=< IPATCH-1 )/ 16*32+1 
IREC»ILIN 

oooa 

0009 


INDXl»l 

IPIX*MOD< ( IPATCH-1 ), 16)*32+70 


c 

c 

DO LOOP TO READ 32 LINES OF PATCH 

0010 

c 

DO 100 1=1,32 

0011 

c 

READ < LUN 'IREC .) IBUF 


c 

c 

DO LOOP TO MOVE 32 PIXELS TO OUTPUT BUFFER CONVERTING TO 1*2 

0012 

0013 


DO 50 J=l, 32 

BWORD( 1 )=BDUF< IPIX+J) 

0014 

0015 


UUF< INDXl, 1 )=IWURD 
DWORD! 1 )=DBUF( IPIX+U+512) 

0016 

0017 

- 

DUF( INDXl, 2)=IW0RD 
DWORD! 1 )=BBUF! IP I X+J+10r'.4 > 

0013 

0019 


BUF! INDXl, 3)=IW0RD 
BWORD! 1 )=BBUF( IPT:<+J+1536> 

0020 

0021 


BUF! INDXl, 4 >-rWORD 
INDX1--INDX1 + 1 

0022 

0023 

50 

CONTINUE 

IREC=IREC+1 

0024 

0025 

100 

CONTINUE 

RETURN 

0026 


END 
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ORIGINAL Fiv 
OF POOR QUA 


FOPTRAN 

IV- 

PLUS V02-51E 14:36 22 30-0CT-80 PAGE 1 

VIN. FTN 


/TH: BLQCKS/WR 


C VIN. FTN 

0001 


SUBROUTINE VIN< ID. N) 

0002 

C 

IMPLICIT INTEGER <A-Z) 


c 

1 GREENNESS 


c 

2 GREEN NUMBER 


c 

3 ASH3URN VEOITATIVE INDEX 


c 

4 PERPENDICULAR VEGITATIDN INDEX 


c 

5 TRANSFORMED VEGITATION INDEX 6 TV I 6 


c 

6 TRANSFORMED VEGITATION INDEX 7 TVI7 


c 

7 LEAF AREA INDEX SIMPLE, FORM 


c 

a LEAF AREA INDEX COMPLEX FORM 


c 

9 BEREKLY 2»CH4/CH2 


c 

10 DVI 


c 

11 RVI 


c 

12 UNUSED 

0003 

c 

COMMON/PATCH, ''CH, KH. FLAG, X, GAR, CLOUD, WATER, GAM. VEG. GIN, lAC, 



1 SAT, SLN, GAMA 

0004 


INTECER»2 lAO, SAT, CH(1024, 4) , KH( 1024, 4) 

0005 


INTEGER.'t2 GAR, CLOUD, WATER, GIN(3) , SLN 

0006 


BYTE FLAG(1024) 

0007 


REA), GAM, X(4) .VEG'<12), GAMA 

0008 


REAL. VEON 

0009 

0 

VEG (I D > =0 

0010 


icn't=o 

OOl 1 


IF(ID. LT. 1 .OR. ID. GT. 12) GO TO 333 

0012 

c 

GO TO <1,2,3,4,5,6,7,8,9,10,11,12) ID 

0013 

1 

CONTINUE 

0014 


DO 15 1=1, N 

0015 


IF(FLAG<I), NE. UGO TO 15 

0016 


VEGN=-, 2B3*CH( I, 1 )-. 660-»CH< I, 2)+. 577*CH( I. 3)+. 3Sa»CH< I, 4) 

0017 


VEG( 1 )=VEO< 1 )+VEGN 

OOIB 


TCNT=ICNT+1 

0019 

15 

CONTINUE 

0020 


GO TO 200 

0021 

2 

CONTINUE 

0022 


DO 20 1=1, N 

0023 


IF(FLAO< I ) . NE. i;G0 TO 20 

0024 


VEGN=-. 283i^CH( I, 1 ) -, 660#CH< I, 2)+. 577»CH( I, 3)+. 3S3^CH( I, 4) 

0020 


VEGN=VEGN+5 

0026 


IF(V£GN. LT 0)VEPN=0 

0027 


VE0<2)=VEG(2)+-V£GN 

0020 


ICNT=ICNT+1 

0029 

20 

CONTINUE 

0030 


GO TO 200 

0031 

3 

CONTINUE 

0032 


DO 30 1=1, N 

0033 


i F ( FLAG < I ) . NE. 1) GO TO 30 

0034 


VEGN=2»CH( I, 4)-CH< I, 2) 

0035 


IF<VEGN . LT. 0. ) VEGN=0. 

0036 


VEG(3)=VEG<3)+VEGN 

0037 


ICNT=ICNT+1 

0030 

30 

CONTINUE 
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oRiGir'iAi. r;. 

OF POOR QUA ' . 


i FORTRAN IV-PLUS V02-51E 14 36,22 30-aCT-B0 PAGE 2 

VIN. FTN /TR' dLoCKS/WR ~ ' 


0039 



GO TO 200 

0040 

4 


CONTINUE 

0041 



DO 40 I-l.N 

0042 



IF<FLAG< I). NE. 1 )G0 TO 40 

0043 



VEGN»ABS< 38425* <2. 4*CH( I, 4)-CH< I, 2) ) ) 

0044 



IF<VEGN, LT, 0)VEGN«0 

0045 



UEG<4)=VFG<4)+VEON 

0046 



ICNT=ICNT*1 

0047 

40 


CONTINUE 

0049 



GO TO 200 

0049 

5 


CONTINUE 

0050 



DO 50 I-1,N 

0051 



IF(FLAG< I) . NE. 1 )G0 TO 50 

0052 



VEGN=SQRT ( ABS ( ( CH < 1 , 3 ) -CH < 1 , 2 ) ) / ( CH ( 1 , 3 ) +CH ( 1 . 2 ) ) +0 5 > ) 

0053 



VEG( 5)=VEG<5)+VEGN 

0054 



ICNT=ICNT+1 

0055 

50 


CONTINUE _ . 

0056 



GO TO 200 

0057 

6 


CONTINUE 

0050 



DO 60 I»l, N 

0059 



IF ( FLAG < I ), NE. 1 )G0 TO 60 

0060 



VEGN=SQRT(AQS( (2*CH< I, 4)-CH< 1.2) ) /(2*CH< I, 4)+CH< I, 2) )*0, 5) > 

0061 



IF( (2*CH< 1, 4)-CH( I, 2) ). LT. 0)VEGN=0 

0062 



VEG(6)=VEG<6)+VEGN 

0063 



ICNT=ICNT+1 

0064 

60 


CONTINUE 

0065 



GO TO 200 

0066 

7 


CONTINUE 

0067 



DO 70 1 = 1, N 

0069 



IF<FLAG< I ). NE 1)00 TO 70 

0069 



VEGN=1. 653^CH<I, 1)/CH( I, 2)-l. 693«-CH<I. 1)/CH< I, 3)-t-, 093 

0070 



VEG(7)=VEG(7)+UEGN 

0071 



ICNT=ICNT+1 

0072 

70 


CONTINUE 

0073 



GO TO 200 

0074 

8 


CONTINUE 

0075 



DO 80 I--=l, N 

0076 



TF<FLAG{ I ), NE 1 )G0 TO SO 

0077 



VEGN=2. 677 - 3, 694 :tCH ( 1 , 1 ) /CH < 1 , 2 ) 




- 2. 309tt-CH< I, 1 )/CH( I, 3) 




+ 2. 375.^CH< 1,1) /CH( I, 4) 




+ 0. 043«CH< I, 2)/CH( I, 3) 



♦ 

- 1. 346^i•CH( I, 2) /CH( 1,4) 




+3. 017*( (CHt 1, 1 )/CH(I, 2) )-<CH( I, 1 )/(2*CH(I, 4) ) ) > 



* 

«CH( I, 1 )/CH< 1,2) 

0079 



VEG(8)=VEG(0)+UEGN 

0079 



ICNT=ICNT+1 

0030 

00 


CONTINUE 

0031 



GO TO 200 

0082 

9 


CONTINUE 

0083 



DO 90 1 = 1, N 

0084 



IF ( FLAG < I ). NE. 1 )G0 TO 70 

0035 



VEGN=2»CH( I, 4)/CH( I, 2) ' 

0086 



VE0(9)=UEG(9)+VEGN 

0037 



ICNT=ICNT+1 

0088 

90 


CONTINUE 
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POOR 


FORTRAN 

IV-PLUS 

V02-51E 14:36:22 30-0CT~a0 PAGE 3 

VIN. F TN 


/TR: BLOCKS/WR 

S089 


go" to 200 

0090 

10 

CONTINUE 

009 i 


00 100 I-l.N 

0092 


IFCFLAGU ). NE. 1 )0C TO 100 

5093 


VECN«CH<r, 2>-2, 4*CH< I. 4) 

0094 


VEG< 10>»VEG( 10)+VEGN 

5093 


ICNT-lCNT+1 

009A 

100 

CONTINUE 

0097 


GO TO 200 

0090 

11 

CONTINUE 



DO 110 I=»1,N 

OlOO 


IF ( FLAG < 1 ). NE. 1 )GO TO 110 

0101 


VEGN=CH(I. 2)/CH( I, 4) 

0102 


VEG( 1 1 )*VEG< 1 1 )-<-VEGN 

0103 


ICNT*ICNT+1 

0104 

no 

CONTINUE 

0105 


GO TO 200 

0106 

12 

CONTINUE 

0107 


GO TO 200 • " ■" 

wmmmm 

333 

WRITE(6. 334) ID 

0109 

334 

FORMAT < ' DAD ID ='>18,' SUBROUTINE VIN') 

Olio 


VEG( ID) =9999. 

0111 


RETURN 

01 12 

o 

o 

IFHCNT. NE. 0)VEG< ID)=VEG( ID)/ICNT 

0113 


IFdCNT. EQ. 0)VEG< ID)=0 

01 14 


RETURN 

G1 1 5 


END 
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ctions and subroutines referenced 



FORTRAN IV-PLUS V02-51E 14:36:22 30-0CT-Q0 RASE 5 

VIN. FTN /TR BLOCKS/WR 

TOTAL SPACE ALLOCATED - 045334 9582 

. TT17 C320< 4Q3A. B/-SP»SYQ: C320. 403V1N 



OF POOH 


* . 




+VIUI } 


\ 

1 FORTRAN 

IV- 

PLUS V02-51E 14 36:42 30-OCT-80 PACE 1 

GAMMA. FTN 

/TR, BLOCKS/WR 


c 

c 

GAMMA, FTN 

0001 

0002 


SUBROUTINE GAMMA 

COMMON/PATCH/CH. KH, flag. X, CAR, CLOUD. UATER, CAM, VEC. GIN, lAC. 

0003 


1 SAT, SUN, GAMA 

INTEGER*2 lAG. SAT. CH( 1024, 4), KH( 1024, 4) 

0004 

0005 


I NTEGER »2 GAR , CLOUD, WATER , G I N < 3 ) , SLN 
BYTE FLAG (1024) 

0006 

c 

REAL GAM, X<4), YEO( 12), GAMA 

gamma computes the xstar optical depth parameter gamma 


c 

c 

lAG IS SUN ELEVATION ANGLE IN DEGREES 
GMA IS THE HAZE PARAMETER OUTPUT 

0007 

c. . 

REAL AA, BD, CC, DD, XX(4), Y<5). Z;(4), ANO, S 


c. . 



D 

),RITE(6. 16) lAG, SAT 

oooa 

0009 

16 

FORMAT!' GAMA IAG»',I5, ' SAT=',I5) 

IFdAG . LT. 15 OR lAG . GT, 75) GO TO 333 

0010 
001 1 

1 

IF (SAT . LT. -2 .OR. SAT . GT 3) GO TO 333 
ANO-IAG »3. 14159 /ISO. 0 

0012 

0013 


CR*. 77713/SIN(AN0) 
LK=3AT 

0014 

0015 


CALL CALI3(LK. 1024) 
CALL SCRALL( 1024) 

0016 

0017 


CALL MEAN! 1, 1024, NN, S) 
CALL MEAN(2, 1024, NN, S) 

OOlS 

0019 


CALL MEAN(3, 1024, NN.S) 
CALL MEAN(4, 1024, NN, S) 

0020 

0021 


GO TO 567 

ENTRY KGMA(Y) 

0022 

0023 


X( 1 )=Y( 1 ) 
X(2)=Y(2> 

0024 

0025 


X(3)=Y(3) 

X(4)=Y(4) 

0026 

0027 

CR -Y(5) 
567 CONTINUE 

0028 

D 

17 

WRITE(6, 17) X 

FORMAT! ' GAMMA MEANS 4F8. 3) 

0029 

0030 


Z!l)=CRi^X! 1) 
Z!2)=CR*X!2) 

0031 

0032 


Z!3)=CR*X!3) 

Z!4)=CR*X!4) 

0023 

oorj'f 


AA=1. 6078^KZ! 1 )-61, 9)#(-0. b995) 
AA-AA>I. 0910*! Z !2)-66. 2)*0. 4283 

0035 

0036 


AA=Am+0. 8358*! 2 !3) -S3. 2)*0, 0759 
AA=AA“0. 5901*(2!4)-33 ?)*0. 0408 

0037 

0038 


BB--1, 2680*! Z ! 1 )-61. 9)*0. 8995 
EB = I3I3 + 1. 0445*!Z(2)-66 2)*0. 4283 

0039 

0040 


BB=I3D+0 9142*!Z!3)-83. 2)*0, 0759 
DB--13D-0. 773'l*!Z!4)-33. 9)*0 0408 

0041 

0042 


CC--2! 1 )»0. S9'i’5+2!2)*0. 4283*Z ! 3 ) *0. 0757-Z!4)*0. 0408*11, 2032 
DD-BD#DD-2. 0*AA*CC 

0043 

0044 


IF!DD . LT. 0. 0) DD=0. 0 
GAMA=!BD-S0RT'DD) )/AA 

0045 

D 

WRITE! 6, 2020) GAMA 
GAM=GAMA* 1 000+300 
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ORiG'Ti;*'. r.’.’ •. . 
OF FQO;-^ / 


FORTRAN 1^^-PLU5 V0 2 -51E 14: 3A: 42 30-QCT-90 PAGE 2 

GAMMA. FTN /TR DLOCKS/WR 


0046 

0047 

2020 

FORMA r(//» ' 
RETURN 

gamma - F8. 4. / / ) 


0048 

0049 

333 

8 

TYPE a. lAO. SAT. NN 

FORMAT(' PROBLEMS rAC-',I6, ' SAT*', 16, ' 

NN* ', 16. 

OOSO 


* /, ' 
ACCEPT 9, I ITT 

S TO STOP CR TO COrJTINUE', lOX. 

'GAMMA ' ) 

0051 

0052 

9 

FORMAT (A1 ) 

IF< IITT. EQ 'S' 

> STOP 


0053 

0054 


GO TO i 
END 







OF POOR QUALITY 










1 


1 


fW 

O ^ 






















(A 



if. 

r« 

rj < 






















in 



(Al O 

o g 






















UJ 



u 

ru 

(M C 






















X 



X 

«r 

^ 0 






















S 



2 

o 

O Q 






















g 



c 

1 

1 T 











(A 











< 



< 













CA 



























Ui 

9 


























X 

9 










Ui 



III 













a 











0. 



(L 

TiJ 

CM ♦ 











q 











> 



> 

« 

9 9 











< 















M. 

M« C 


























a 












-I 











w 
















UJ 














c 

Q 

Qf 











m 

* 










2 



3 

u 

V 











< 

h. 













Z 

u 

a (A 











mJ 















o 

•A «0 






















(A 



UI 
























(A 



UI 

o 

Q 






















Ui 



lij 

o 

r< g 






















s 



X 

o 

*t 0 






















Q 



Q 

p 

O 0 






















o 



Q 

T 

1 1 











(A 











< 



< 

•t 

•A ♦ 











(A 



























Ui 

9 


























¥ 

9 










Ui 



UJ 













a 











a. 



a 

«r 

<r (M 











Q 











> 



> 

* 

9 9 











< 











K 



h 

(X 

X 






































J 











UJ 





< 











Ui 











c 



f. 


c 











n 

% 










5 



<. 

o 

< ^ 











< 

A 










Z 



z 

o 

o z 











-1 


n 














•A ♦ 

♦ 





















(A 



cn 

o 

rj N 

rj 












Ui 









(A 



U) 

o 

o o 

o 












C9 









Ui 



UJ 

p 

(M O 

CM 












< 









X 



X 

o 

^ O 













CL 









Q 



o 

o 

O O 

p 











O 










Q 



Q 

1 

1 1 

T 










CA 

n 










< 



< 


■A f 

A 










(A 

«4 


























Ui 

9 * 


























K 

n 9 










UJ 



liJ 













a 

o 










0. 



0. 


W CM 

N 










q 

1 










>• 



>• 

« 

9 9 

9 










5 

n 










H- 



h 

a: 

X ^ 




























Z 










,i 


O 









(i/ 



Ui 



UJ 










Ui 

o 

CD 





-i -j 

-I -J 



r 



i: 


r 

1- 










m 

Tij 

j 




CJ u 

u u 

C C3 



< 



< 

a 

< X 

< 










< 

' n 




Cfl 


J mJ 

J O 



z 



z 

CB 

O .J 

3 


(f) 









0 ^ ry 

U 



Ui 
























o 



K 

z Z 

Z z 

Z X 











u 










1 



D 

o o 

o o 

o > 














«r 







o 



Q 

(J u 

u u 

o o 




*0 



o 

«r o 



cn 

, 









n 



M 

. . 

, . 

» « 



cn 

n 


tft 

>0 

- o 

fJ 


z 













or 

« Q 

Q Q 

a Q 



cn 

Cii 


U) 

o 

o 



UJ 

rtf 

rj 

CM ^ 










K 

« » 


« ^ 



Ui 

o 


Ui 

o 

o o 

Tii 


E 

OlO ^ 

o rj 

ox. >M, 

MM |Mt 








h 

3 3 

3 3 

3 3 



L« 

o 


X 

c 

o o 




rt 

w \ / 



in ^ 








< 

tr CE 

(X cc 

q: (£ 



Q 

o 


Ct 

o 

o o 

O 


c 



W *M> 


W 




o *t 

fV 









a 

1 


Q 

1 

I ( 

1 










(A 

o n 

•r 










MM 


< 

•y 

♦ 

A 










cn 

O n 


























UJ 

o o 


























X 

o o 

n 









UJ 



lit 






Aioj n 

A <r 

CO 03 

o m 



a 

o o 










0. 



a. 


«r CM 

CM 



o 

r- 

Ch CM 


mM 



Q 

1 1 










> 



> 

« 

r 9 

9 



oi»n 

o 





< 

tn 













i" 

X 

(C 






























Ui 























































OlO A 

o o 

o o 

«T O 






S 








UJ 

< 


Ui 


K- 


cn 

o 

o o 

o A.fj ^*dirj rj 



UJ 



3 



rs o 

nO n 

CM O 



i: 

c 


i: 

o 

h-lZ 


o o c 

o o 

o o 

O O 



ca 



s. 




^ n 




< 

o 


•t 


o 

-1 


o rj o 

O O'C o 

o o 



< 

' A 





n 





z 



z 

< 

a 

in 1 


“J o o 

o o 0‘0 o 



mJ 

m in 







03 











OlO o 

O 0,0 o. 

o o 






L) 


Ui 












1 











UJIQ 


fsi 












1 






9 





•^1 

-/ 


MM 

-0 n- 

O (M O (M 



o 



o 

O OfCM 1 


0.0 rJ 

, 

o rji 

o o 

rj o 





Xi 

a 


cn 

•o ^ 

n o 

n n 


l(n 

o 


CD 

o 

rN ralra 

cn 

o 

o 

o n o rji 

O Mf 





1 .. 



ru -irj — 

o -• 


'tn 

o 


CO 

o 

O 

mM 


bi 

o 

o 

o o 

o o 

o o 





tv a: 



-• o 

o olo rj 


!uj 

o 


u 

o 

o rjiCJ 


»*j 

o 

O (u 

o rj rtf oio o 





o 




c o 

o o 

o «r 


,(£ 

o 


X 

o 

O T 

«T 


a* 


y *r 

CJ *r 

*r C'- 

o o 

. 







o o 

O Olo o 


l<5 

o 


a 

o 

o o 

o 


Q( 

0,0 o 

o o 

n o 

o o 




rtf «t 



tn 






,.4 

1 


a 

1 

1 1 

1 


a 


i 1 

\ 1 

' ' ; 

t 1 



01 

N 

U) 


z 








OT« 


< 

▼ 

^ A 

A 


< 

A] 

A A 

A Al 

A ^ 

X <J 



'JJ 

O -* 



□ 


MM < 

< 

cn 














1 

1 



ILU 1 




M 

111 

UJ h 

K (n 

iX I 












; 



j 





o o 

a>, 



c 

a *t 

< OTIC u 


01 

III 



UJ 





UJ 




1 

, 



a 

o o 

1 


u 

*? 

cj a 

n < 

Ui K 


h— 

0. 



a 

•r 

<T rj 

ni 


a 

rj’ 

-M CJ 

rj 

•r it! 

*7 *r 



a 

1 1 



XI 

z 

U IL 

2> 



z 

> 



> 

« 

9 * 

■ 


> 

■9 1 

« 9 

• 9| 

9 9J • *1 



< 

mM tM 


z 



<A «» 

# «» 

a 


MM 

1- 


fn 


q: 

rr *-M* 





-J 

~ ir 

ft X 

Z X 






•• 







o 



Ui 

, 















r: 

X 

z: 






0. 












i 







T 


< 

(t 







•z 

r-4 

■ 




cn 


I 


J 




tn 

-J 


c 

<c 

x 

UJ 





> 

UJ 

z 

< 

u 



1 

> 

iji 

j 






,j 

u 


K'x: ' 


ca 





ir 

j: 

r 

MM 



C* 

K 

•j: 

i: 

j 

•t Z 

<•5 




UJ 


n 

Si 


O 

Z 

-• rtf 

3 

4 




< 

< 

ir 

•( 

< 

X < 

< 

z 

/J. 

2, 


X [Ij 

X 



u 

< 

n 



X 

D 





2 


o 

<1 

2 

< 

o ^ 

0) 

X 

T* 

O 

X O 

i. > 

X xi 

> fM 


< 


MM n 

u. 

“ 1 

0. 

Z 





Ui 



->1 

1 





< 








-J 




A-22 


FUNCTIONS AND SUDNOUTINES »IFFERENCEO 


ORIGINAL PAGf: IJ 
OF POOR QUALITY 


FORTRAN IV-PL'J5 W2-5 1E 14:36 42 30-0CT->90 PACE 4 

OAMMA FTN /TR. DLOCKS/WR 


TOTAL space ALLOCATED - 044146 9567 

. TT17: C320, 403A. B/-SP«5Y0; C320, 4030AMMA 
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FUHCTIONS AND SUDROUTINES REFERENCED 



ORIGINAL PAGE IS 
OF POOR QUALITY 


^CRTRAN IV;;2PLUS V02-51E, 

J^AuTH FTn ”/TR aUDCKS/WR 


14 37 46 


DO-OCT-BO 


KAUTH COMPUTES THE KAUTH THOMAS TRANSFORMATION 

SUD RQ UTINE KAUTH(N) 

IMpXIcIT INfECERtS (A-S) 

COMMON/PATCH/CH. KH. FLAG- X- GAR, CLOUD. WATER. GAM, VEC, GIN, lAC. 


1 SAT, SLN, GAMA 

INTEGER«2 lAC, SAT, CH( 1024, 4), KH( 1024, 4 ) 


INTEGER*2 CAR, CLOUD, WATER. CINO), SLN 
BYTE FLAG < 1024) 


REAL CAN, X ( 4 ) , VEG < 1 2 ) , GAMA 
REAL KHC(4) 


REAL R(4, 4) 

DATA R / 332, 603, 676, 263, * 


- 233, - 660, 577, 388, 

- 900^423, 076, - 040, 

C -"016*, 13ir-“45'2, 882/ 

WRITE (6, 66) R 






ORIGINAL rV/.L. . 
OF POOR QUALITY 



1 
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TOTAL SPACE ALLOCATED » 0426ia B901 


ORIGINAL PAGE IS 
OF POOR QUALITY 


i FORTRAN 

lY- 

-PLUS U02-51E 14.40,07 30-0CT-B0 PAGE 1 

3CRALL. FTN 

/TR; BLOCKS/WR 


C 

C 

SUBROUTINE SCR ALL KAUTI-' TRANSFORMS AND SCREENS TAfA 


c 

c 

COMPILE WITH THE CO. 90 SWITCH. 

0001 


SUBROUTINE SCRALL(N) 

0002 


IMPLICIT INTEGER (A-Z) 

0003 


COMMON/PATCH/CH, KH. FLAG. X, GAR, CLOUD. WATER, GAM. VEO. GIN. IA6. 
1 SAT. SLN. GAMA 


0004 1NTE-CER*2 lAG. SAT. CH< 1024. 4 ) > KH( 1024. 4 ) 

0009 INTEQER*2 QAR. CLOUD. WATER. QINO). SLN 


BYTE FLAG (1024) 

REAL GAM. X<4). VEQ( 12). GAM A 

REAL R<4. 4). ZZ.C. Z(4. 19). RR<4. 4). Y<4),CC(19> 

INTEGER 11(19). JJ(17).LL(21). XX(4). T 

INTEGER*4 14 

INTEGER»2 55N 

DATA RR/. 33231. . 60316. . 675B1. . 26278. 

C "■ 2831 7. 66006. . 57739. . 38933. 

C 89992. . 42830. . 07592. 04080. 

C 01594, , 1306B. 45187. . 88232/ 

0013 DATA Z/0. 0. 0. 0. 0. 0. 1. 0. 

C 0. 0. 0. 0.0, 0. -1. 0. 

C -.09379.0.-0.1.0.0.0. 

C -■ 1879. 0. 0. -1. 0. 0. 0. 

C -0. 1. -1. 0. 0. 0. 0. 0. 

C 0. 999996. 1. 0. 0. 0. 0. 0. 

C -0. 83333333. 1.0.0. 0. 0. 0. 

C 1. 0. 0. 0. 0. 0. 0. 0. 

C 1.0. 0. -1. 0. 0. 0. 

C 1. 0. 0. 0, 0. 0. 0. 0. 

C -0.142857,0.0,-1.0.0.0. 

C -1.0.0.0.00.0.0. 

C -.0679,-1.0,0.0.0.0. 

C 0 0. 0. 0, 0. 0. -1. 0. 

C 0. 0. -1. 0, 0. 0. -1. 0. 

C -. 9. -1. 0. -1. 0. -0. 5. 


C -. 4. 1. 0. 0. 0. 0. 0. 

C -0. 4. 1. 0. -0. 6. -0. 6. 



0014 


C -1. 0. 0. 4, 0. 0. 0. 0/ 

DATA CC/16. , 12. . -4. . 14. . 20. . 156. , -8. , 100. , 7. 5. 69. . 3. 29. 




C -79. . 0. 5. -1. 9. 4. 5. -10. , -12. 2, -9. . -37. 79/ 

• 

0019 


DATA II/2, 3. 4. 9. 6. 7. 8. 10. 10. 12. 12. 18. 18. 18, 18. 18. 60. 20. 20/ 


0016 


DATA JJ/60.60, 60. 60.60, 60. 60. 9.60. 11.60. 13. 14. 19. 16. 17.21. 




C 19,60/ 


0017 


DATA LL/47. 46. 49. 44, 43. 42. 41. 34. 33. 32. 31. 26. 25. 24. 23. 22. 21. 




» 13. 12. 1. 11/ 



c 

4« GARBLED; 34.33 CLOUD; 32.31 HAZE 



C 

2« WATER; 1 * SHADOW, 11 OVER WATER; 1 GOOD. 

13 


C 


II _ 


D 

TYPE 808 

10 

0018 

808 

FORMAT (' SCRALL') 

• 

0019 


SSN«30 

• 

0020 


GAR-0 

7 

0021 


CLOUD-0 

C 

0022 


WATER -0 

t_ 

0023 


0IN(l)-0 


0006 

0007 

0008 

0009 

0010 
0011 
0012 
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ORIGINAL Fi ' - 




OF 


tAl.lTY 


1 FORTRAN 

IV-PLUS V02-51E 14:40:07 30-0CT-BO PACE 2 

SCR ALL. FTN 

/TR. BLQCKS/WR 

0024 


G1N<2)-6 

0025 


0IN(3>«0 

0026 


lCNT-0 

0027 

600 

IF<IAC. LT. 15. OR. I AG. OT. 80 » TYPE 101. 8EG. DATE 

0028 

101 

FORMAT < ' INPUT SUNANCLF FOR M6. 2t3. ' SCRALL') 

0029 


IF( XAG. LT. 15. OR. I AO. OT. ACCEPT 201 . lAO 

0030 

201 

FORMAT! 12) 

0031 


Z2*»IAO/57. 3 

0032 


IF<2Z. LT. 0. 1) ZZ-6. «^07 

0033 


C»0. 77713/SIN(7.Z> 


0 

TYPE 990. lAO. C 


D 

WRITERS. 990) lAO.C 

0034 


IF (C. LT. 0. J ' TYPE 990. lAO. C 

0035 


IF <C. LT, 0. i ) RETURN 

0036 

990 

FORMAT!' C M4. FIO. 4) 

0037 


DO 100 K-l.N 

0038 


03 5 1-1.4 

0039 

5 

Y!i)-0 

0040 


XX<1)-CH(K. 1> 

0041 


XX!2)>CH!K. 2) 

0042 


XX!3)-CH!K. 3) 

0043 


XX!4)-CH!K. 4) 


0 

IF!K. EG. 1) TYPE 995, X 

0044 

995 

FORMAT!' X '.415) 

0045 


00 6 1-1.4 

0046 


DO 6 J-1.4 

0047 

6 

Y ! I ) -Y ! 1 ) +XX ( J ) *RR ! J, I ) * 

0048 


ISN-Y!2)-»-. 5 

0049 


IF! ISN. LT. SSN)SSN»ISN 

0050 

98 

DO 66 1-1.4 


66 

Y(I)»Y<I)*C ! SUN ANGLE CORR^CTiON 

0052 


L-1 

0053 

1 

CONTINUE 


D 

IF!K. LT. 3) TYPE 996, L. T. U. CC!L) 


D 

IF(K. LT. 20) WRITEO. 996) L. T. U. CC !L). Y 

0054 

994 

FORMAT!' 14. 14, 5F10. 4. FIO. 4) 

0055 


IF(L. LE. 21)T»LL!L) 

0056 


IF!L . OE. 20) GO TO 99 

0057 


U-0. 

0058 


DO 7 1-1.4 

0059 

7 

U-U^-Y! I )*Z(I.L) 

0060 


IF!U . LT. CC!L) ) GO TO 3 

0061 


L-JJ!L) 

0062 


GO TO 1 

0063 

3 

CONTINUE 

0064 


L-II!L) 

0065 


GO TO 1 

0066 

99 

FLAO<K)-T 

0067 


IF!T. EQ. DICNT-ICNT+l 

0068 


IF!T. OE. 40. AND. T. LE. 49)GAR»GAR-H 

0069 


IFiT. GE. 30. AND T. LE. 39)CL0UD*CL0UD+1 

0070 


IF(T. GE. 10. AND. T. LE. 29) WATER -WATER +1 

0071 


IF!T. NE. UGO TO 100 

0072 


IF<Y(2). GE. 25. )0IN(3)«GIN(3)-H 

0073 


IF(Y!2). GE. 10. )GIN(2)-GIN(2)>1 

0074 
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ORtQH'iAL PAGE i- 

OF POOR QUALITY 


FORTRAN IV-PLUS V02-9IE 14; 40:07 

30-0CT-80 

PACE 3 

SCRAL FTN 

/TR; BLOCHS/ WR 



0075 • 100 

0076 

CONTINUE 

SLN-SSN 



0077 

0079 

IF<SLN. LT. -10>8LN— 10 
XFdCNT. LE. 0) RETURN 

. 


0079 

0090 

I4-0XN<1) 

14»X4*100 



0091 

0092 

0IN(1)«I4/1CNT 

14*GIN(2) 



0093 

0084 

14-14*100 

0IN(2)-I4/ICNT 



0099 

0086 

I4-0IN(3) 

14-14*100 



0087 

D 

0IN(3)-I4/ICNT 
TYPE 991. FLAG 


• 

0088 991 

FORMAT (' FLAGM9I3) 



0089 

0090 

RETURN 

END 







• 
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FORffiAfl IV-PU'S V02-51E 14:40-07 3O-OCT-B0 PACE 4 

SCR.il.L FTN /TR DLOCKS/UR 


UKIG1N/\L I*. . - 

OF i^..Anrv 
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LISTING OF SKIP PROCESSOR 



B-1 


Of POOR QUALITY 


1 FORTRAN 

IV- 

•PLUS V02-S1E 14:30:53 30-0CT-80 PAGE 2 

8kT?TTTR 

/fR. duSSkITwR 

002S 

read ( 5. lOOd>PILCN 
1000 FORMAT (4A1) 


C 

CREATING T iE NAME OF THE FILES TO OE WRITTEN 

0024 

T“ 

00 IS K-1.6 

00^7 

0020 


FLNIMH(K)-FiLEN(K) 

FLNIMDCKlxFILENCK) 

555? 

c 

15 CONTINUE 


c 

c 

READ START LINE AND START PIXEL 

SSso 

0031 

READ<5. IIOOIUXNST 
1100 FORMAT! 112) 

0032 

C 

READ(5> 1100)IPIX8T 


t 

c 

THIS RbUTINE WILL REaD THE hEa6£r THE IMAOE ^aAE 

0633 

0034 


CALL <bE'fAbR<HbRA6R« 
CALL QI0(“2400iLUNl) 

0035 

c 



T“ 

c 

CHEC:;iN0 FOR ERRORS DURING TAPE READ 
CHECKING FOR END OF FILE 

0036 

T“ 

IF(IE. NE. 0)G0T0 40 


TT" 

c 

WRITE THE HEADER -IN IMDACS FORMAT 

0037 

T" 

CALL HDRFLN(LUN2, FLNIMH. UIFHR. SUNEL. FLAT. FLONG) 


T“ 

c 

OPEN THE IMAGE DATA FILE 

0038 

"5“ 

0PEK(UNIT-LUN2. NAME-FLNIMD, TYPE- 'NEW'. FORM- 'UNFORMATTED 


c 

* ACCESS- 'DIRECT'. RECORDSIZE-640) 


"c 

c 

CALL SKPRD PROGRAM T^'CE FOR EACH CCT 6TRIP 

0039 

0040 


DO 100 ICCT-1.4 

CALL 8KPRD<LUN1. LUN2. OUF. OFSZ. EOF. IE. SUFFER. IFRM. HDRADR. ICCT. 

0041 


1 LINST. IPIXST) 
IF<IE. NE. 0)G0 TO 40 

55?5 

0043 


IF (ICCT. EG. 4)00 TO 50 
IF( ICCT. NE. 2)00 TO 60 

C 

C 


c 

c 

INFORM THE OPERATOR TO MOUNT THE NEXT TAPE IF IT EXISTS 

0044 

c 

CALL QIO< "2400. LUND 

0045 

0046 

35 

WRITE<3. 6) 
READ <3. 7)ANS 

0047 

004S 


IF(AMS. NE. 'C')00T0 35 
CALL Q10(’‘2400. LUNl ) 


"c 

READ HEADER RECORDS FOR NEXT STRIP 


B-2 



ORIG!?^AL PAGE IS 
OF POOR QUALITY 


FORTRAN 

1V-FLU8 


SKXF. FTN 

/TR; BLOCK!/ WR 

\ . 

e 


004? 

40 

CALL HDREAD(LUN1. IFRH> BUFFER. EOF. tNft. H&f^AbA» 

1 188) 

5555 


1F(IE. NE. 0)00 TO 40 

0051 

K K*X*mHlii 



c 

JS 

error . 


C 


0052 40 CONTINUE 


0053 WRITE(6. 4>IE. NAM£ 

0054 QQTD >499 


c 

£ 

0055 


50 CONTINUE 


-iL 



c 

CLOSE THE IMAOE DATA FILE 


—5. 


0054 


CLOSE < UNI T«LUN2) 



CALL 0J0< •• 2400. LUND . 


c 




WRITE THE HISTOGRAMS TO THE HEADER FILE 


c 


0058 


CALL HDRHIS<LUN2.FLNIMH.UIFHR) 




0059 

9999 STOP 

55S5 


EnB 
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FORTRAN IV-PLUS V02-51E 

14 31 11 

30-0CT-B0 

PACE i 


SKPRD. PTN /TR; BLCCKS/UR 


oo(Jr susi^ouTiNE skpr6(luni. LUN&. DUF. BFsZi IE. buFFer. Ifru. H^ ftA&kT 


1 ICCT. UNST. IPIXST) 



C 

c 

THIS SUBROUTINE MILL READ THE LINES OF DATA FROM THE LAND8AT TAPE 


c 

c 

AND Select tHE pixels Td urite to the FIle. 


t 

dALL AR(^UMENT8: 


c 

LUNl LOGICAL UNIT OF TAPE DRIVE 


c 

LUN2 LOGICAL UNIT OF IMAGE DATA FIlS 


c 

BUF BUFFER FOR DATA RETURNED FROM LDCODE 


c 

BFSZ SIZE OF BUF 


c 

EOF END OF FILE INDICATOR FROM LDCODE »1 FOR EOF 


c 

IE ERROR INDICATOR -1 FOR IkRROR 


c 

BUFFER BUFFER FOR LDCODE TO READ TAPE RECORD INTO 


c 

IFRM FORMAT OF TAPE »2 FOR LANDSAT 


c 

HDRADR ADDRESS OF BUFFER 


c 

ICCT STRIP NUMBER CURRENTLY PROCESSING 


c 

LINST START LINE FDR SKIPPED IMAGE 


c 

c 

IPIXST START PIXEL FOR SKIPPED IMAGE 

0002 


IMPLICIT INTEGER (A-Z) 

0003 


REAL FLAT. FLONO. NLAT. NLONO. H4. H5. H6. H7 

0004 


PYTE BUF(l). BUFFER (l).SCNID 

OOOS 


DIMENSION LAD0(2). IPADD<4) 

o5o6 


COMMON /HIST/H4(256). H5(256). H6(256).H7(256) 

0007 


COMMON/COMH/ST. SE. LS. LE. NOCHAN. NOSAMP. SCNIDl 12). 

1 SUNELi SUNAZ. SNSHDi HREV. EXYR. EXOAY> OAV, HR. nZN< 

2 SEC. SEQN0<2). FLAT. FLONO. NLAT, NLONO 

COOS 


DATA IPADD/6. 6. 6. 7/ 

0009 


DATA LADD/4. 5/ 

5CT5 



0011 


EOF-0 



-g-g ^ 

0013 


ILINE-LINST 

CoT5 


LS^«1 

0015 


LSE-SE-ST+1 ' 



•. ' 

0017 


IFCICCT. NE. 1)00 TO 10 

SoTS 


IPIXl-IPIXST 

0019 


ISTRT-1 


c 

c 

DO LOOP FOR READING LINES FROM TAPE 

0020 

c 

10 

DO 100 I>1.256 

0021 


DO 100 11-1.2 

0022 


LINE-LINE+1 

0023 


DO 80 lCHAN-1.4 

0024 


CALL LDCODE (LUNl. IFRM. BUFFER. BFSZ. ILINE. ICHAN. LSS. LSE. 



1 BUF. EOF. PC. IE. INIT) 

0025 


IF (EOF. EG. 1 >G0 TO 900 

0026 


IF(1E. EG. 0)GO TO 20 

0027 


WRITE(6. 1200) ILINE. ICCT 

0028 

1200 

FORMATdX. 'ERROR ON TAPE READ OF DATA FOR LINE'. 15. 
1 ' FOR CCT STRIP '. 13) 

0029 


GO TO 90 

0030 

20 

INDX-0 
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1 FORTRAN 

IV-PLU8 

V02-51E 14:31:11 

30-0CT-80 

PAGE 3 

SKPRD. FTN 

/TR:OLOCK8/UR 



n 

o 

o 

c 

ipix*ipixi 




C 

c 

DO LOOP TO SELECT PIXELS 

FRON LINE 


0033 

0033 


00 65 J"l«33 
DO 60 NPIX-ISTRT. 4 




0034 INDX«INDX-<>1 

0035 BUF( IIMPX)»BUF< IP IX) 


0036 1P1X-IPIX^IPA00<NPIX) 

0037 IFdPlX. QT. LSE)QO TO 70 


0038 

0039 

60 

CONTINUE 

ISTRT-1 


0040 

65 

C 

CONTINUE 



C 

C 

MOVED ALL DATA TO BE USED FROM THIS LINE AND CHANNEL 

0041 

0043 

70 

80 

CALL SKPURT(LUNS«BUF< INOX* ICHAN. ICCT*LINE>. 
CONTINUE 



C 

C 

INCREMENT LINE NUMBER 


0043 

C 

90 

ZLlNE»ILINE-cLADO( II ) 

• 

0044 

100 

C 

CONTINUE 



C 

C 

FINISHED ALL DATA FOR THIS STRIP 


0045 

0046 


IPIXl-IPIX-LSE 

ISTRT-NPIX+1 


0047 

0048 


IFdSTRT. OT. 4)ISTRT-4 
ILINE-3000 


504? 

0050 


ICHAN*1 

CALL LOCODE(LUN1. IFRM* BUFFER. BFSZ. ILINE. ICHAN. LSS 

. LSE. BUF. 

0051 


1 EOF. PC. IE. INIT) 

IF ( EOF. NE. 1)00 TO 910 


0053 

C 

RETURN 



C 

C 

ERROR ENCOUNTERED EOF IN DAYa SET 


0053 

0054 

900 

1000 

URITIE(6. 1000) I CCT. ILINE 

FORMATdX. 'ERROR: END OF FILE IN STRIP'. 15. ' AT 

Lir.-E'. 15) 

0055 

0056 


IE"1 

RETURN 



C 

C 

ERROR DID NOT FIND END OF FILE AT END OF DATA 


0057 

C 

910 

IE-1 


0058 

0059 

1100 

WRITE(6. 1100) ZCCT 

FORMATdX. 'ERROR; NO END OF FILE FOUND FOR STRIP' 

. 15) 

0060 

0061 


RETURN 

END 
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SI 


II 
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8KPRD. FTN /TR: BLOCKS/WR 



TT17: C320i 403A. B/»SP»SYO: C320> 403SKPRD 









ORIGINAL PAGE IS 
OF POOR QUAimr 


FORTRAN IV-PLUS V02-S1E 

14; 31. 31 

30-0CT-30 

PAGE 1 


SKPWRT. 7tRTDLO<tkS/Wrt 




Wdi” suBftOuTtKte ' SKfiMyr^uui^ts, Bi^ 7T Ki g yr! ’ gT<m !g eT;Lii><e^ 


C 


— c 
c 

this routine will fill Th6 Output bopfer and waite yhc 

DATA TO THE IMAGE FILE FOR THE SKIP IMAGE PROCESSOR. 


t 

c 

CALL ARGUMENTS; 


c 

c 

ETun5 logical unIt for thE ^mag^ 

BUF BUFFER CONTAINING DATA TO INPUT TO BUFFER 


c 

c 

INDX COUNT OF NUMBER OF ENTRIES IN BUF TO PRO^E^? 

I CHAN CHANNEL NUMBER TO PROCESS 


c 

c 

ICCT STRIP NUMBER 

LINE LINE NUMBER TO WRITE 

0002 

c 

IMPLICIT INTEGER (A-Z) 

0003 

0004 


BYTE BUF(1).0BUF<2960).BW0R0<2) 
EQUIVALENCE<BWORD. IWQRD) 


c 

c 

CHECK TO SEE IF FIRST CCT SO INITIALIZE ANCILLARY DATA 

0009 

c 

IF<ICCT. EG. 4>INDX1-S12-ICNT 

0006 

0007 


IF<ICCT. EG. 4. AND. INDXl. LT. INDX ) INDX-INDX 1 
IFdCCT. NE. DGO TO 10 

0008 

0009 


IrdCHAN. NE. 1)00 TO 20 
ICNT-0 


MHI oouk<ii)-i 

0011 IWORD-LINE 


■S5T2 dQUF(6<^>«OWOkD<i) 

0013 0BUF(70)-QM0RD(1) 


0614 

C 

GO TO ^0 


C 

C 

11^ NflT CCT ^TRiP 1 ftEAD £)ATA WRltT^N ^OR 1 

53T3 

10 

C 

ff^dCHAN. &G. i)R£Ab(LUNi!'LlN£)OOUF 


C 

C 

DO LOOP FOR INSERTING DATA INTO PROfiei^ PART OF 6BU<*' 

0016 

0017 

26 

06 100 I-l. INDX 

IND«70+< ICHAN-1)*912+ICNT-H 

0018 

0019 

100 

OBUF(IND)-BUFd ) 
CONTINUE 


6 

c 

IF CHANNEL IS 4 WRITE DATA BACK TO FILE 

0020 

6 

IFdCHAN. EG. 4>WRITE{LUN2'LINE)0BUF 


6 

c 

IF LINE IS LAST LINE UPDATE ICNT 

0021 

6 

1F(LINE. EG. 512. AND. ICHAN, EG 4 ) ICNT-ICNT-dNDX 


c 

c 

IF LAST CCT AND CHAN 4 COMPUTE HISTOGRAM 

0022 

c 

IFdCHAN. EG. 4. AND. ICCT, EG. 4)CALL SKPHST<OBUF> 

0023 

0024 


RETURN 

END 
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1 FORTRAN 

IV-PLUS 

V02-SIE 14:31;S2 30-0CT-B0 

PAGE 1 

SKPHST, FTN 



0^0 i 

C 

'SkPH5T<0SUPI 



c 

c 




c 

c 

CALL ARGUMENTS; 



c 

c 

bUPPEA CilNTAllMlt«4A P1)(BU~5ATA 

TO COUWT 

0002 

0003 


6vtE 0DUF<1> 

C0MM0N/HIST/H4(2S6). H5(296). H6<296). H7(296) 



c 

c 

ADD TO PROPER HISTOGRAM 


0004 

t 

DO 100 1-1,912 


53os 

0006 


INDX1-0BUF(70+I). AN6. '*S77 
INDX2-0BUF<982+1).AND. *'377 


5557 

OOOB 


1nDX3-0BUP< 1094+1). AND. '*377 
IN0X4-0BUF( 1606+1 >. AND. “377 


5559 

0010 


H4( INDXl+1 )-H4(lNDkl+i)+i 
H9( INDX2+1 )-H5(INDX2+i )+l 


55TI 

0012 


H6( IN0X3+1 )-H6( INDX:!3+1 )+i 
H7( INDX4+1 )-H7( INDX4+1 )+l 

• 

55T3 

0014 

100 

CONTINUE 

RETURN 


55T3 


InB 


^ 
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APPENDIX C 


LISTING OF READTP PROCESSOR 



ORIGINAL PAGE » 
or POOR QUALITY 


fortran IV-PLUS V02-51E 14:2021 


READTP, FTN /TR: L..OCKS/UR 


30-0CT-B0 


PAGE 1 


r.T*T*TI 


0002 


0003 


0004 


0005 

0006 


0007 

0008 


0009 

0010 


C 

C THIS PROGRAM WILL READ A TAPE IN LANDSATI 1. 2)i UNIVERSAL. 
C AND LARSVfi FORMAT. WHILE THE TAPE IS READ. THE DATA 18 


XTR ACTED IN ' 

C SIXTEEN PATCHES PER FULL FRAME LINE. OR A TOTAL OF 206 




NOMINAL POINT. 



r^5!B:W TMiTJrJ ?i?4Tr^1 


C LOCAL VARIADLES; 


C OFSZ- DUFFER SIZE IN BYTES 


BUF> IMAGE DATA QU 
C DELTAL- THE CHANGE IN CENTER POINT LINE 


E 

C IE- ERROR FLAG 


INCL- STARTING LINE INCREMENT FOR 
INCP- STARTING PIXEL INCREMENT FOR .EACH ROW OF PATCHES 


C INITL- STARTING LINE NUMBER TO 
C INITP- !^TARTING PIXEL NUMBER TO ^E USED FOR EXTRACTINw IMAGERY DATA 


C NOML- NOMINAL CENTER POINT LINE 

J NOMP - NOMINAL C ENTER POINT PIXEL 

C NPOINT- (ARRAY) NOMINAL CENTER POINT IN DEO. 

C PSIZE- THE NUMBER OF PIXELS ON A SIDE OF PATCH (32.64.126.296) 


C STARTL- THE LINE TO START EXTRACTING FROM 
C 8TARTP- THE FIRST PIXEL TO EXTRACT DURING PATCH 


EXTRACTION FROM FULL FRAME 
C 8T0PP- THE LAST PIXEL TO EXTRACT DURING PATCH 


EXTRACTION FROM FULL FRAME 
C 8TRIPC- THE VALUE "SEQUENCE NUMBER" OF CCT 


IMPLICIT INTEGER (A-Z) 


REAL FLAT. FLONG. NLAT. NLONO. DELTAL. DEL TAP. P 
♦ H4, H9. H6. H7 


LOGICAL*! BUFFER ( 32B0 ). SCNID. BW0RD(2) . PIXBUF(21 18) . ANS. 

* *■ FILEN(6).FLNIMD(16).FLNIMH(i6).MISNO(2).BUF(4000) 


LOGICAL*! BEL1.(4>. ISTATB(4) 

DIMENSION NP0INT(2). IPRM(6) . 0CT(2). JIFHR( 1936) . ISTAT(2) 


EQUIVALENCE (OCTAL. B WORD) . ( ISTATB. ISTAT) 




SUNAZ. SNSHD. NREV. EXYR. EXDAY DAY. HR. MIN. 
SEC. SEGNO ( 2 ) > FLAT. FLONG. NLA T. NLONO 


COMMON/LADEL/NAMEO) . MISNO. £ 

COMMON/H I ST/H4 ( 296 ) . H9 ( 296 ) . H6 ( 296 ) . H7 ( 296 ) 


bATA UIFHR/1. 912. 1. 912. 1. 912. 1. 9l2. 3. 4. 1.2. 1. 1 


0012 


<•19729. 19861. -19729. 19861. -19729. 
-J* 19861.-19729.19961.-15729.1.427.1. 


* 3*0.4.297. 297. 10*0. 512. 1099. 1.2. 1459*0/ 

DATA BFSZ/4000/. WLB/"0400/.LUNl/l/.LUN2/2/, INIT/0/. 


* PIXBUF/2118*0/. B£LL/"007. "007, "007, "007/ 

DATA FLNIMH/' ' . ' , ' I '. 'M', 'H 


♦ * i * » *»**$**» 0/ » , 

* FLNIMD/" "I ", "M". "D ' , 


C -1 



ORIGINAL PAGE IS 
OF POOR QUALITY 



C FOMnAT LIITINO; 


0015 


1 F0RHAT<17X, 15. f, 18X. 15, /. 12X. 13. /. lOX. AAl, /, 11X713; 

9X.14. /. 14X. !4. /, SOX. 14) 


4 FORMAK/. IX. 10( ERROR- M5. 'AT- 
♦ 3AS. /. IX. 10( '*') 


5 F0RMAT<I3> 

FORMAT( IX. 10< /. TX. 4A1. '* REPLACE MOUNTED TAPE WITH NEXT 


'TAPE IN SEQUENCE'. /. IX. tO( '»')) 
a FORMAT! IX. ' » DAD PATCH SIZE') 


NITIALIZINO HISTOGRAMS TO 


% OOSO 




0024 


H4(L)-0 




H6<L)-0 
H7<D- 


10 CONTINUE 


8TRIPC-0 



READING THE NOMINAL POINT TO BE USED IN EXTRACTION OF DATA AS THE 
REFERENCE POINT FROM THE PARAMETER FILE. 


0PEN<U!)IT-3. NAME- 'RPARAM. DAT '. TYPE- 'OLD ' ) 


READO. 1. ERR-45, END-9999 )NP 01 NT. UERSN. FILEN. PSIZE. 

INCL. INITL, INCP. INITP 


CLOSE (UN IT-3) 


CHECK FOR PATCH SIZE TO BE USED FOR PATCH IMAGE CONSTRUCTION 


IFIPSIZE. EQ. 32. OR. PSIZE. EQ. 64)G0T0 12 
IF(PSI2E. EQ..12S. OR. PSIZE. EQ. 256)G0TQ 12 


WRirE(6. 8) 
GOTO 9999 . 


12 CONTINUE 



0035 

0034 

0037 


C CREATING THE NAME OF THE FILES TO BE WRITTEN 
C 

DO 13 K-l.A 

FLNIMH<K)-FILEN(K) 

FLNIMD(K)-FILEN<K) 

IS CONTINUE 


C THIS ROUTINE WILL READ THE HEADER OF THE IMAGE TAPE 

s 

0038 CALL GETADR(HDR ADR. BUFFER) 

0039 CALL QIO<”240Q. LUN l > 

0040 CALL HDREAD(LUNI. ZFRM. BUFFER. BFSZ. EOF. PC. IE, BUFFER. INS. HDRADR. ISS) 

s 

C CHECKING FOR ERRORS DURING TAPE READ 

C CHECKING FOR END OF FILE 


C-2 
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FORTRAN 

IV-PLU6 U02-91E 14: SO' 81 30-0CT-80 PAGE 3 

RiAfitF. FTN 

/TR; BLOCiaVUR 

0041 


IF<IC. NC. 0)00TQ 40 


IT 

C 

BUILD THE HEADER FILE NAME. AND THE IMAGE DATA FILC NAME 

0048 

T' 

CALL INTOCT(EXDAV. OCTAL) 

0043 

0044 


FLNIMH(12)>QU0Rb(i) 

FLNIMH<13)-OUORD(1) 

0049 

0046 


FLNIMDi 12)«0W0RD(3) 
FLNIMD<13)-BU0RD<1) 

0047 

004B 


CALL INTOCT ( VERSN. OCTAL > 
FLNIMH(14)-DW0RD<2) . 

0049 
. 0090 


FL.>ilMH(19)«QM0R0(l) 
FLNIMD( 14)-BH0RD<2) 

0091 

c 

FLNIM0<19)-DW0RD(1) 


c 

c 

WRITE THE HEADER IN IMDAC8 FORMAT 

• 

0092 


CALL H0RFLN(LUN2. FLNIMH. UIFHR. 8UNEL. FLAt. FLONb) 

- 

c 

c 

OPEN THE IMAGE DATA FILE 

0093 


0PEN(UNIT-LUN2. NAME-FLNIMO. TYPE- 'NEW'. WRM- 'UNFORMATTED 
* ACCESS- 'DIRECT '. RECORDSlZE-640) 


t 

c 

FILL THE ENTIRE FILE WITH ZEROS 

0094 

.c 

DvT i M-1. 912 

0093 
' 0096 


WRl rt:(LUN2'M)PIXBUF 
16 CONTINUE 


"2~ 

c 

DETERMINE IF THE IMAGE TAPE IS DIFFERENT FROM THE HDT FULL FRAME 

0097 

c 

IF<IN8. NE. 3948)G0TP 30 

0098 

c 

20 CONTINUE 

0099 

c 

STRIPC-STRIPC+1 


c 

c 

THIS ROUTINE WILL BUILD A PATCH FRAME IMAGE USING HDT FULL FRAME 
SIZE IMAGES 

0060 

c 

CALL HDT(LUN1. LUN2, DUF. BFSZ. NPOINT. FLNIMD. EOF. IE. 

f BUFFER. IKRM, ISS. !NIT, STRIPC. PSIZE. INCL. INITL. 

♦ INCP, INITP) 

0061 

0068 


IP<IE. NE. 0>G0T0 40 
IF<STRIP':. EO. 4)00T0 90 

0063 


IF<EOF. EQ. DCALL HDREAD(LUN1. IFRM. Suf^F^R. BFSZ. fiOF. 

PC. IE. DUFFER. INS. HDRADR. ISS) 

0064 

0069 


IF< IE. NE. 0)G0T0 40 
ASSIGN 20 TO PROCES 

0066 

0067 


IF < EOF. EO. 1 )GOTO 39 
INIT-0 

0068 

c 

GO 10 20 

0069 

0070 

. 

30 CONTINUE 

STRIPC-STRIPC+1 
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FORTRAN 1V-PLU8 V02-9II 

■«A6t#f?tN 7tR7iL( 


14 aO;gi 30-0CT~e0 


THIS RO U TINE MILL BUILD A PATCH FR/ 
1 fcgFULL FR^E SIZf iMAoS 

CALL LAN8AT<LUN1 . LUN2< DUF. R^IZ. 1 


INCP, INITP) 
IF<IC.NE.O)OOTO 40 
IF<8TRIPC. EO. 4>00T0 SO 


DUFFER. INS. HORAOR. 188) 


IF( IE. NE. 0)C0T0 40 
ASSIGN 30 TO PR0CE8 
I F<eQF^EO, UGOTO 3S. 
INIT-0 

GOTO 30 


CALL HTQIO( ‘'2340. LUNl. 1. . ISTAT. AST) 

IFdSTATDd). LT. 0)OOTO 40 

MRITE(4. 6)DELL 

37 CONTINUE 

CHECK TO SEE IF DATA TAPE IS READY 

CALL UTQIO( ”5520. LUNl. 1, . ISTAT. IDS) 

IF< 1STATD(4). NE. "044>00T0 Z7 

CALL HAIT <18000. 1.11) 

INIT-0 

CALL HDREAD(L<;N1, IFRM, BUFFER. BF5Z. EOF. PC. IE. 
* BUFFER. INS. HDRADR. ISS> 

GOTO PROCES 

40 CONTINUE 

klRITE(6. 4)IE. NAME 

GOTO 9999 ] 

45 CONTINU E 

GOTO 99?9 

50 CONTINUE 

CLOSE THE IMAGE DATA FILE 
CL0SE<UNIT-I.UN2) 

CALL QIO< "2400, LUNi > 


C HRITE THE HISTOGRAMS TO THE HEADER FILl 


_F0RTRAN IV-PLUS 1 E 

READTP. FTN >^TR: QLOCKS/UR 


fS 

QUAUty 


_i4: 20; gl 3Q-QCT~BQ 
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1 FORTRAN IV-PLU5 VOa-SlE- 14:21:00 3Q»0CT-fl0 PAGE 1 



lansat. ftn 

/TR: BLOCKS/UR 


0001 


SUBROUTINE LANSAT<LUNl,LUN2,bOF. BFSZ, NPOINT, i^LNIMD. e6f, IE. 

IFRM, SS. INIT, STRIPC. PSIZE, INCL. INITL, INCP. 




- ♦ INI TP) 



c 

THIS ROUTINE HILL EXTRACT PATCH IMAGES FROM THE FULL FRAME IMAOEt 
PROCESSING TAKING PLACE ONE STRIP OF CCT AT A TIME. IT WILL 



c 

FURTHER BUILD A PATCH FRAME IMAGE 16 PATCHES BY 16 PATCHES. EACH 
PATCH CONSISTING OF 32 BY 32 PIXELS. BANDS C4-73 . THE ROUTINE 



c 

WILL USE LANOSAT 1 AND 2 FULL FRAME SIZE IMAGES. 



c 

LOCAL VARIABLES: 

•• 


• 

c 

c 

BLOCKN- THE BLOCK NUMBER 

ENDCCT- A PARAMETER FOR READING THE END OF THE FILE 



c 

c 

LCOUNT- THE LINE COUNTER FOR THE PATCH IMAGE 
LEND- THE ENDING LINE WITHIN THE FULL FRAME 



c 

c 

LSTART- THE STARTING LINE WITHIN THE FULL FRAME 
NPATCH- THE NUMBER OF PATCHES ON A SIDE OF A PATCH. IMAGE 



c 

c 

PATCH- THE CURRENT PATCH NUMBER 
PIXELS- THE PIXEL LOCATION IN PIXBUF 



c 

c 

PTCH- A FLAG TO INDICATE THAT THE PATCH COUNTER HAS 
SAVEP- THE LAST PATCH PROCESSED ON PREVIOUS CCT 



c 

c 

STARTP- THE STARTING PIXEL LOCATION IN ''ATCH EXTRACTION 
STOPP- THE LAST PIXEL LOCATION IN PATCH EXTRACTION 



c 

c 

BEEN READJUSTED 


0002 

0003 


IMPLICIT INTEGER <A-Z) 
LOGICAL NEXTL, NEXTP, PTCH 


0004 

0009 


REAL FLAT, FLONG, NLAT, NLONO, DELTAL, DEL TAP, PO 
LOGICAL*! BUF( 1 ). PIXDUF(2118), FLNIMD(14),SCNID, MISN0<2) 


0006 

0007 


LOGICAL*! BUFFER! 1), BLOCKS (2) 
DIMENSION NP0INT(2), ISTAT(2) 


. oooa 

0009 


EQUIVALENCE (BLOCKS, BLOCKN) 

■ COMMON/COMH/ST, SE, LS, LE. NOCHAN, NOSAMP, SCNIDC 12), 

SUNEL. SUNA2, SNSHD. h4r“EV, EXYR. EXDAV, DAY, 

HR, MIN, SEC. SEQfJ0(2). FLAT, FLQNO. NLAT. NLONO 


0010 

0011 


COMMQN/LABEL/NAME ( 3 ) , M I SNO, E 
DATA SAVEP/0/, T3E/0/, T3T/0/ 



c 

c 

INITIALIZING PARAMETERS 


0012 

c 

LSS-1 


0013 

0014 


IFdFRM. NE. 3)TSE»SE 
IF< IFRM. EQ. 3)TST-TSE+1 


0019 

0016 


IF! IFRM. EQ. 3 ) TSE-TSE+SE 
IFdFRM. NE. 3)TST«ST 


0017 

0010 


LSE*SE-ST*1 

EOF-0 

2 

t 

0019 

0020 


DO 5 M-71, 2110 
PIXDUF(M)-0 

0 

• 

0021 

0022 


9 CONTINUE 

NPATCH-912/PSIZE 

1 

7 


c 

c 

CHECKING FOR THE CASE WHERE ALL THE PATCHES HAVE BEEN 

r 


c 

EXTRACTED 


C 
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FORTRAN 

XV-PLUS V02-91E 14:21:00 30-0CT-80 PAGE 2 

UANSAT. FTN 

/TR: BLOCKS/UR 

5553 

0024 


IF<SAVEP. NPATCH)GOTd 99^4 

PCOUNT-SAVEP 


C 

SET START AND STOP PIXEL FOR READING EACH STRIP PROM CCT 


C 

C 

FST START LINit FOR EXTRACTION FROH CCT 


C 

C 

PO— 1*<. 299*0. 00363»FLAT) 
OELTAL«(NPOXNT( 1 )-FLAT)»23. 073 


C 

C 

OELTAP- ( NPOINT ( 2 ) -FUONO ) *COS ( FLAT > *32. 99 
NOMP-1919, ♦DELTAP+DELTAL*PO 

0029 

C 

NOML-1170. ♦DELTAL 
URITE<6. IIIDDELTAL. DELTAP 

002d 

• 

1111 FORMAK IX. 'DELTAL-'. IFIO. 9. 'DELTAP*'. IFIO. 9> 
C 

0027 

0028 


DO 90 M-l.NPATCH 
LCOUNT-0 

5559 

0030 


LINE«n*PSIZE > 

STAR TL« I N I TL+ I NC L» < M- 1 ) 


“5 

c 

THIS ROUTINE WILL CHECK FOR THE CASE WHERE THE STARTING 


c 

c 

AND ENDING LINE OF PATCHES ARE OUT OF THE , FULL FSaHE 

5531 

0032 


CACl LOORDR (SYaRTL. NEXTL. LS. LE. LENS. LSTAf^t; PSl 2E) 
IF(NEXTL)OOTO SO 


T” 

c 

PROCESS LINES FOR EACH PATCH 

0033 

"5“ 

DO 70 K-LSTART, LEND 

553i 

' 0039 


L50UNTiLCQUNT-*-l 

BLOCKN-L I NE-PS I ZE+ ( LSTART- 1 ) •►LCOUNT 


T~ 

c 

PROCESSING ONE BAND AT A TIME FOR EACH LINE IN PATCH 

>0036 

"5“ 

DO 60 1-1.4 

0037 

0038 


SAVEP -PATCH-1 
PATCH-PCOUNT 

5559 

0040 


RCHAN-I 
PTCH-. TRUE. 


c 

c 

THIS ROUTINE WILL READ ONE LINE OF IMAGE DATA ONE 


~r~ 

c 

CHANNEL AT A TIME 

5541 


CALL LDCODECLUNl. IFRM. BUFFER. BFSZ. STARTL. RCHAN. LSS. 
♦ LSE. BUF. EOF, PC. IE. INIT) 


c 

c 

CHECKING FOR READ ERROR 

0042 

c 

1F<E0F. EQ. DIE-l 

0043 

C 

IF< IE. NE. 0)G0T0 9999 

0044 

0049. 


10 CONTINUE 

PATCH-PATCH* 1 


c 

c 

CALCULATING THE NUMBER OF PIXELS PROCESSED PER PATCH 


c 

c 

CALCULATING START A.NO STOP PIXELS PER PATCH 
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FORTRAN IV-PLU8 V02-916 . t4:21:00 30-0CT-B0 RACE 3 


-r^MSAT 

FTN 


/TR. DLOCKS/uM 

0046 

0047 



P 1 X^LS- < PATCH- 1 > *(^S I iE* ( 70+9 1 4* < t - 1 ) ) «•*- 

STARTP«lNlTPfINCP*<PATCH-l ) 

004a 

0049 


- 

IF<STARTP. OT. TSE. dR. PATCH. OT. NPATCH>OOTO 6^ 
8T0PP-STARTP4.P9XZE-1 ' 


C 


THIS ROUTINE WILL CHECK THE DOUNDARV OF THE PATCH ... 


c 


TO DETERMINE IF XT'S IN THE FULL FRAME 

r090 

0091 



CALL PB0RDR<STaRTP. STOPP. STOP. NEXTP, TSE. tSt, PSIZE) 
IF(NEXTP)COTO 10 


c 

c 


THIS ROUTINE MILL CHECK FOR THE CASE WHERE A PATCH 

• 

«F 

c 


IS DIVIDED OY TWO CCT'S 

0092 



CALL PARTP<PTCH, STARTP, STOPP. STOP. STRIPC. TSE. T8T. 
PIXELS. I. PATCH. PSIZE) 


c 

c 


READ PIXEL ARRAY FROM . IMD 'FILE 

0093 

c 


RE AD ( 2 ' DLOCKN ) P I DUF 


c 

c 


PLACE THE DESIRED PIXELS IN THE PIXEL ARRAY 

0094 

c 


L-0 

0099 
' ■ 0096 



DO 90 J-STAR TP. STOPP 

L-L>1 - . 

0097 

0098 



N-STARTP-TST+L 

PIXCUF<PIXELS-i-L)>GUF(N) 

. 0099 

c 

90 

CONTINUE 


c~ 

c 


PROCESS DATA EXTRACTED FOR HISTOGRAM 

0060 

c 


CALL HIST<DUF. STARTP. STOPP. RCHAN. TST) 


c 

c 


WRITE THE PIXEL ARRAY BACK TO . IMD FILE 

0061 

0062 



PIXBUF(2>»1 
P I XDUF < 69 > -BLOCKS ( 2 > 

0063 

0064 



' PIXDUF<70)»BL0CKD(1) 

WRITE<2'BL0CKN)PIXBUF . 

0069 

c 


GOTO 10 

0066 

c 

60 

CONTINUE 

0067 

c 

70 

STARTL-STARTL^l 

0063 

c 


GOTO 90 

0069 

c 

80 

CONTINUE 

1 » 

c 

c 


WRITING ZERO FILLED PATCH LINES 

1 0070 

0071 



DO 85 J-71.2118 
PIXBUF< J)«0 

0072 

c 

85 

CONTINUE 
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FORTRAN 

IV-PLUS V02-5IE ' 14; 21. 00 30-0CT-80 PACE 4 

.ANSAT. FTN 

/TR; BLOCKS/UR 

. 

0074 


00 49 j-i.psize 

' - BL0CKN«BL0CKN4>1 

OO?? 

007& 


’P!XBUF(4)-1 ■> 

P I X BUF ( 69 ) XBLOCKO ( 2 ) 

5397 

0078 


P 1 X DUF < 70 r-BL0dHii“( 1 ) 
WRITE! LUN2 ' BLOCKN > P I XDUF 


C 

SJ cdNYiiviue 

5600 

C 

■i^5”55NtiNue 


C 

C 

f^EAD UNYiL ENlTbP CCT fetRtP ‘ 

0081 

0082 


gND55T=3555 

CALL LDCOOE(LUN1 . IFRh. BUFFER. OFSZ. ENDCCT, RCHAN. LSS. LSE. 


C 

? bufTeBPTpcTTITTniT) 

# 

5505 

0084 
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1 FORTRAN 

IV 

-PLUS V02-91E • 142198 30-0CT-B0 PACEl 

HOT. FTN 


/TR BLOCKS/ WR " 

0001 


SUBROUTINE HDT<LUNl . LUN2. BUF. BFSZ. FLNIhO. ^OF. IE. DUFFER. 

• ♦ 88. INIT.STRIPC.PSIZE. INCL. INITL, INCP. INITP» . 


c 

THIS ROUTINE WILL EXTRACT PATCH INAOES FROM THE FULL FRAME 


c 

JL_ 

PROCESSING TAKING PLACE ONE STRIP OF CCT AT A TIME IT WILL 

further DUltO A RATCH FRAME IMAGE 16 PATCHES BY 16 PATCHES. FACH 


c 

JL. 

PATCH CONSISTING OF 32 BY 32 PIXELS. BANOS C4-7J . THE ROUTINE 
WILL USE HOT FULL FRAME SIZE "MAGES. 


c 

LOCAL VARIABLES: 


c 

c 

Bt.OCKN- THE BLOCK NUMBER 

• 

•A 

c 

c 

ENDCCT- A parameter FOR READING THE END OF THE FIL^ 
LCOUNT- THE LIME COUNTER FOR THE PATCH IMAGE 


c 

c 

LEND- THE ENDING LlNt WITHIN THE FULL FRAME 
LSTART- THE STARTING LINE WITHIN THE FULL FRAME 


c 

c 

NPATCH- the NUMBER OF PATCHES ON A SIDE OF A PATCH IMAGE 
PATCH- THE CURRENT PATCH NUMBER 


c 

c 

PIXELS® THE PIXEL LOCATION IN PIXBUF 

PTCH- A FLAG TO INDICATE THAT THE PATCH COUNTER HAS 


c 

c 

SAVEP- THE LAST PATCH PROCESSED ON PREVIOUS CCT 
STARTP- THE STARTING PIXEL LOCATION IN PATCH EXTRACTION 


c 

c 

STOPP- THE LAST PIXEL LOCATION IN PATCH EXTRACTION 
BEEN READJUSTED 

' • 0002 

T" 

IMPLICIT INTEGER <A-Z) 

0003 

0004 


LOGICAL NEXTL. NEXTP. PTCH ’ . - 

REAL FLAT, FLONO, NLAT. NLONG 

’ 0009 
0004 


LOGICAL#! B'.F( 1 ). PIXBUF(21 18). FLNIMD< 14), 8CNID, MI^Nb'(2) 
LOGICAL#! DUFFER < 1). BLOCKS (2) 

0007 

0008 


DIMENSION ISTAT(2) 

EQUIVALENCE (BLOCKS. BLOCKN) 

.0009 


COMMON/COMH/ST. 3E. LS. LE. NOCHAN. NOSAMP. S(^NID< 12). 

■+ SUNEL. SUNAZ. SNSHD. NREV. EXYR, EXDAY, DAY. 

0010 


# HR. MIN. SEC. SEQNO',2). FLAT. FLONC. NLAT. NLONG 

COMMON/LA0EL/NAME(3). MISNO. E 

0011 

c 

DATA SAVEP/O/.TSE/O/.TST/O/ 


c 

c 

INITIALIZING PARAMETERS 

0012 

0013 


LSS-1 

IF( TFRM, NE, 3)TSE-SE 

0014 

0019 


IF< IFRM. EG. 3>TST«TSE+1 
IF< IFRM EG. 3)TSE-TSE+SE 

0016 

0017 


IF< IFRM. NE 3)TST»ST 
LSE-SE-ST+1 

0018 

0019 


EOF-0 

DO 9 M-71.2118 

0020 

0021 


PIXBUF<M)«0 
9 CONTINUE 

0022 

0023* 


NPATCH-M2/PSIZE 

PCOUNT-SAVEP 


c 

SET START AND STOP PIXEL FOR READING EACH STRIP FROM CCT 


c 

SET START LINE FOR EXTRACTION FROM CCT 

! 
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PORTRAN 

IV- 

-PLUS V03<'3tE- 14 31:38 30*0CT>80 PAGES 



/TR OLOCK^/UA 

ooss 


————————— 

5S2S 

0037 


LiNE-n*i»dt2e ; 

8TARTL«INITL-*-INCL*<«-l ) 


T~ 

c 

THIS ROUTINE UILL CHECK FOR THE CASE WHERE THE STARTING 


T~ 

c 

AND UNblNG LINE 6^ PAtCHE^ ARE buf OP YhE I^ULL/ I^AahC 

5o5i 

0039 




IT 

c 

PROCESS? LINES POP EACH PATCH 

• 

•&030 

T~ 

DO 70 K«LST ART. LEND 

553T 

0033 


LdbUNT-LtOUNt-^l 

DLOCKN-L INE-PS I ZE-^ < LSTART- 1 ) ♦LCOUNT 


c 

PROCESSING ON^iS SAND AT A TIME FOR EACH LINE IN PATCH 

0033 

IT 

DO 60 1-1.4 

55331 

0033 


SAUef»-t‘AUH-« 

PATCH-PCOUNT 

553i 

0037 


ITfiHAFPT 
PTCH-. TRUE. 


T" 

c 

THIS ROUTINE WILL READ ONE LINE OP IMAGE DATA ONE 


T~ 

c 

CHANNEL At A TTmE 

' 6539 




T" 

c 

CHECKING FOR READ ERROR 

■*0037 

c 

IF(EQF. EQ. 1)IE-*1 

0040 

c 

IF< IE, NE. O)0OTO 9999 

0041 

0043 


10 CONTINUE 

PATCH-PATCH^l 


c 

c 

CALCULATING THE NUMBER OF PIXELS PROCESSED PER PATCH 


c 

c 

CALCULATING START AND STOP PIXELS PER PATCH 

0043 

0044 


PIXELS-<BATCH-1 )*PSIZEf(70-*-513*< I-l ) > 
STARTP»1NITP>INCP*(PATCH-1 > 

0043 

0046 


IF<STARTP OT, TSE)GQTQ 60 
ST0PP-STARTP-*-PSIZE-1 


c 

c 

THIS ROUTINE WILL CHECK THE BOUNDARY OF THE PATCH 


c 

c 

TO DETERMINE IF IT'S IN THE FULL FRAME 

0047 

0048' 


CALL PBORDR<STARTP, STOPP. STOP. NEXTP. TEE. TST, PSIZE) 
IF(NEXTP)COTO 10 


c 

c 

THIS ROUTINE WILL CHECK FOR THE CASE WHERE A PATCH 


c 

c 

IS DIVIDED BY TWO CCT'S 
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i FORTRAN IV-PLU« VOa-SiC 14.21:59 30-0CT-e0 l»40g 3 


HOT FTN 



/TR . DLOCKft/WI 

004V 


i 



T" 

. 

READ PIXEL ARRAY FROH . IHO FILE 

0090 



✓ 

HEAD ( 2 ' BLOCKN ) A I XICF 


T" 


PLACE THE DESIRED PIXELS IN THE PIXEL ARRAY 

0091 



L«0 

5o52 

0093 




0094 

0099 



RStttARtF-TST*L 

PIXBUF(PIXELS>L>«DUF(N) 

0«94 

j£_ 


CONTINUE 


c 

c 


MOCESS 6aTA extracted for Hl^TOORAH 

0 

53?’? 

c 


CALL HlSf <8UF. STARTP, sVOPP. RdHAN/t6t> 


T“ 

c 


URiti THE PIXEL ARRAY SACK TO IMO VlL^ 

5o5S 

0099 



piXBUF(i)-l 

PIXOUF(69)-8LOCHD(2) 

5oS3 

0041 



Ft'xBuPrf6)»aL5ck6< i') 

UR I TE < 2 ' BLOCKN > P I XBUF 

'0042 

T" 


GOTO 10 

. 0063 

c 

60 

CONTINUE 

0064 

c 

70 

STARTL-STARTL^l 

006^ 

c 


GOTO 90 

0066 

c 

80 

CONTINUE 


c 

c 


WRITING ZERO FILLED PATCH LINES 

0067 

0068 



DO 89 J-71.2118 
PIXBUF<J)-0 . 

0069 

S. 

89 

CONTINUE 

0070 

0071 



DO 87 J-l.PSIZE 
OLOCKN-BLOCKN+l 

0078 

0073 



PIX8UF(2)«1 

PIXBUF(69>-8L0CKB<2) 

0074 

0079 



PIXBUF<70>-BL0CK8< 1 ) 
WRITE(LUN2'8L0CKN)PIXBUF 

0076 

c 

87 

CONTINUE 

0077 

c 

90 

CONTINUE 

• 

c 

c 

READ UNTIL END OF CCT STRIP 

0070 

0079 



ENDCCT-3000 

CALL LDC0DE(LUN1, IFRM, BUFFER. BFSZ. ENDCCT. RCHAN. LSS. LSE. 
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I 

I FORTRAN IV-PLUS VOa-SlE- "’»:i4:2a:19 30-0CT-80 PAQE X 


POOR OR 

FTN 


/TR. BLOCKS/ WR 

000 1 

c 


SUBROUTINE PBOROR (6VaRTP. stOPP. STOP. NEXfP. fS£. T^t. PSli^i 


c 

THIS ROUTINE WILL DETERMINE WHEN IMAGERY DATA DESIRED FROM 
THE FULL FRAME IN NOT AVAILADLEi PATCH PIXELS ARE OUTSIDE 


c 

-JL_ 

OF 

THE FULL FRAME 


c 

LOCAL VARIABLES; 


c 

NEXTP* A FLAG WHICH INDICATES) PROCESS THE NEXT PATCH 

OOOS 

0003 



IMPLICIT INTEGER <a»Z) 
LOGICAL NEXTP 

• 

r* 

c 

c 

CHECKING FOR THE CASE WHERE PATCH PIXELS ARE NOT WITHIN 


c 

c 

THE FULL FRAME 

0 

0004 

0005 



NEXTP*. FALSE. 

IF<STARTP, OE. TST. AND. STARTP. LE. (TSE-PSIZE-1 ) )GOTO 9999 

0006 

0007 



IFISTARTP.OT TST ) 0OTO"4O ' 

IF< (STARTP+PSIZE-1 >. OT. TST)GOTO 10 

oooa 

0009 



NEXTP*. TRUE. 
GOTO 9999 

0010 

c 

10 

continue 


■"c ' 
c 

CALCULATE THE NUMBER OF PATCH PIXELS WHXCH ARE OUT OF 


c 

c 

THE LEFT BORDER OF THE FULL FRAME 

' oou 
0012 



DO 20 K*1.PSIZE 

IF< (STARTP+K-l ). EQ. TST)00T0 30 

0013 

c 

20 

CONTINUE 

.0014 

0015 


30 

CONTINUE 

STARTP*TST 

0016 

0017 

; 


STOPP*STARTP-*-PSIZE-K 
NEXTP*, FALSE. 

0018 

c 


GOTO 9999 

0019 

c 

40 

CONTINUE 


c 

c 

CALCULATE THE NUMBER OF PATCH PIXELS WHICH ARE OUT OF 
THE RIGHT BORDER OF THE FULL FRAME 

0020 

c 


IF(9TARTP. CT. TSE)GOTO 60 

0021 

0022 



DO SO K-t.PSIZE 

IF( (STARTP-t-K-l). EQ, TSE)GOTO 70 

0023 

c 

30 

CONTINUE 

0024 

0025 


60 

CONTINUE 
NEXTP* TRUE. 

0026, 



GOTO 9999 

0027 

0020 


70 

CONTINUE 

STOP-STOPP 

0029 

2232 



STOPP-K+STARTP-1 
NEXTP*. FALSE. 
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0031 


9999 RETURN 
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t FORTRAN 

19- 

-PLUS U02-51E 14.22:37 30-0CT-B0 PAGE 1 

LOOROR. 

Ptn 

/TR: BLOCKS/UR 


C 



C 



C 

IS NOT AVAILABLE; PATCH LINES ARE OUTSIDE OF THE FULL FRAME. 


“3“ 



c 

LOCAL VARIABLES: 


"5“ 

c 

LEND- LINE END 


c 

LSTART- LINE START 


c 

NEXTL- A FLAG WHICH INDICATES; PROCESS THE NEXT FULL FRAME LINE 

0002 

c 

IMPLICIT INTEGER (A-Z) 

0003 

c 

LOGICAL NEXTL 


c 

c 

INITIALIZING PARAMETERS 

<5o04 


LEMD-PSIZE 

0009 


LSTART-1 


"3” 

c 

f 

CHECKING FOR THE CASE WHERE PATCH LINES ARE NOT WITHIN THE FULL FRAME 

0006 

IT 

NEXTL-. FALSE. 



IFCSTARTL. CE. LS. AND. StARTL. Lfe. <LE-PSlife-l ) )(iOTO 

0008 


IF(STARTL. OT. LSIOOTO 90 




0010 


NEXTL-. TRUE. 

66ii 

c 


mmiiii 

c 

20 CONTINUE ^ 


*3“ 

mm Ui.r — — — 


c 

TOP BORDER OF THE FULL FRAME 

0013 

~T~ 

DO 30 K-l.PSIZE 

' 0014 


. IF( (STARTL+K-l ). EQ. LSIOOTO 40 

0019 


30 CONTINUE 

0016 

c 

40 CONTINUE 

0017 


LSTART-K 

0019 


NEXTL-. FALSE. 

0019 

c 

GOTO 9999 

(5320 

c 

50 CdlNTINU^ 


c 

CALCULATE THE NUMBER OF PATCH LINES WHICH ARE OUT OF THE 


c 

BOTTOM BORDER OF THE FULL FRAME 

0021 

c 

IF { STARTL. GT. LE) GOTO 70 

0022 


DO 60 K-l.PSIZE 

0023 


IF( (STARTL-^K-l ). EQ. LElGOTO 80 

> 0024 

c 

60 CONTINUE 

> 002» 


70 CONTINUE 

' 0026 


NEXTL-. TRUE. 

1 0027 


GOTO 9999 


* c 
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PAGE 2 

LDOROR. FTN /TR: Bt.0CK8/WR 

OOSB 

0059 

80 CONTINUE 

LEND-K+LSTART-l 




0030 

NEXTL*. FALSE. 

c 



*• 

0031 

0032 

9999 RETURN 

gNfi 
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ORIGINAL PAGS 13 
OF POOR QUALITY 


I Pen TRAN I V-PLUS V02-»tE 14 S i M 30- OCT-BO JP^EJ 

MDRHI8 FTN /TR DLOCKS/wi 


0001 

5 

§UB*«buTINi MDAhi 4(LUN2. FLNtHM. UI^ME) ■ . 


c 

THIS ROUTINE WILL WRITE THE HISTOCRAHS FOR "EACH lAND C4>7] 
TO _THC I MDACS . HEADER . F ILE _ 

0002 

c 

IMRLICIT INTEOER <A-Z) 

0003 
_ 0004 _ 


REAL H4.H9. HA. H7 
LOOICAL*! PLNIMMU) 

0009 

0004 


DInENStON UIFHR(193A> 

COMHON/H IST/H4 ( 296) . H9< 29A ) . HA< 296) . H7 ( 294) 


c 

c 

OFEN THE IMPACS HEADER FILE 


• 

OOPy 0REN(WNIT» LUN 2. NAME-PX N ir« H, TYPE- 'OLD FORM- 'UNFOWWATTED ^ 

♦ ACCESS-'DtHECT'i RECDRDSl ZE»y6a > 

S 



c 

WRITE THE HISTOGRAnS FOR BANDS C4-7] 

0008 


WRITE(LUN2'3)H4.H9.H4 


M— M 



c 

C THE WEWAININQ BAND HISTOGRAMS ARE NOT USED IN THIS APPLICATION 


0010 DO 10 I»5. to 


0011 

0012 


WR1TE(I.UN2'I)UIFHR 
10 CONTINUE 


c 

CLOSE THE U10ACS HEADER FILE 

0013 

c 

CL0SE<UNIT“LUN2) ^ 

0014 

c 

RETURN 

0019 


END 
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ORIGINAL PAG2 !S 
OF POOR QUALITY 


1 PQRTRAN 

IV-PLU8 002-51E 14.2252 30-0CT-80 PACE 1 


HTlrTTN 

/TP. OLOCKI/WR 



C 

SUBROUTINE HIST (BUF*. WhAn, TBT) " ” 



T~ 

tHis Routine wIll ^uilo histoorahs i^or each Channel, usino 

THE VALUE FROM THE FULL FRAME DATA 


oooa 


IMPLICIT INTEGER <A-2> 


0003 




REAL H4.H5.H6.H7 
LOGICAL*! DUF( 1 >, MI8N0(2) 


o o 
o o 
oo 


COMKON/H I ST/H4 ( 256 > > H5 < 256 ) . H6 ( 256 ) < H7 ( 256 ) 
C0MN0N/LABEL/NAME(3). MISNO. E 



c 

ROUTINE ID 


. 0007 

c 

NAME < 1 )«'HI ' 


oooe 

0009 


NAME<2»-'8T' 
NAME<3)-' ' 



c 

JL. 

BUILDING HISTOGRAM FOR EACH CHANNEL 


0010 

c 

GOTO( 10. 20. 30. 40)RCHAN 


0011 

c 

10 CONTINUE . 


OOlS 

0013 


K-0 

DO 15 M-STARTP. STOPP 


0014 

0015 


K-K*l 

N-STARTP-TST *-K 


0016 

0017 


PIXEL-BUF(N). AND. "377 
H4<PIXEL+l )-H4(PIXEL+l)*l 


0018 

0019 


15 CONTINUE 
GOTO 9999 


0020 

c 

20 CONTINUE 


0021 

0022 


k«6 

DO 25 M-STARTP. STOPP 


0023 

0024 


K-K+1 

N-STARTP-TST+K 


0025 

0026 


PIXEL«flUF(N). AND. "377 
H5<PIXEL>1 )-H5(PIXEL-fl )♦! 


0027 

0028 


25 CONTINUE • 
GOTO 9999 


0029 

c 

30 CONTINUE 


0030 

0031 


K»0 

DO 35 M-STARTP. STOPP 


0032 

0033 


K>K*1 

N»STARTP-TSTfK 


0034 

0035 


PIXEL«BUF<N). AND. "377 
H6<PIXEL*-n«H6<PIXEL-^l )♦! 


0036 

0037 


35 CONTINUE 
GOTO 9999 


0038 

c 

40 CONTINUE 


0039 

0040 


K»0 - ■ 

DO 45 M-STARTP. STOPP 


0041 

0042 


K-K+1 

N-STARTP-TST+K 
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fortran IV-PLUS V02-51E 14:22:52 3O-OCT-e0 PAOC 3 

MIST FTN /TR: BLOCKS/UR 


ohiqinal paoc is 

OF POOR QUALITY 
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TT17 1350. 403A B/-SP«=SVO: C320. 401HIST 



ORIGINAL PAG?- 53 
OF POOR QUALITY 


i FORTRAN IV-PLU5 V05-51E M: 29- 12 30-0CT-B0 PAGE 1 


PARTP 

FTN 

/TR: BL0CK5/WR 

”o6o: 


Subroutine PARTPcptcH. STAftt^, sTOpp.stbl^. stripc. tse. tst. pixels. 

♦ 1. PATCH. PSIZE) 


C 

THIS ROUTINE WILL ADJUST THE PARAMETERS IF A PATCH IS DIVKikO BETWEEN 


C 

C 

WO ccT Strips. 


C 

l6£al vaSiablSS; 


c 

c 

CCtNO- tHE CCT StRIP NUMBSr'LaSY S'TORSO 

piece- an array CONTAINING THE PARAMETER VALUES TO USE FOR THE 


c 

c 

REMAINDER OF THE RIGHT MOST PATCH ON NEXT CCT STRIP 

0002 

0002 


IMPLICIT INTEGER (A-Z) 
LOGICAL PTCH 

0004 

0005 


DIMENSION PIECE<2. 4) 

DATA CCTNO/O/. F.ECE/8»0/ 

0006 

c 

IFISTRIPC. EG. 1. AND STOP GT. TSE)00T0 20 

0007 

0008 


IF(STRIPC, EG, 1 )GOTO 9999 ' 

IF(. NOT. PTCHIGOTO 10 

0009 

c 

1F<PIECE(STRIPC-1. 2). NE, PATCH-1 )00T0 10 


c 

c 

THE CURRENT PATCH IS THE PORTION OF THE UNPROCESSED RIGHT MOST PATCH 
FROM PREVIOUS CCT STRIP 

0010 

c 

PTCH-. FALSE. 

0011 

0012 


STARTP-TST 

ST0PP«TST-*-(PSI2E-PIECE<STRIPC-l, 1 ) )-l 

0013 

0014 


NPIXEL-(PIECECSTRIPC-1. 2) -1 ) *PS I ZEX 70+512* ( I-l ) ) 
PIXELS-NPIXEL-^PIECE<STRIPC-1, 1) 

0015 

c 

PATCH-PATCH-1 

0016 

0017 


10 CONTINUE 

IF < STOP. LE. TSE) GOTO 9999 

0018 

c 

20 CONTINUE 

0019 

c 

1F<CC1N0. EQ. STRIPOGQTO 9999 


c 

c 

SAVE THE NUMBER OF -PIXELS PROCESSED FOR THE RIGHT MOST PATCH ON 
CURRENT CCT STRIP 

0020 

c 

PIECE<STRIPC. l)-<3T0PP-STARTP)-«-: 

0021 

0022 


PIECECSTRIPC. 2)=PATCH 
CCTNO-STRIPC 

0023 

c 

9999 RETURN 


0024 END 
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PARTP. FTN /TR: BLOCKS/MR 
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FORTRAN IV-PLUS VQ2-51E 


14. gS. 15 


30-0CT-a0 


PACE 3 



ffiv 


1 

1 FORTRAN 

IV- 

PLUS V02-91E 14 26. 10 30-0CT-Q0 

PAGE 1 

INTOCT. FTN 

/TR: BLOCKS/UR 


0001 


SUBROUT I NE I NTOC T ( 1 NTEG. OC TAL ) 



c 

JL— 

THIS PROORAM WILL CONVERT A DECIMAL NUMBER 
IT CAN HANDLE VALUES IN DECIMAL FROM 1 TO 

TO OCTAL 
4099 


c 

DECIMAL NUMBER LESS THAN 4096 



c 

INTEGER OCTAL 

* 

0003 

0004 


IFdNTEG LT. 512)GO TO 10 
K-MOD< INTEG. 912> 


0009 

0006 


L-MOD<K. 64) 
M»MOD<L, 9) 


0007 

oooa 


I0CT*INTEC/912*1000+K/64*100+L/8*10-^M 
CO TO 100 



c 

c 

DECIMAL NUMBER LESS THAN 912 GE 64 


0009 

c 

10 

IF< INTEG. LT. 64 >G0 TO 20 


0010 

0011 


K>MOD( INTEG. 64) 
L«MOD<K. 8) 


0012 

0013 


IOCT»1NTEG/64*1OO+K/0*1O+L 
GO TO 100 



c 

c 

DECIMAL NUMBER LT 64 GE 8 


0014 

c 

20 

IF< INTEG. LT. 8)G0 TO 30 


0019 

0016 


K-M0D< INTEG. 8) 
I0CT»INTEG/8*10-*-K 


0017 

c 

GO TO 100 



c 

c 

DECIMAL NUMBER LT 9 


0018 

0019 

30 

100 

lOCT-INTEG 

OCTAL-IOCT 


0020 

0021 


RETURN 

END 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


I FORTRAN IV-PLUS V02-91E 14:26:27 30-0CT-B0 I^E I 


UNVRSL. FTN 

/TR; BLOCKS/WR 

0001 

C 

SUBROUTINE UNVRSL (LUN2. SUNEL> FLAT. FLONO) ^ 


C 

<? 

THIS ROUTINE WILL BUILD AN SPU-UNI VERSAL HEADER FOR IMAGERY FILES. 

0002 

0003 


IMPLICIT INTEGER <A-Z) 

REAL FLAT.FLONO. ILAT. ILONG, RLAT. RLONO 

0004 


LOGICAL*! UHEADR(3072). TPFORM( 32 > . BUORO< 2) > M1SN0(2> . 
* ITIME(8>. ASCI JD<3>. ASCIIV(2>. 0(2). 

0009 


* BWLAT(2). BWL0NG<2) 

E(JUI VALENCE (P I XLN, DWORD) , (SUNELA. B>, < BWLAT, DLAT ) , 

0006 


* (BWLONG. DLONO) 

CaMM0N/LADEL/NAME(3>. MISNO. E 

0007 

0008 


COMMON/LATLON/LATD. LONOD 

DATA UHEADR/3072*0/, TPF0RM/32»0/ 

0009 


DATA TPFORM/'S'. 'P'. 'U', ' 'U'. 'N'. 'I'. 'V', 'E'. 

♦ 'R'. 'S'. 'A'. 'L'. ' '. 'F'. 'O'. 'R'. 'M'. 'A'. 


c 

♦ 'T', ' ' '. ' '. ' ', ' ', ' '. ' '/ 


c 

c 

FORMAT LISTINGS 

0010 

0011 


1 FORMAT! 13) 

2 FORMAT! 12) 


c 

c 

FIRST LINE OF UNIVERSAL HEADER 

0012 

c 

DO 9 K«l. 32 

0013 

0014 


UHEADR ( K ) -TPFORM ! K ) 
9 CONTINUE 


c 

c 

THESE SYSTEM ROUTINES ACQUIRE THE PRESENT DATE AND TIME 

0019 

c 

CALL TIMEdTIME) 

0016 

c 

CALL IDATE! IM. ID. lY) 


c 

THIS ROUTINE WILL CONVERT GREGORIAN CALENDAR DATE TO JULIAN 

0017 

c 

JDAY-JULIAN! IM. ID. lY) 


c 

c 

CHANGE JULIAN DATE. AND YEAR TO ASCII FORM 

0018 

0019 


ENCODE (3.1. A5C I JO ) JDAY 
UHEADR ( 39 ) -ASC I JD ! 1 ) 

0020 

0021 


UHEADR ( 36 > =ASC I JD ! 2 ) 
UHEADR ! 37 ) «ASC I JD ! 3 ) 

0022 

0023 


ENCODE (2.2. ASC 1 1 Y ) I Y 
UHEADR!33)-ASCIIY!1) 

0024 

c 

UHEADR ! 34 ) «A5C 1 1 Y ( 2 ) 


c 

c 

PLACE TIME (HOURS. MINUTES) INTO ARRAY 

0029, 

0026 


J»0 

DO 10 1-38.42 

0027 

0028 


\ 

UHEADR( I )«ITIME( J) 

0029 

c 

10 CONTINUE 



ORIGINAL PAGE tS 
OF POOR OUAL’.TY 


i FORTRAN 

IV- 

-PLUS V02-51E 14:26.27 30-0CT-B0 PAGE 2 

UNVftSL. FTN 

/TR OLOCKS/Wft 


c 

PLACfe THE MlSstClN NUtlfi^A INTO HEADER 

0^33 

c 

UHEA0R(66)>^ 


"5“ 

c 

^LAClb ACTIVE CHANNEL FLAOSi 1 ■ ACTIVE, 6«N0T ACtlVC 

oo5I 

0032 


UHEA0R<81 )«240 
UHEADR(90>-4 


c 

c 

PLACE PIXEL SIZE 

0033 

c 

UHEADR(91 )«a 


c 

c 

PLACE VIDEO DATA STARTi INDICATING DATA FOLLOUS HEADER 

0034 

c 

UHEADR(93)«1 


c 

c 

PLACE PIXELS PER LINE PER BAND 

0035 

c 

UHEADR(96)-2 

0036 

c 

UHEADR<97)«0 


c 

c 

PHYSICAL RECORD SI ZE< MULTIPLE OF 180 
ACTUAL SIZE IS ( 16«32»4«2048 BYTES) 

0037 

c 

PIXLN«2560 

0038 

0039 


UHEADR (100) -DWORD < 2 ) 
UHEADR<101)-DUORD(1) 


c 

c 

PLACE THE DATA SET SIZE; NUMBER OF RECORDS PER LINE 

0040 

c 

UHEADR ( 104) -1 


c 

c 

ANCILLARY BLOCK SIZE 

0041 

c 

UHEADR < 106 )-70 


c 

c 

BINARY START PIXEL 

0042 

c 

UHEADR < 109) -1 


c 

c 

BINARY STOP PIXEL 

0043 

c 

UHEADR ( 110)-2 

0044 

c 

UHEADR < 1 1 1 )=-0 


c 

c 

GENERATION COMPUTER WORD SIZE 

0045 

c 

UHEADR <753) -16 


c 

c 

DATA SETS PER PHYSICAL RECORD 

0046 

c 

UHEADR < 1778) -1 


c 

c 

BANDS IN PRIMARY RECORD 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


FORTRAN IV-PLUS VQ2-51E 14: 27 3O~0CT-flO PAGE 3 


UNVRSL 

Tfsr 

/trTlockITwR 

5047 


UHEADR< 1786>«4 - 

' 

c 

.. . 


c 

— C - 

BYTES PER LINE PER BAND 

' 0048 


UHEA0R(17e7)-8 

0049 

- 

UHEA0R< i7ea>«o 


c 



--C 

PIXEL SKIP FACTOR 


c 


0050 

- 

UHEADRl 1790)-1 


c 



c 

LINE SKIP FACTOR 


c 



— 

UHeADR< 1792)»1 


c 



-C 

MISSION NUMBER 


c 


. 0058 

_ 

UHEAOR( 1817)*E 


c 



c 



c 


0053 

- 

UHEADR( iai9)>2 

0054 


UHEADR(1820)>0 



-c 



c 

PIXELS PER LINE PER BAND 


--C.. 

- .. .. 

0055 


UHEADR< 1821)-2 

0056 - 

- — 

UHCADR(1822)-0 


c 





LINE SKIP FACTOR 


c 


0037 


UHEADR(iaa3)»l 


c 


.... 

- c 

PIXEL SKIP FACTOR 


c 


0058 



UHEA0R(1824)-l _ . 


c 


- 

— c 

ACTIVE BAND FLAGS 


c 


0059 • 

. - . 

UHEADR<l825)-''360 


c 



c 

CLOUD COVER 


c 


0060 

- 

UHEADR< 1827)>1 


c 



c 

INGEST BAND STATUS 


c 


0061 - 


UHEADR( 1828)-2 

0062 


UHEAUR( 1829)»2 

0063 


UHEADRC 1S30)«2 

0064 


UHEA0R(1831)>2 

- 

c 



c 

SUN ELEVATION ANGLE 


c 


0065 


SUNELA*SUNEL 

0066 


UHEADR( 1033)*B(2) 
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I FORTRAN !V~PLUS V 0 5-SIE 14:24:27 30-0CT-B0 PACE 4 


UNVRSL. ^TN 

/TR: DLOCKS/WR - - 

0067 


UHEADR( 18341-0(1)- - ■ 


c 

c 

LATITUDEi CENTER OF FULL PRAME 

TOSS 

0069 


ILAT-FLAT ’ 

DLAT-FLAT 



i^LAt-< < lL'At-DL“Af")*i60)/"60 

0071 


UHEAOR< t989)-t.AT0 

o575 

0073 


UHEADR< 1986)-PULAT(4) 
UHF.ADR < 1 987 ) -BVJLAT ( 1 ) 

0074 

c 

UHEADR< 19B9)-RLAT -- 

- 

c 

c 

LONGlTUDEi CENTER OF FULL FRAME — 

0079 

0076 


lUONO-FLONO 

DLONC-FLONG 

03T7 

007B 


RLONG-(<ILONG*OLONG>»100)/6a ' 
UHEADR ( 1 990 ) -LONGD 

0079 

0080 


UHEADR(1991 )-BWL0N0<2> — - ” 

UHEADR ( 1 992 ) -BULCNG ( 1 ) 

0081 

c 

UHEADR <1994) -RLONG 


c 

c 

WRITE THE BUFFER TO THE HEADER FILE 

6664 

“c 

. ,.WRITE<LUN2'2)UHEADR - — - 

56S3 

0084 ■ 

— 

RETURN - -- - - - 

END 

— 
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F0RTRAN_IV_-PLUS V02-51E 14 24:27 30-0CT-B0 PAGE 4 

UfJvRSLTFTN /TR: 0L0CK8/UII 


f6TAi. 8PAC^ AuuocAT^'fi - iA7A 

. tT17: C320. 403A. 8/-8P-8Y0: t'2i6. 462'JNVRSl 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


I FORTR AN IV-PLU3 VOS-glE 14 gg 47 30»0CT"a0 RAQI I 

FRONT /TR BLOCKi/VM 




T< 

C 

Z 

c 

T 

c 

T 

c 

T 

c 

T 

c_ 

c 

c 

T 

c 

T 

jl 

c 

c 

c 

c 

c 

c 

T 

c 

T 

c 

T 

_c 

c 

_c. 

c 

_c 

c 

_c 

c 

_c 

c 

c 

_c 

c 

_c 

c 

c 

_c 

c 

_c 

c 

c 


LOCKHEED ELECTRONICS COHRANY. INC. 

-ifiMflTe ssns m A^ikLiyATi dNi "Ofr " 

< R 8 A L ) 


NAME; LEFT JUSTIFY AStll ^HARAdTCft STRING (►AbNY) 

PRtiCUf^^D FR(HM; NASA. j6hNSON ^AaC£ '^^nTEM 

DLPQ. 17. DATA TECHNIQUES LAB 

ON: 01-MA^-^S 


REVISED BY: 

ooooo< 

B. L. TAYLOR 


ON: PURPOSE FOR REVISION: 

o< <><.-;<>oooooooo<!K>o<^<i><!J><><ii><i 

a2-MAY-7B DR I NO TO RSAL STANDARD 


PR OGRAM DESCRIPTION: 

oooooooooooooooooooooooooooooooc 

LEFT JUSTIFIES AN ASCII CHAfMCTER StRINO REMOVES ANY LEADING 
BLANKS AND BLANK FILLS THE ARRAY AFTER THE LAST CHARACTER 18 
SHIFTED. 

ENTRY POINTS: 

00 <> 00000000000000 <><><> 000 <><> 0000 <><>< 


FRONT 


PROGRAM INITIATION OR CALL SEQUENCE: 

0'::>oo<>ooooo<>oooo«:'.;000o<!>o<l><>«d>oo<!><i5^<di>^^^^ 


CALL FRONT < I ARRAY. NUM) 

CALL ARGUMENTS; ^ FORMATS; 

OOOOOOO 0<>00000<.'^000000<>00<>0< 

• 

IARRAY ^ BYTE STRING (ARRAY) TO BE LEFT JUSTIFIED. 

NUM MAXIMUM NUMBER OF BYTES IN THE STRING (ARRAY). 


INPUTS: FORMATS. 

OOOOOOO ooooooooc^ooooooooooc 


NONE 


OUTPUTS: 

OOOOOOO 


FORMATS 

OOOOOOOK 


>000000 ;>oooo< 


NONE 


SPECIAL REMARKS: 

ooooooooooooooooooooooooooooooo< 


NONE 
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Oh.C.rMI, fUCL ,'1 
OF PuOi^ QU#iLl f y 


i FO>tTRAN W->»m8 V08-9tg 14 gg. 47 30>0CT-e0 PAQM 8 


?Sof5T. FTN 

/TM; 8LbCK8/UR 


X 

MOUTINES CALLED: 


c 

X 

<>0<><><><>0<><>00<3 K><><><>000<>0<t><>0<iK><lK!>^>0<i 


c 

X 

NONE 

C 



c 


X 

LEFT JUSTIFY A CHARACTER STRING 

0001 

0002 


SUBROUTINE FRONT! I »N) 
- IVTC 1 

0003 

X 

DIMENSION I(N) 

0004 

0009 


00 1 J-l.N 

IF(KJ) NE. “40)00 TO 2 

0006 

0007 

1 

CONTINUE 

RETURN 

OOOB 

c 

X 

K-1 - 

0009 

0010 


00 3 J-J. N 
I(K>-I(J> 

0011 

001? 

3 

K«K^1 

IF(K. OT. N) RETURN 

0013 

C 

DO 4 J>K. N 

0014 

0019 

4 

I ( J>-“40 
RETURN 

0016 

C 

END 
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ORiC;n/'.’. PAGE '3 
OF POCR QUALITY 


i FORTRAN IV-PLUS 

uoa-siE 

14:28:33 30-0CT-80 PAGE 1 

JULIAN, FTN 

/TR: 8L0CKS/WI* 



S 

C 

c 

LOCK^«ED ELECTRONICS COMPANY. INC. 

C 

REMOTE SENSING APPLICATION^ LAB. 

c 


( 

R a A L ) 

c 

c 




c 




c 

DEVELOPED DY: 

TED KELL 


c 

c 

ON: 



c 




c 

REVISED OY; 

ON: 

PURPOSE FOR REVISION: 

c 

<><><><><>< 

o< 

o<x>o<? <x>o«:x>oo<.-> ooocxxx 

c 

S OBRIEN 

6/78 

ADD COMMENTS FOR STANDARDIZATION 


PROGRAM DESCRIPTION: 

oooooooooooooooooooooooooooooooc 


c 

c 


THIS FUNCTION COMPUTES JULIAN DAV INTO INTEGER FORMAT 


C 

C 


ENTRY POINTS: 

ooooooooooooooooooooooooOcSo-dxxxi^c 

NONE 


c 

c 


PROGRAM INITIATION OR CALL SEQUENCE.: 


C 

C 


000000000<><.>00000<>0000000<X>00<XXX 
WAR-w»JLIAN(MQN. IDAY. lYR) 



c 

CALL ARGUMENTS: 

FORMATS: 


c 

ooooooo 

o<>«:>o«;xxx.-5 «:xx>ooooo^:xxx 


c 

MON 

MONTH NUMDER IN INTEGER 


c 

IDAY 

DAY NUMDER IN INTEGER' 


c 

lYR 

YEAR NUMBER IN INTEGER 


c 

c 

INPUTS: 

FORMATS' 

2 

c 

oooo<xx> 

oo«:>oooo<>ooooooooooo< 

1 

c 

MON 


0 

c 

IDAY 


> 

c 

lYR 



c 



7 

c 


' 

« 

c 

OUTPUTS; 

FORMATS: 

% 

c 

<-..>000 

«:x>ooooo<>oo<.x'.><:;><.‘>o<.x:'<>o< 
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FORTRAN TV-PLUS VOg-5iE l» gB 33 30-0CT-B0 PACE g 

JULIAN. FTN /TR: BLOCKS. tJR 



0001 INTEGER FUNCTION JULIAN<N. D, V) 








FORTR AN IV-PLUS V0 2-5tE 14:28:33 30-OCT-BO PAPE 3 

JULIAN FTN ~ /TR: DLDCKS/UR 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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ORlGif^AL PAGE SS 

OF POOR QUALITY 


i FORTRAN IV-PLUS V0 2-S1E 14 58:1 A 30-0CT-B0 P4G E I 

HPROS, FTN /TR BLQCKS/WR 


C+1 

c 

c 

c 


LOCKHEED ELECTRONICS COHPANY. INC. 


REMOTE SENSING APPLICATIONS LAB. 
< W S A L ) 


C 

C 


NAME; SUBROUTINE HPROS 


DEVELOPED BY: TED KELL 
ON: 


REVISED BY: 

ooooo< 


ON; 

o< 


PURPOSE FOR REVISION: 

OOP :>ooooooo<>oo<><><><><>< 


SUSANNE OBRIEN 6/7S 


STANDARDIZATION 


C PROGRAM DESCRIPTION; 

_c oo«poooooooooooooooooooooooooo<><>< 


THIS ROUTINE WILL BUILD A COMMON BLOCK OF TAPE HEADER DATA FOR 
SUBSEQUENT PROCRAM USE 


ENTRY POINTS: 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOCX 


PROGRAM INITIATION OR CALL SEQUENCE; 

ooooooo<>ooooo«::!::>oooooooooo<>oo<^ oo< 


CALL HPROS < LUN, BUF, FMT, EOF, PR TV ) 


CALL ARGUMENTS: 

ooooooo 


FOR.IATS: 

ooooooooooooooocxxxx 


LUN 

BUF 


n LOGICAL UNIT NUMBER 

BYTE STRING BUFFER CONTAINING HEADER DATA 


FMT 


II FORMAT OF HEADER 
1*UN1 VERBAL 


2>LANDSAT 

3»LARSYS 


EOF 

PRTY 


n END OF FILE INDICATOR 
n *'^RITY COUNT 


INPUTS: 

OOOOOOO 


FORMATS; 


OUTPUTS; 

OOOOOOO 


FORMATS: 


HCOM 


COMMON AREA OF HEADER DATA BUILT BY HPROS 


SPECIAL REMARKS; 

OOOOOOOOOOOOOOO 


IXXXXXXXXXXXXX 


C ROUTINES CALLED; 

c oooooooo<xx>oooo<xxs ooo ,xx: ;:>oooooo< 


IBYTE 

MTREAD 


CONVERTS DATA FROM ASCII TO INTEGER 
MAG TAPE READ ROUTINE 
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ORiGINAL PACE *3 

OF POOR QUALITY 


I FORTRAN IV-PLUS V05-51E 14:28:16 30-QCT-B0 PAGE 2 


HPROS, FTN 

/TR: BUOCKS/WR 


C 


0001 

0002 


SUBROUTINE HPROSILUNI, HEADER. FMT. EOF. PRTY, IPRM) 
IMPLlCITINTE0ER<A-2) 


5o03 

0004 


INTEGER HC6nB<26) ! 

DIMENSION ISTAT<2>. IPRM(6) 


0005 

0006 


BYTE HEADER( 1 ). SCNID<12>. MISNC(2). ISTATD<2) 
COMMON/HCOM/SS. SE. LS, LE. NRPDS. NDSPR, NCPR. NRPC. ANCL. NC. NS. 


0007 


1 NBIT, DOI, NCAR. SVD, RSIZ, PSKIP, HSI2. CALP. CERR. NSPR 
C0MM0N/LABEL/NAME(3). MISNO. E 


0008 

0009 


EQUIVALENCEISS. KCOMBd > ). < ISTAT< 1 ). ISTATB) 
DATA RLB/"1000/. IEE0F/-10/ 


0010 

0011 


B DUF ( I ) - 1 B YTE < 0. HEADER < I ) ) *256+ 1 BYTE < 0. HEADER < I +1 ) > 
NAME< 1 )*'HP' 


0012 

0013 


NAME(2)»'R0' 
NAME<3)«'S ' 


0014 

0015 


IF (CERR. NE. 0)RETURN 

IF<FMT. EQ. 1) CO TO 1 (FORMAT-UNIVERSAL 


0016 

0017 


IK(FMT. EQ. 2) GO TO 2 ! FORMAT-LANDSAT 

IF<FmT.EQ. 3) CO TO 3 ! FORMAT-LARUYS 


OOil 

0019 


CERR - 9 
RETURN 


0020 

C 

1 

BUILD HCOM FbR UNIVERSAL' FORMAT TAPE 
SS - BBUFdOS) 


0021 

0022 


SE-BBUFdld) 
LE - 0 


!S623 

0024 


NRPDS-HEADEHd04) 
NDSPR-HEADER d 778 ) 


SSSs 

0026 


NCiPR-HfADEfr < ioT) 
NRPC-HEADERd03) 


0027 

0028 


ANCL»BBUF( 105) 
NC-HEADER(90) 


0029 

0030 


NS-BBUFd787) 
NDIT»HEADER(91 ) 


oo5I 

0032 


DOI ■ HEADER (107) 
NCAR»BBUFd785) 


0033 

0034 


SVD - BBUF('?2l . 

IF (SVD. EQ, 0)SVD-1 


0035 

0036 


RSIZ » BBUFdOO) 
PSKIP « 0 


0037 

0038 


HSIZ » 1530 
CALP - 0 


0039 

0040 


NPSL-B21 (HEADER. 96. 97) 
GOTO 9999 



c 

c 

BUILD HCOM FOR LARSYS FORMAT TAPE 


0041 

c 

3 

SS - 1 



• 

SE - BBUF(23) 
LE - DBUF(39) 


0044 

0045 


NRPDS -1 
NDSPR - 1 



« 0046 NCPR » BBUF<19) 

» 0047 NRPC * 1 
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OF POOR QUAUTY 


FORTRAN 

IV 

-PLUS V02-51E 14; 29; 16 30-OCT-80 


PACE 3 

HPROS. FTN 

/TR. DLOCKS/UR 



0046 

0049 


ANCL - 4 
NC > NCPR 



0050 

0051 


NS - SE 
NOIT > 8 



0052 

0053 


001 • 0 
NCAR - NCPR 



0054 

0055 


SVD s 1 

RSIZ » NC*NS«4 



0056 

0057 


NS-N3-6 

SE«SE'i. 



0056 

0059 


PSKIP - 0 
HSIZ • 400 



0060 

0061 


CALP ■ 6 
GOTO 9999 




c 

c 

BUILD HCOM FOR LANDSAT FORMAT TAPE 



0062 

c 

2 

NS - DBUF<39>/4 



0063 

0064 


I>HEADER{ 14) 

SS > NS»< IANO( I. 15>-1 >-«-l 



0065 

0066 


SE - SS-^NS-l 
LS - 1 



0067 

0066 


LE - 2340 
NRPDS • 1 



0069 

0070 


NDSPR « 1 
NCPR » 4 



0071 

0072 


NRPC - 1 
ANCL - 0 



0073 

0074 


NC - 4 
MBIT - 8 



0075 

0076 


DOI > 2 
NCAR > 4 



0077 

0076 


SVD =» 1 

RSIZ - BBUFU7) 



0079 

OOSO 


PSKIP =» 9 
HSIZ » 338 



0081 

c 

CALP » 0 



0082 

c 

4 CONTINUE 




c 

c 

READ ANCILLARY RECORD S. CONCATENATE TO 
TO FORM HEADER RECORD FOR LANDGAT TAPE 

ID RECORD 


c 

c 

PUTTING STARTING ADDRESS OF HEADER ARRAY INTO IPRM 

ARRAY 

C 

C 


c 

c 

READING FROM TAPE 



0033 

0034 


IPRtK 1 ) = IPRM( n+40 
IPRM(2)»624 



0085 

c 

CALL WTQIOCRLB, LUNl, 1, . ISTAT. IPRM. ISW) 




c 

c 

ERROR CHECK 



0036 

0037 


5 IF< ISTATDU ). GE. *‘000)G6 TC 999C 'READ 

IFUSTATOd), EO. IEE0F)E0F»1 'END OF 

OKAY a" 
FILE=" 

001 

366«-l0 


1 
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ORIGINAL PAGE [S 


I FORTRAN IV-PLUS V08-51E t4; 27. 54 30-OCT-fiO RAOC I 

HCiRlACnSTFJ /TR: 8LO&KS/WR 


z 


u 

t1 

10 

f 

0 

t 


C 


5 

c 

LOCKHEED 

ELECTRONICS COMPANY. INC. 

c 

c 

REhOTE sensing APPLlCAftONS LAfi. 

< R S A L > 

d 

c 


- 

c 

c 

NAME; H0REa6 


c 

c 

DEVELOPED OY; SUSANNE 
ON: 6/7B 

OBRIEN 

e 

£_ _ . 

c 

c 

REVISED OY: ON; 

ooooo< o< 

PURPCiE FOR REVISION: 

OOOOOOOOOOOOOOOOOOOC 

c 

c 

c 

c 

PROGRAM DESCRIPTION: 

OOOOOOOOOOOOOOOOOOOOOC 

c 

c 

THIS routine will READ HEADER DATA RECORDS FROM A FOREIGN 
IMAGERY TAPE IN ONE OF THE FOLLOWING FORMATS: UNIVERSAL. 

c 

c 

LANDSAT 1 OR 2. LARSYS 2 OR 3 
ENTRY POINTS: 

c 

c 

oooooooooooooooooooooooooooooooc 

NONE 

c 

c 

PROGRAM INITIATION OR CALL SEQUENCE; 

c 

c 

oooooooooo 

CALL HDREAD(LUN, IFRM 

0000<.50000000000000000< 
. HEADER, BFSZ, EOF, PC. IE, BUFFER) 

c 

c 

CALL ARGUMENTS; 

FORMATS; 

c 

c 

ooooooo 

LUN 

0000000000<?00000000< 
11 LOGICAL UNIT NUMBER 

c 

c 

IFRM 

01 FORMAT OF IMAGE TAPE 
1 -UNIVERSAL 

c 

c 


2- LANOSAT 

3- LAR5YS 

c 

c 

• 

O-ERROR HEADER RECORD NOT 
ACCEPTABLE SIZE 

c 

c 

HEADER 

BFSZ 

BYTE STRING DUFFER FOR HEADER DATA 
11 SIZE OF READ IN BUFFER 'BUFFER' 

c 

c 

EOF 

PC 

11 END OF FILE 1-END 
PARITY COUNT 

c 

c 

IE 

01 ERROR CODE 
IE-0 MEANS NO ERROR 

c 

c 

BUFFER 

INPUTS; 

BYTE STRING DUFFER TO READ DATA INTO 
FORMATS: 

c 

c 

OOOOOOO 

OUTPUTS: 

<:-'0 <;<>ooooooooooooooo< 

FORMATS: 

c 

c 

OOOOOOO 

HCOM 

000<>0000<50000000000< 

COMMON DATA AREA BUILT BY HPROS (CALLED) 

c 

c 

SPECIAL REMARKS. 


c 

OOOOOOOOOO 

oooo'w' ooo<;? oo oo<>oooooo< 


C 
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1 FORTRAN 

IV- 

■PLUS V02-91E 142754 30-OCT-SO PP 51 3 

HDRCAO PTN 


/TR DLOCKS/WR 

0034 



CALL HPWitLUNl. HEADER, I#RM. fiOP. F5, "IPRH^ ''SulLfi W5 (Sm 
IF<CERR NE. 0)00 TO 100 

oo55 




0034 

0037 



IFCNRPC. LE. 1)00 TO 40 
NSPR-NS/NRPC ICALC ••IX PER REC 

0038 
003^ 



IF(M00',NS. NRPC). EQ. 0)00 TO 40 ‘EVEN DIVlfttdN 66 tO 4 

N8PR-NSPR4>1 

0040 

0041 

40 


0 ■ RSIZ 

IFIRBIZ. OT. DFSZ) CERR-12 

0042 
_ 0043 



IFICERR. EQ 12)G0 TO 100 ‘REC TOO DIO FOR BUFFER 
IPRM2<2)-R8IZ 


C 

c 


CALL HDCODE FOR HEADER DECODE OF HUFFER 

0044 

c 


CALL HDCODE < IFRM. HEADER. IE, LVL) 


c 

PUTTING STARTING ADDRESS OF BUFFER INTO IPRM2 ARRAY 

0049 

c 


CALL GETADRi IPRM2, DUFFER) 


c 

c 

READING IMAGE DATA FROM TAPE 

0044 

c 


CALL UTOIO<RLO. LUNl. 1. , ISTAT, IPRM2, I8W) 


TT" 

c 

ERROR CHECK 

0047 

c 

49 

IF( ISTATD( 1 ) OE. '‘OOOCOTO 90 !READ OKAV-1 

0049 



IF< ISTATDi 1 ). EQ. lEEOFJEOF-1 !EN6 OF FILE? ^ 

IF (EOF. EQ. 1 ) RETURN 

oo3o 

0091 



i)^( isfAtD< 1 ). EQ. -4)P(i-pc+i ! Parity cn^k-^ ~ - -- — ■■■ 

IF( ISTATDC 1 >. EQ, -4)00 TO 90 

0093 

50 


GO t 6 ll6 ‘READ ERROR ^ ■ 

IF ( IFRM. EQ. 2)00 TO 60 ILANDSAT TAPE 

0094 

T" 


CALCULATE LINE START NUMBER 
LS-LSLOC< IFRM) 

5553 

c 


LS«tDYTE(0. DUFFER(LS) ) »296+IBYTE(6. STiF^F^Eft (L6+1 )7 ^ • — 

VERIFY DATA 

0056 

T” 

60 


CHECK DATA ORDER INDICAT(jR • . 

IFCDOI EQ. 0 0R.*< IFRM EQ, 2 AND. DOI. EQ. 2) ) GO TO 70 

009? 

0053 



CERR -19 
GO TO 100 

0059 

0060 

75" 


if^(NBlt. EQ. 8) CO t6' so ! CHECK 4) bltS PER BYTE - - - 

CERR-16 

0061 

0062 

00 


GO TO 100 

IF(NRPC. OE. 1) GO TO 90 !# RECS PER CHAN 

0063 

0064 



NRPC-1 

NCPR-0 !# CHANS PER REC 

0069 

c 


CHECK # RECS PER DATA SET AND RESET NCPR IF NEEDED 
IFCNRPDS. OT 1) NCPP«(NC-NCAR)/(NRPDS-1 ) 

0066 

0067 

90 


IFCNDSPR, CE DCD TO 100 •» DATA SETS PER REC 
NDSPR-1 fNDSPR CANNOT -0 IN CALCS 

0063 

0069 

100 

lE-CERR !SET ERROR CODE 

RETURN 

0070 

0071 


Ic=*ISTATD< 1 ) !SET READ ERROR CODE 

RETURN 
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i 


! 


i 


>Z 



\ \ 


1 






li 

UJ 


I 



I 


I 


<1 



C,! 

SI 


<r r»i 
n '0. 

- pj 


o‘ 

s= 


•f ♦ 

o •» 

c* r, 


' -I 
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ORIGirJAI. PACa 
OF POOk QUAUTY 


FORTRAN IV-PLU3 V02-91E 


"HtiSEAD. FTN 


TtOCScksTwr' 


14; 27: S4 


30-0CT-80 


PAGE 6 


15n 3 Afi&~"S0dA'auTlNE5~REf eAeNeeS"’ 
"5ETTOB RBraBE RFTOS IfiVTfe riTAIfi 


TOTAL SPACE ALLOCATED ■ 001724 490 

NO FPP INSTRUCTIONS GENERATED 

. TT17: C320. 4Q3A. B/-SP«SYO: COHO. 4Q3HDREAD 
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FORTRAN IV-PLUS 

V02-51E 

14: 27. 

29 30>0CT-a0 PAGE I 

HDCODE. FTN 

/TR: BLOCKS/UR 



C 




c 




c 

LOCKHEED ELECTRONICS COMPANY. INC. 

c 

REMOTE SENSING APPLICATIONS LAB. 

c 



< R 8 A L ) 

c 

c 

NAME; HDCODE 

HEADER 

DECODE 

c 

c 

DEVELOPED BY: 

TED KELL 


c 

c 

ON: 



c 




c 

REVISED BY; 

ON; 

PURPOSE FOR REVISION: 

c 

ooooo< 

<>< 


c 

S. OBRIEN 

6/9/78 

ADD COMMENTS FOR STANDARDIZATION 

c 

S. OBRIEN 

8/3/ 78 

ADD PROCESSING FOR LANDSAT LEVEL 2 

c 



AND LEVEL 3 

c 

c 

c 

c 

PROGRAM DESCRIPTION: 



THIS ROUTINE MILL DECODE SELECTED PORTIONS OF TAPE 
"header RECORDS IN ONE OF THE *^OLLOWINO FORMATS 


UNIVERSAL 


LAr4DSAT I OR 2 
LARSYS 1 OR 2 


C ENTRY POINTS^ 

c oooooooooooooooooooooooo ooooooo< 


NONE 


r 

c 




c 




c 

PROGRAM INITIATION OR CALL SEQUENCE: 


c 

ooooooo< 

X. ‘.'OOOOOC n:; OOOOOOOO 0<><>0000< 


w 

c 

CALL HDCODE (IF.. 

H, IE, LVL) 

c 

c 


c 

CALL ARGUMENTS; 

FORMATS; 


c 

ooooooo 

oooo<>o<>ooooo<:>oooooo< 


c 

IF 

FORMAT IF HEADER • ?“> 


c 


l=UN I VERBAL 


c 


2=LAND3AT 


c 


3=LARSYS 

a 

c 

H 

HEADER RECORD EXrtCTLY AS READ FROM TAPE 

1 

c 

IE 

ERROR CODE 

0 

c 


0=N0 ERROR 

( 

c 


1 = IF NOT 1. 2. 3 

• 

c 

LVL 

II LANDSAT LEVEL INDICATOR 

7 

c 


1*LEVEL 1 

S 

c 


2-LEVEL 2 


c 


3-=LEVEL 3 
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1 FORTRAN 

IV 

-PLUS V02-51E 14:27:25 30-0CT-80 PAGE 2 

HDCOOE. FTN 

/TR; DLOCKS/WR 


c 

InI^uTS: FORflATSi . ' - .:V.. 

ooooooo oooo<>ooooooooo<o«:>o<><>< 


T~ 

c 

, - -- 


nr" 

c 

Outputs • PO^HaTs: ^ ^ 

ooooooo <>oooooooooooooooooo< 


T~ 

c 

CZSmH AREA FILLEb WITH DATA 

EXTRACTED FROM TAPE FOR USER 


IT 

c 

SPECIAL REMARKS: 


"T" 

c 



TT 

c 

_ — -• - .. — — — ' 


TT 

C 

I^OUTINES CALLES; 

000000000000<>0<;>0<>00000000000000< 


T" 

c 

^TAB. A2I. S2I. JOlIAN. ' LATLON 


IT 

c 

- • - 

C 

c 


c 

c 

‘ 


c 

. ... - — • ■ — — — - 


c 

C— 


COOl 

0002 


3UBP0UTINE HDCODEUF. H, lEi LVL) - 

REAL LATLON 

0003 

0004 


INTEGER DY, E. LArD#J-ONOD- - - 

BYTE HU) 

0005 

0006 


BYTE MONO, 12) .. - 

BYTE HB 

0007 

0008 


NTECER HK32) . — - 

REAL HR (16) 

0009 

0010 


INTEGER A21.B2I - 

INTEGER LVL 

0011 

0012 


INTEGER HOLDl - • 

INTEGER«4 H0LD2 

0013 

0014 


BYTE HLDB(4) _ - 

BYTE HLDA< 12), EXDAY<2) , SCNIDU2) , MISN0(2) 

0015 

0016 


C0MM0N/C0MH/Hfl(64) - 

COMMON/LATLON/LATD, LONGD 

0017 

c 

C0MMQN/LABEL/NAME(3), MISNO, E - -- 

HB IS ARRAY NAME FOR HEADER DATA ADDRESSED IN BYTES 


c 

c 

HR IS ARRAY NAME FOR HEADEr. DATA ADDRESSED IN REAL FORM 
HI IS ARRAY NAME FOR HEADER DATA ADDRESSED IN INTEGER FORM 

0018 

c 

HB, HR, AND HI ALL ADDRESS DATA IN ARRAY H 

EQUIVALENCE(Hn( 1 ), HI U ), HR< 1 ) ) 

0019 

0020 


EQUIVAL£NCE(HLDA< 1 ), H3U3) ) -- ■ 

EGUIVALENCE(HLDD< 1 ). H0LD2) 

• 

c 

c 

MON • ARRAY OF 3 LETTER EBCDIC CHARACTER ABBREVIATIONS 
FOR MONTHS OF THE YEAR NAMES 

0021 


DATA MON/ "321, "301, "225, "306, "205, "302, "',24, "201, "331, 
1 "301, "327, "331, "324, "301, "350, "321, "344, "325. 
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1 

i FORTRAN 

IV- 

PLUS V02-51E 14:27: 

25 30-0CT-a0 PAGE 3 

HDCODE. FIN 

/TR: BLOCKS/UR 




2 "321. "344, "323. "301. 

3 "326. "303. "343. "325. 

"544. "30^, "342. "305. "3S7. 
"326. "345. "304. "305. "303/ 

0022 

0023 


NAME(1)> 'HO' 
NAME<2>«'C0' 


0024 

0025 


NArlE<3>» 'DE' 
IE - 0 

*SET ERROR CODE NO ERROR 

0026 

c 

GO TO (1.2. 3). IF 

SET ERROR CODE TO 1 AND 1 ETURN TO MAIN PROGRAM 

0027 

0028 


IE - 1 
RETURN 



c 

c 

PROCESSING FOR LARSVS 

FORMAT TAPE 

0029 

c 

3 

HI<1>»1 ! ST (SCAN START) FIELD IN COMMON 

0030 

0031 


HI (2>-D2I(H. 23. 24) 
HI(2)*HI(2)-6 

! SE (SCAN END) FIELD IN COMMON 

0032 

0033 


HI(3)-1 ! LS (LINE START) 

HI(4)-2340 

0034 

0035 


HI <5)«Q2I (H. 19. 20) 
HI(6)-Q2I(H.23. 24) 

! NOCHAN (# CHANNELS ON TAPE) 

! NOSAMP (# SAMPLES PER SCAN ) 

0036 

C 

HI (6)-HI (6)-6 

BUILD SCNID (SCENE ID) 


0037 

0038 

31 

DO 31 I - 9. 12 
HLDB( I-8)«H( I ) 


0039 

0040 


HD ( 13) -"040 
DO 131 1-13.23 

- 

0041 

0042 

131 

HD( I-fl )-HB( I ) 
CONTINUE 


0043 

0044 

300 

ENCODE ( 12. 300. HLDA)HQL02 
FORMAT( 1 112) 

0045 

0046 


CALL FRONT (HLDA. 12) 

HR(7)-0 !SUNEL AND SUNA2 

0047 

C 

SUNEL (SUN ELEVATION) 
HI ( 15)-D2I (H. 63. 64) 

SUNAZ (SUN AZIMUTH) 

! SNSHD (SENSOR HEADING) 

0048 

0049 


HI ( 16)=0 

HI ( 17)-B2I (H. 51. 52)- 

!NREV (ORBIT REVOLUTION NUMBER) 
!EXYR (EXPOSURE YEAR) 

0050 

0051 


DY-Q2I (H. 47. 48) 
M0-D2I <H, 43. 44) 


0052 

0053 


HI ( 18)-UULIAN(MC. DY, HI( 17) ) lEXDAY (DAY-UULIAN) 

HI(19)»DY !DAY (SCENE ID DAY) 

0054 

0055 


HI (20)»A2I (H. 53. 54) 
HI (21 )*A2I (H. 55. 56) 

!HR (SCENE ID HOUR) 
!MIN (SCE.'JE ID MIN) 

0056 

0057 


HI (22)-0 
HI (23)-l 

!SEC (SCENE ID SECOND) 
'SEGNO(l) (TAPE SEQUENCE 4) 

0058 

0059 


HI (24)-l 
HR( 13)«0 

!SEQN0(2) (TAPE SEQUENCE 4) 
!FLAT (FORMAT CENTER LATITUDE) 

006L 

0061 


HR( 14)-0 
HR{ 15)*0 

IFLONO (FORMAT CENTER LONGITUDE) 
!NLAT (NADIR LATITUDE) 

0062 

0063 


HR( 16)-0 
RETURN 

INLONG (NADIR LONGITUDE) 


c 

c 

**«*♦«»**#«» 


• 

c 

c 

PROCESSING FOR LANDSAT 

«***•»*••*•* 

FORMAT TAPE 


c 

C MISSION NUMER 


ORIGIIViAL PAGE IS 
OF POOR QUALITY 


FORTRAN IV-PLUS V02-51E 14, 27: 3O-OCT-0O PACE 4 

hMOBS TTW TTft'DLfiBKff/wft r-^: 


2 

t 

0 

a 

m 

s 


0064 2 E«A21(H. 1, 1) 

666» - "■-" --iP/e, eq. — — ^ 

0066 IF(E. EQ. 2. OR. E. EQ. 6)E»2 

0667— - — XF<E. NE. 1, AND. E. NE.-2)E*3 - ^ TT 

0068 N»H<14> 

-IL.-".. ..':-TBKigff7Tr'^TOE"'SES0ENCE ii' ~ 

0070 WRITE(6. 13U)HI(23) 

0071 - ran FORMAT <TX. 'TAPE SEQUEN£6 NO. «'.12> I "" ! TT”- 

0072 N«H<16) 

0073 hT < 24 ) -I’AND < n 7 1 5 > ! SEGNO ( £ )' ( t APE SE6 uENCE N ) ^ 

0074 HI (6)-D2I(H, 39, 40)/4 INOSAHP <W SAMPLES PER SCAN) 

0075 - «I< 1 )»HK6)'*(HI<23)-1 )■*•! i^T (SCAN StAflT) “ 

0076 HI(2>»HI(1)-*‘HK6)“1 !SE (SCAN END) 

0077 - HK3)»l !LS (LINE START) 

0078 H: (4)«2340 'LE (LINE END) 

0079 HI(5)»4 ! NOCHAN (# CHANNELS ON TAPE) 

0080 HI ( 17)»A2I(H. 46. 47) !EXYR (EXPOSURE YEAR) 

0081 -- DY-AEKH, 41,42) - - - ^ 

C »*»•*** 

C • — THIS ROUTINE CONVERTS EBCDIC MONTH NAMES TO 
C INTEGER MONTH NUMBER 

C ■•*♦***** — • 

0082 DO 22 I - 1, 12 

0083 D0 2S-U- 1,3 - -- - 

0084 1F(H(J»42). NE. MON(J, I) ) GO TO 22 

0089 29 CONTINUE 

0086 GO TO 23 

0007 22 - CONTINUE 

0QB8 H( 18)»0 ! EXDAY IF NO DATA IN HEADER 

0009 GO TO 24 - - -• — - — 

0090 23 HI ( 1S)»-JULIAN( I, DY, HI ( 17) ) ! EX DAY (DAY-JULIAN) 

0091 - CALL ETAD(H(51), LATD, 1> 

0092 2 4 HR( 13)»LATL0N(H. 51. 96) ^ LAT (FORMAT CENTER LATITUDE) 

0093 HR( 14)-LATL0N(H. 59. 64) IFLONG (FORMAT CENTER LONGITUDE) 

0094 CALL ETAD(H(58). LONGD, 2) 

- ^ TEST FOR CORRECT LANDSAT LEVEL ■■ 

0099 IF ( (H( 139). AND. ”377). EQ. ■'3Q9)LVL«1 ! LEVEL 1 

3096 - IF (LVL. EQ. DGO TO 1000 

0097 IF ( (H(6>. AND. '‘377). EQ. ''140)LVL«3 ! LEVEL 3 

0098 -- IF (LVL. EQ. 3)G0 TO 3000 - - -- 

0090 LVL»2 

C ■ PROCESS LEVEL I 

C ETAB CONVERTS EBCDIC TO ASCII BYTE BY BYTE 

C BUILDS SCNID FIELD (SCENE ID) 

0100 1000 CALL ETAD(H( 1 ). HB( 13., 4) i 

0101 CALL ETAfl(H(6),HB(17),2) 

0102 CALL ETAB(H( 1 1 ). HB( 19), 2) 

0103 1100 DO 100 I»21,24 

0104 HD( 1)»''04 0 

0109 100 CONTINUE 

0106 IF (LVL. NE. 1)00 TO 2000 

0107 HI ( 13)»A2I (H. 101, 102) 'SUNEL (SUN ELEVATION) 

OlQg HI ( 14)^A2I(H, 106, 108) '5UNAZ (SUN AZIMUTH) 

0-109 - - HI ( 15)-A2I (H. 1 10. 112) ! SNSHD (SENSOR HEADING) 

0110 HI(16)»A2I(H, 114, 117) !NREV (ORBIT REVOLUTION #) 
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1 FORTRAN 

1 .'-PLUS V02-91E 14,27:29 

30-0CT-B0 PAGE 9 

HOCODE. FTN 

/TR: BLOCKS/UR 



om 


Hl(l9>-A2t(H. 14i, 14T5 


!^AY (SCENE ID 6aY) 

0112 


NX (20>«A2I(H. 147. 148) 


■HR (SCENE ID HOUR) 

0113 


HI(21)*A2I(H. 149. 190) 


•niN (SCENE I 6 MINUTE) 

0114 . 


Hl(22)-H<191) 



0119 


HI(22>-IAN0<HX(22>. 19)«10 


'SEC (SCENE ID SE'SUnST 

0116.. . 


HR ( 1 9 ) -LATLON < H. 68. 73 ) 


!NLAT (NADIR LATITUDE) 

0117 


HR( 16>«LATL0N<H. 79. 8' ) 


•NLONG (NADIR LONGITUDE) 

one 


RETURN 


- . 


C 

PROCESS LEVEL 2 



0119 

2000 

HI(13>»A2I(H. 108. 109) ! 

SUNEL 


0120 


HI(14)»A21<H. 112. 114) 


! SUNAZ 

0121 


HI(19)-A2I<H. 67, 69) 


! SNSHD 

0122 


HI ( 16)«0 


,‘NREV 

0123 - 


HX<19)«A2I(H. 144, 147) i 

DAY 


0124 


HI(20)»A2I(H. 149, 190) 


•HR 

0129 


Hn21)-A2I(H. 151, 152) 


•MIN 

0126 


HI (22>-H( 193) 



-- 0127 


HI<22>-IAND(HI<22), 19)»10 


: SEC 

0128 


HR( 15)«LATLQN(H. 77, 82) 


! NLAT 

0129 


HR< 16)*LATL0N{H, 84, 90) ! 

NLONG 


0130 


RETURN 




C 

PROCESS LEVEL 3 




C 

BUILD SCNID 



0131 

3000 

CALL ETAD<H< 1 ). HB< 13), 9) 


. ... 

0132 


CALL £TAB(H<6), H0(18), 2> 



0133 


CALL ETAQ(H<12>.HB(20), 1) 


, 

0134 


GO TO 1100 




C 

«**««*« 




C 

PROCESSING FOR UNIVERSAL 

FORMAT 

TAPE 


C 



. . „ . 

0139 

1 

H1(1>»1 !ST (START SCAN) 


0136 


HI(3)*0 !LS (LINE 

START) 


0137 


HI(4)-0 !LE (LINE 

END) 


0138 


HI(9)-H(90) ! 

NOCHAN 

(# CHANNELS ON TAPE) 

0139 


HI(6)=B2KH. 1787, 1788) 


INOSAMP (4 SAMPLES PER SCAN) 

0140 


HI(2)»HI(6) . . • ) 

SE (SCAN END) 

0141 


MISNO( 1 )-H(65) 



0142 


MISN0(2>*H(66) 



0143 


H0LD1>B2I (H, 67, 68) 


•SCNID (SCENE ID) 

0144 . 


DO 111 r=13. 24- 


_ 

0149 


HB( I )«“040 



0146 

111 

CONTINUE 



0147 


ENCODE (6, 200. HLDA)H0LD1 



0148 

200 

FORMAT! 1 16) 



0149 


CALL FRONT (HLDA. 12) 



0150 


ENCODE (6, 200, SCNID )H0LD1 



3T5I 


HI ( 13):^A21 (H. 2163, 2164) 


! SUNEL ( SUN ELEVA '' N ) 

0192 


HI ( 14)»A2I (H, 2169, 2170) 


! SUNAZ (SUN AZIMUTH) 

0193 


ril(15)*0 ! 

SNSHD (SENSOR HEADING) 

0194 


HI ( 16)»B2I (H, 71, 72) 


'NREV (ORBIT REVOLUTION 4) 

0199 


HI(17)*H(80) ! 

EXYR (EXPOSURE YEAR) 

0196 


DY=*H(79) 



0197- 


ENCODE (2, 250. EXDAY )DY 



0198 

250 

FORMAT! 12) 



0159 


SCNID(7)»EXDAY( 1 ) 



0160 


SCNID(8)-EXDAY(2) 
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original page 'S 



000 1 


•T.T.y: 


0003 




0003 


■.TST.7? 


0007 


INTEQgW FUNCTION Oai(H. BliB2) 
6 vTE HU > 

INTEGER Bli82, C. D 

0*H(B3> 


0SI-IAN0(C. ^S5)*S 


*9^0 tTV 


RETURN 
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At 


1 FORTRAN IV-PLUS V02-9IE 14-27:48 30-0CT-80 PAGE li 


— RW66E.TTN 

/tA: DLbCK§/W|i| 


C »«*»4 

C THIS 

»**•** 

ROUTINE CONVERTS DATA FROM BCD TO INTEGER 

0001 

c »*••< 

»***«• 

INTEGER FUNCTION A2KH. B1.B2) 

5oo5 

0003 


BYTE H<1> 
INTEGER B1.B2. C 

0004 

0009 


A2I-0 

DO 1 I • B1.D2 

0006 

0007 


C-H(I> 

C«IAND(C. 19) 

0008 

0009 

1 

IF(C. GT 9)C-0 
A2I-A2I*10-^C 

0010 

0011 


RETURN 

END 




OR'iGINA!- “ 
OF POOR v’-- 
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ORWINM. PA« » 
OF POOR QUALITY 


i FORTRAN IV-PLU5 V02-91E 14.27.49 30~0CT-a0 • PAQg 13 ■ 


HDCOOE. FTN 

“TtrTIlockSTwr = - - ~ 



C ••••< 
c ••••< 




C THIS ROUTINE COfIPUTES LATITUDE AND LONGITUDE 


0001 

0002 


REAL FUNCTION LATLONIH. 01. B2> 
IMPLICIT INTEGER (A-Z) 


0003 

0004 


BYTE HB - - . - .... 

BYTE H<1) 


0003 

0004 


BYTE N. S. E. H. D - 

DATA N. S. E. W. D/213. 226. 197. 230. 96/ 

• 

0007 

0008 


LL « 0 
SIGN • 0 


0009 

0010 


DO 1 I - Bl. 02 , - 

HO«H< I) 

“■ 

0011 

0012 


IF(HD. EG. N) CO TO 1 * 

»F(HB. EG. S) GO TO 3 


0013 

0014 


iF(HB.EG. E) GO TO 1 
IF(HB. EG. W) GO TO 3 

• 

0013 

0016 


IF(HB. NE. D) GO TO 4 * " • 

GO TO 1 


0017 

0018 

3 

SIGN - 1 
GO TO 1 


0019 

0020 

4 

J - HD 

J - IAND(0. 15) 


5351 

0022 

1 

LL ■ LL*10+J - - 

CONTINUE 


5353 

0024 


LM » LL - (LL / lOa > * 100 - fMINUTES CALCULATION 

LL « LL / 100 ! DEGREES CALCULATION 

5353 

C 

LATLON - FLOAT < LM ) / 60. -p FL(jAT( LL ) 

LATLQN CONVERTS DATA TO REAL VARIABLE FOR RETURN 

0026 

C 

SIGN IS neoativs tor s or w direction 

IF< SIGN . NE. 0 ) LATLON « - LATLON 

. .. . 

5557 

0028 


REtURN ^ 

END 













• 
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I FORTRAN IV-PL^S yoa-9lE 14 2* 41 30-0CT-80 PACE 1 

CBCOOrTTN ■ /TR DL6cRI7wfl 


la 

II 

Id 

» 


£ 

c 

i 

c 

£ 

c 

c 

T 

c 

£ 

c 

c. 

c 

c 

c 

c_ 

c 

c. 

c 

c 

c“ 

c. 

c 

£, 

c 

c 

? 

c 

c 

c. 

c 

c_ 

c 

c 

c. 

c 

c 

c 

c_ 

c 

c 


LOCF.HS ED electro n I C5 COMPANY, INC. 

RENOTE SENSINR APjPLlCATlONi LAB. 
< R S A L ) 


NAN£: LOCOOE 

DEVELOPED BY: SUSANNE OBRIEN 
ON: 6/78 


REVISED DY ON' PURPOSE FOR REVISION: 

o<>ooo< o< ooooooooooooooooooo< 


PROORAM DESCRIPTION: 

<sooooo«;:>ooo<x:!0<>0v'^ooo<:!000oooooo<><>< 

THIS ROUTINE WILL DECODE VIDEO DATA FPOH A FOREIGN IMAGERY 
TAPE IN ONE OF THE FOLLOWING FORMATS: UNIVERSAL. LANDSAT 1 OR 2. 
OR LARSYS 2 OR 3 : 

ENTRY POINTS: 

<><><>0<>00<>00000<X><><><>0<>00<><>0000<>00< 


PROGRAM INITIATION OR CALL SE Q UENCE' 

oooo<><>o<><>o<><;>oo<:Toooo<>o<><><>oo'o<.><><><>< 

CALL LDCODE(LUN. IF, BUFFER. DFSZ. L. RCHAN, LSS. LSE, BUF, EOF, 

1 PC, IE. INlf) ■ 


CALL ARGUMENTS: FORMATS. 

ooooc^oo oo»:x><>ooooo<>oo<>o<>ooo< 

LUN . II LOGICAL UNIT NUMOER 

IF IF FORHAT OF THE INPUT TAPE 


BUFFER 

iFSl 

L 

RCHAN 

l^S 

LSE 

BUF 

EOF 

PC 

IE 

INIT 

'INPUTS: 

o<^oo<^<>o 

HCOM 

OUTPUTS: 


1»UNIVERS AL 

S-LANDSAT 

3*LARSY5 

BYTE'sTRING BUFFER FOR INTERNAL USE 
OF LDCODE 

11 SIZE OF DUF IN BYTES ' 

n REOUES1ED LIME NUMBER 

II REQUESTED CHANNEL NUMBER 
II FIRST PIXEL TO MOVE 
‘n L/Cf PIXEL IN MOVE 
BYTE STRING BUPCER FOR DECODED PIXELS 
END OF FILE l=EOF 

11 PARITY COU NT 

II ERROR CODE ' 

11 FIRST PASS SWITCH Q«FIRST PASS 

FOKMATS 

COMMON BLOCK OF DATA FROM HPROS 

formats 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


1 FORTRAN 

IV^PLUS V02-51E 14:26:41 30-0CT-e0 PAGE 

2 

LWOCe. F^N 

/TR: BLOCKS/WR 



C 



c 

c 

SPECIAL REMAf^KS: 

ooooooooooooooo oooooocxxxxxxxxxx 


c 

c 

1. THIS PROORAH ASSUMES THE MAIN ROUTINE WILL ALWAYS 
FOR DATA BEGINNING WITH LINE 1 AND PIXELS BEGINNING 

ASK 


T" 

c 

WITH PIXEL I AND CHANNELS BEGINNING WITH 

CHANNEL 1 REGARDLESS OF WHAT THE ACTUAL VALUES FOR THE START 


c 

c 

OF THOSE VARIABLES ARE. 

2. THE PROGRAM HOREAD MUST BE CALLED BEFORE 



c 

c 

LDCOOE OR ELSE THE COMMON AREA WILL NOT BE SET UP. 



c 

c 

ROUTINES CALLED; 

OOOOOOOC- 'XXXXXXXXXXXXXXXXXXXXXXX 


c 

c 

MTREAD READS MAG TAPE RECORDS 

IBYTE converts DATA FROM ASCII TO INTEGER 


C 

C 


c-< 

c 




0001 


SUBROUTINE LDC0DE(LUN1, IF, BUFFER. BFSZ, DLIN, RCHAN. LSS. 
IE. INIT) 

LSE, BUF. EOF. PC. 

5552 

0003 


IMPLICit INTEGER <A-Z) 
INTEGER LSLOCO) 


0004 

0009 


BYTE BUFFER < 1 ). ISTATB<2) 

BYTE BUF( 1 ), SCNID, M1SN0(2) 


0006 

0007 


INTEGER IB(64>, JR«64> 
DIMENSION I3TAT(2). i:’RM<6> 


5553 

0009 


EQUIVALENCE ( iSTAT< 1 ) . ISTATD ) 

COMMON /HCOM/SS. SE, LS, LE, NRPDS. NDSPR, NCPR. NRPC. ANCL. NC. NS. 

0010 


1 NBIT, DOI, NCAR. SVD, RSIZ, PSKIP, HSIZ, CALP. CFRR. NSPR 
C0MM0N/LASEL/NAME(3). MISNO, E 


5oTi 

0012 


DATA LSLOC/71. 0, 1/, RLB/"1000/. IPRM(2) /32B0/ 
lEEOF— 10 


55T5 

0014 


IEPC«-4 
NAME ( 1 )«'LD' 



MT5 >JAME<2)-'C0' 

0016 MAME<3>-'DE- 


05l7 1F(LSE. Gt. SE)LSE»SE !CHECK SCAN END 

00 1 8 IF( INIT. NE. 0> GO TO 1 -CHECK FIRST PASS 


0019 

0020 


INIT ■ 1 
PRS2»(RSIZ*1 )/2 

ISET FIRSi PASS INDICATOR OFF 
!DATA SIZE IN WORDS 


C 

PUTTING STARTING ADDRESS OF BUFFER INTO IPRM ARRAY 

0021 

c 

CALL GETADR< IPRM, DUFFER ; 

0022 

c 

DSL»ANCL-*-NC*(NS-^CALP) 'DATA SET LENGTH 

0023 

0024 


ADD-=0 

ONEST*! !SET FIRST RECORD 

0029 

0026 


CL«1 'CURRENT 

CCA,N«1 'CURRENT 

LINE 4 
CHAN N 

0027 

0028 

9 

ANC ■ ANCL-^GVD 
IF< IF, EQ. 1 ) ANC « 

! TOTAL ancillary BYTES 
ANC 2 'FORMAT »UNIVERSAL 
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OF POOR QUALITY 


POBTR/,^ 1V-PLU8 V02-51E 


LDCOO'.. PTN /TR: BLOCKS/WR 


30-0CT-B0 


RAOC 3 



FUF EQ. S) GO TO 1 'FORMAT aLANOSA 

THIS ROUTINE BUILDS ARRAYS IB AND JR — 


IB<l)»OFFSET BYTES TO CHAN I 

JR<I)»RECORD » WHERE CHAN I IS FOUND 


IF <NRPC. OT. 1)00 TO 100 'RRECS/CHAN >1 

F<NCAR. QT. 0)00 TO 10 'DATA IN ANCIL REC 00 


THIS ROUTINE FOR NO DATA IN ANCL REC FIRST CHAN IN REC 

«< 1 > 


I 0 < 1)»2 
0 140 1«2. N 


JR<I)aJR<l-l)t>l 
IB<I)a2 


CONTINUE 

00 TO 1 ” - 


FC - 1 
LC - NCAR 


• 1 

IR • 1 


DO 12 CHAN ■ l.NC 

CONTINUE - - _ 


IFdR. OT. 1) ANC • 2-fSUD 

IF<CHAN. OE. FC. AND. CHAN. LE. UC5 GO TO 14 


IF<CHAN. LE. LC. OR. IR. OE. NRPOS) GO TO 19 

THIS ROUTINE RESETS VARIABLES FOR NEU RECORD NEEDED 


- LC ♦ 1 

LC • LC ♦ NCPR -• - . 


IR - 1 
CO TO 13 


THIS ROUTINE CALCULATES OFFSET BYTES ANDRECORD « 
IB<I)-(CHAN -‘FC)*(NS+CALP)+ANC-1 


JR< 1)«IR 

I » I ■*. 1 


CONTINUE 
00 TO 1 


CONTINUE 

CERR « 4 lERI^OR TOO MANY CHANS OR RECS 


GO TO SO 

THIS ROUTINE SETS UP ARRAYS FOR MULTIPLE RECS/CHAN 


IFINCAR. GT. 0)00 TO 120 
JR<1)«2 


IQ( 1 )«2 
DO 130 I>2>NC 


JR<I)-JR<I-l)-*-NRPC 

IB<I)-2 


CONTINUE 
00 TO 1 


JR< 1 )«1 

ID< I )»ANCL>SVD-l 


GO TO no 

START IF NOT FIRST PASS OR IS LANDSAT 


IFIRCHAN. LE NOGO TO 31 ! REQ CHAN TOO BIO? 
CONTINUE 


CERR>6 !YES ERROR 
GO TO SO 


CONTINUE 

IF<CL-DLIN)29, 17, 18 iCURRENT LINE»DESIRED LINE 








ORIGir.’Al , ^ r? 

OF POOR OUALiTY 


I FORTRAN IV-PLU8 V02-51C U: 8A: 4i 30-0CT»B0 PMt 4 


LDCODE. FTN 

/TR BLOCHS/UR 


0071 

li 

6&nTinue 

URITE(6. llinCL, DLXN 

- 1- ’ 

5ot5 

ooao 

iUi 

FORMAKIX, 'CL*'. 19. 'OLIN*'. tS) 
CERR-7 !N0 ERROR CL>DL 

• 

jmiijmii 

c 

66 T(S it 

ROUTINE FOR NEED NEXT LINE CL<DL 


0062 

t 

29 

» o/^Ya r£66ro < or >i. 

IF <NOSPR. LE 1 )G0 TO 43 



6 

C 

ROUTINE FOR » DAtA SETS f^&R RECORI 
CALC FIRST LINE OF DESIRI’D REC 

Ti 

5313 

0084 


flin»dlin-modTdlxnhone1TTn61P1T 

ADD«<DLIN-FLIN)»DSL 


5315 

0066 


L1M»<FLIN-CL)/NDSPR 

IFUIODC (FLIN-CL). ND5PR). NE. 0)LIM>LII1-t>l 

0087 

0088 


IF<LIt1. EQ. 0)00 TO 190 
GO TO 44 


0089 

t 

43 

CONTINUE 


55?T5 

C 

LIM-<Di.iN-CL-l )»Nf^Fi6s-*-i 



C REAi!>INd DAtA FROM TAPE REAl!) LIM RE(^0A1>E 
C 



6091 

0092 

44 

00 119 I»1.L1M 

CALL UTQIOIRLB.LUNI. 1. . ISTAT. IPRM 

ISU) 


0093 

6 

IF<1STATQ(1).EQ. "001)00 TO 119 

•READ 0KAY>1 


0094 

0095 


Ff < 1 stATB i 1 ) ; Eg." I eeOE ) eoP- i 
IF<EOF, EQ. 1 )RETURN 

(End of «^ile-" 366— lo 


”35^6 

0097 


IF< ISTATBl 1 ). EQ. lEf^C )PC»PC-*-l 
1F< ISTATBl 1 ). EQ. lEPOGOTO il9 

•I^ARltY 6HECK-"374— 4 


“3591 

0099 

119 

CONTINUE 




C 

C 

SET CL TO DESIRED LINE AND CCAN 

TO 1 


0100 

c 

190 

CL»DLIN 

• 


0101 

C 

CCAN » 1 

THIS ROUTINE FOR CURRENT LINE » 

DESIRED LINE CL-DL 


0102 

17 

C 

IFINRPDS. LE. 1 ) GO TO 23 

IF NO THEN DETERMINE IF REQ REC 

(PRECS/DATA SET <-l 
IN DUFFER 


0103 

C THIS CHECK FOR NO DATA IN ANCIL REC 

IF(RCHAN NE. DGO TO 170 'REQ CHAN NOT ■! CO 


0104 

0109 

170 

IF< <LSS, EO. 1 ). AND. (NCAR.EQ. 0)) 
1 > vIR(RCHAN)-JR(CCAN) 

GO TO 180 


0106 

0107 


CCAN ■ RCHAN 
IF< I. EQ. 0) GO TO 23 

•REQ REC IN BUFFER 


0108 

C 

180 

CONTINUE 


1 


C 

C READING DATA FROM TAPE 


0 

0109 

C 

CALL WTQIOIRLD. LUNl. 1. , ISTAT. IPRM 

ISU) 


0110 

C 

1F< ISTATB< 1 ) EQ "001 >00 TO 23 



01 1 1 
01 12 


IF( 1STATD( 1 > EO. IEEDF)E0F--*1 
IF(EOF EQ DRETURN 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


i 

i FORTB^iN IV-PLU3 V02-91E l<:S4;41 30»0CT-a 0 PAQt 5 


LDCODE. FTN 

/TR: BL0CK8/UR 

0113 

, <?ll* 


lf<ISfATB<2>.^0. IkPOFC-f^Col 
IFdSTATBm. EQ. lEPOfiD TO 23 

0119 

0114 

S3 

GOTO 40 'READ ERROR 

IF(XF £02) 00 TO 29 !LAN08AT7 

0117 

c 

THIS ROUTINE MOVES DATA ^OR UNIVERSAt/LARSYS FORMAT 
. IFINRPC. OT. 1)00 TO 200 'MULTIPLE REC8/CHAN 

one 


B - IB<RCHAN)4> L88 ♦ADO 

B - B^ LSE- L88 

oiao 

0121 


K - 0 

00 24 I • B.E 

01» 

0123 

24 

K ■ K ♦ 1 
BUF<H>>BUFFER(I) 

?T24 

C 

60 tO 90 

THIS ROUTINE MOVES DATA FOR LANDSAT FORMAT 

0124 


B«< < LSi-1 >/3)*8^RCHAN»2~M0D( LSS. 2) 
E- ( < LSE- 1 > /2) •8vRCHAN»2-M00 ( LSEi 2 > 

515? 

0128 


J ■ 1^M0D( LSS. 2> 
K ■ 0 

5T5? 

0130 

id 

5a"5? I - 1,0 

K - K ♦ 1 

0131 

0132 


BUF<K)»BUFFER<D) " " 

B ■ B ♦ 1 

5T33 

0134 

27 

CONTINUE . .. 

IF(B. 01. E)00 TO 90 

“5159 

0134 


J ■ 

B - B ♦ 4 

0137 

c 

00 TO 28 

THIS ROUTINE FOR MULTIPLE RECS/CHAN 

0138 

0139 

200 

IF(L5S. NE. 1)00 TO 210 

Z-NSPR IZ-CONTROL VARIABLE FOR ROUTINE 

o o 
— o 


Y-LSE !Y»C0NTR0L VARIABLE FOR UPPER LIMIT 

IF<LSE. OT, NSPR)V«NSPR 

0142 

0143 


B-1B(RCHAN)^LSS 

E-B+Y-1 

0144 

0149 

230 

K«0 

DO 220 1-B. E 

0144 

0147 


K-K+i 

BUF ( K ) -BUFFER ( I ) 

0148 

0149 

220 

CONTINUE 

IF (LSE. LE. Z)00 TO 90 

C 

C READ DATA FROM TAPE . 

0190 

C 

CALL WTQI0(RLB. LUNl. 1. . ISTAT. IPRM. ISM) 

0191 

0192 


IF ( ISTATBd ). EQ. "001)00 TO 223 
IF(ISTATB(l).EQ. TEE0F)E0F-1 

0193 

0194 


IF(E0F EQ. 1) RETURN 

1F< ISTATB( 1). EQ. lEPC )PC«PC ♦•1 

0199 

0194 


IF < ISTATB< 1). EQ. lEPCXJO TO 223 

GOTO 40 (READ ERROR 

0197 

0198 

223 

Z-Z+NSPR 

IF(LSE. LT. Z)Y-LSE-(Z-NSPR) 

0199 

0140 


D-3 

E»2>Y 

0141 

C 

GO TO 230 

THIS ROUTINE FOR SECOND PASS PROC FOR NRPCM 
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SNOI 1 3~>r- nvyr >tu 





i 

t 


J 

i 
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APPENDIX D 


LISTING OF POSTP PROCESSOR 



ORIGinn ' 

OF POOR QUALIFY 


1 FORTRAN IV-PLUS V02-51E 14-33.33 30-OCT-80 PAGE 1 

POSTPP, ?tN 

/TR DLOCKS/WR 

0001 

c 

PS564aM *^OSTf*t* 


Z YHT¥T*RdtfRAM x^’mf^esrFATCTr’PAKeSS f MS’' WDCeSSOft 

C IT COMPUTES THE SPATIAL RESPONSE FUNCTION BY REORESSIOf. 

Z Srito A QUAbAAttC~SURfA6£ 


c 


— SSCT" 

0003 


"•isyt-g-pn:’NAMT55T:'irfepflff^^ yr.i5 ri"5~) WAiie < 6 > 

DIMENSION I0AR(16. 16). 1CL0UD(16, 16), IU/i7ER(16. 16) 

5554 


DIMENSION VINO, 16, 16). IGIN(3. 16, 16 ) . vInHI 16, 16, 2) . R ( 16. 16) 

0003 


DIMENSION P(6), C<6). X<4) 


c 

c 

READ FILE NAME FROM INPUT CARD 

0006 

d 

READ<3, 1000)FILNAM 

555? 

1000 

(30A1) 

0008 


DO 10 1-1.30 

*5555 


iF(FILNAH( i)7Ne. ' ')G(j tfl l5 . 

0010 


FILNAM(I)»0 

55T1 


QQ Jj 

0012 

10 

CONTINUE 


c 

c 

READ REPORT TYPE REQUEST FROM CARD SHOULD QE LONG OR SHORT 

0013 

0 

13 

RCAD<3. 1200)REPCRT 

0014 

T355 

FoAmaTF4ai) 

0013 


IPASS-1 


C 

c 

OPEN FILE AND READ HEADER AND DATA 

0016 

t 

0PEN<UNIT-1. NAME-FILNAM, TYPE- 'OLD', FORM- 'FORMATTED ' ) 

55T7 


Aead( 1, iooo) ID, iaO, iSAt, la1dA, latd, laYm, lCnDA. lOnD, uonm, NaH^ 

0018 

2000 

FORMATd.X, lOAl. 12, 1 1, Al, 14. IS, Al. 14, 12. 2X.6A1) 


■S3T^ DO sd I«l, 16 

0020 DO 2S J-1,16 


TJ551 READ ( 1 . 1 1 00 HPATCH, <X(K), K»l. 4 > , IOaR ( I , J ) , ICLOUDl t, J)/ 


1 IWATERd. J), GAMA, ISLN, (UINCK, 1. J),K-1. »), i IGIN<K. I, J). K»l. 3) 


— 6o25" 

0023 

1100 

23 

FORMAKIX. 13. 4F6. 2, 314, F6. 3. 13. 3Jif4. 1, 2F5. i. 3l3) 
CONTINUE 


C 

C 

SET UP RN»* CONSTi^NTS REGRESSION CONSTANTS 

0024 

C 

RNOO-16 

5553 

0026 


RNO 1-5440. 5 

RN02-91392. **. 5 

0027 

C 

RNl 1-115600. **. 5 


C 

c 

CALCULATE C VALUES 

0028 

0029 

30 

DO 40 I-;\,6 
C d ) -0 

COSO 

0031 

40 

CONTINUE 

IGOOD-0 



DO 100 I-l, 16 
DO 100 J-1, 16 








ORIGINAL PAGE iS 
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i FORTRAN 1V-PCU8 V02-SIE 14. 33 33 30-OCT-iO PAOC 3 

X ' 'v ■ - 7 ' * ■"' " ■^'■1 " ■■ m m ■. T7 , ■ ■ — ,i — l,l■ll, ,i, — , i^i — ».i t 


POSfpk 

1 A V r WV9 

^TN 

VV» 9AC. «JW WV 1 «W ■ 

/TR DLOCK«i ' "" 

5535 

OP3A 


li^UH-tdAfld. j)-^I5lcIuo< (. J)«lwAr^«^( it. 3) 
IFdSUn OT 128)00 TO 100 

owf 
P03i _ 


iodbo-ioodb^i 

IPdPASS eo 1 )VAR>IOtF‘d. I. Jl 

0039 

C 

fPdPAili 3)VAR>ViNH(l. 1 ) 


t 

c 


5545 

0041 


Pd)«l/RNOO 
P<2)-l/RN01*d-8. 5) 

5543 

0043 


■><3>-l/^N0t*< j-a. 5) 
P<4)«l/RN02*< ( 1-0 5) **2-21. 25) 

5544 

0045 


P(5)»l/RM02*( < J-8 5>**2-21 25) 
P(6)-l/RNll*d“-9 5>*<J-0. 5) 


c 

c 

C<K) RAW ORTHOGONAL C0EFFICIENT3 

0044 

c 

00 SO K«l. 6 

5543 

0048 

50 

C<K)»C<K)*P(K)*VAR 

CONTINUE 

0049 

0050 

100 

CONflNUi 
00 150 K*1.6 

o55i 

0053 

150 

coo*CiK>*'25Zr7T5o3D 

CONTINUE 


C 

c 

CALCULATE ViNHAT. R, SSR. AND MSE 


c 

ViNH UINHAriS r‘H£ ESTIMATfeb VA*^iABLE VALUfi 
SSR IS THE SUM OF THE SQUARES OF THE RESIDUALS 


5 ' 
c 

“Sm5e~ 71“ fHE“ R00t MgAN SQUARE EhROi^ 

5553 

00'.^< 


SUM-5 

SSR»0 

CO 5 5 
C054 


DO 2C0 1*1. 16 
00 200 J-1 16 

0057 

0058 


IFdPASS. EQ 1 )UAR-10IN< 1> I. J) 
IFdPASS. EQ. 2)VAR-VINH< I, J. 1) 

0059 

0060 


Pd)-l/RN00 ■ 1 

P<2)«1/RN01'^< 1-8. 5) 

0061 

0062 


P<3)-l/RN01*(J-a. 5) 
P<4)»1/RN02*< d-8 5)**2-2I ?5) 

0063 

0064 


P<5)«1/RN02*( < J-8. 5)**2~21. 25) 
P<6>»l/F,Nll*d-8. 5>*(J-9 5) 

0065 

0066 


DO 175 K-1. 6 

VINMd, J. IPASS)»UINHd, J. IPASS) ■►C <K)*P<K) 

0067 

0060 

175 

CONTINUE 

Rd. J)»UAR-VINHd, J. IPASS) 

C^e>9 

0070 


IF< IPASS EQ. 2)Q0 TO 178 
1 SUM* I GAR < I. J>-t-ICLOUD( I. J)-*>IWATER( I. J) 

0071 

0072 

17S 

IF< 1SU11. ST. 128)&0 TO 200 
S5R»SSR-*-Rd, J)**2 

1 0073 

» 0074 

200 

SUM»SUM-H<I, J) 

continue 

1 0075 

' 0076 


RMSt'^SS" 'dCOOD-6) 
RMSE*RM-e** 5 


c 

c 

CALCULATE CN»* VALUES COEFFICIENTS FOR THE ORTHOGONAL POLYNOMIALS 
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1 FORTRAN 

1V-FLU9 

V02*9tC 14-33-33 30-0CT-80 FACE 3 

PMWTFTN 

7trt 'iLOOKi/urt 

0077 

C 

CNOO-^< 1 >/RNOO 

007» 


TF35T551TT7IFS3I 

CN10«C(3>/RN01 

5513 

OOtl 


0no^«0 1 4 i /kN08 
CN30"C<9)/RN02 

— 5513 

e 

‘CT7rT=5Tirr7iwn 


c 

c 

C*LCUi;ATE“»r V»UUf § 

COEFFICIENTS FOR SIMPLE POLYNOMIALS TO BE INPUT TO 8A0 

00S3 

t 

8A81-CN00-8 9»<CN01>CN10>*91 *<CN05^CN20)+72. 29*CN1 1 

558^ 

008S 


SASY-CNlO-17. *CN20-B 9*CN11 

0004 

0007 


gA8k^-CN02 

6A8Y2-CN20 

35BB 

c 

SAS>iY--ORn 


c 

c 

URttfe SHffffT'ftEPflAT INPO 

• 

0009 

c 

IFdPASS. EQ. 2)60'T0 600 

d090 

0091 


IAA88-il*ASs+i 
WRITE<6. 1300) ID. NAME 

5W3 

0093 

T355 

^IHl. 4'oy, 'AAtCH FOft ScSnE', ix. lOAl, 2 X. 'FILE; TTATTTTr 

WRITE(6, 1400) lAO. IfaAT. LATDR, LATO, LATM. LONDR. LOND, LONM 

0094 

1406 

■pl«lMAT(iX. 'OUN 'AKlOLfe . Iti. ix, '§ATCLLIT8 NUMSER ; '.'19, 

1 'LATITUDE IX, Al. IX. 14. 'DEGREES', IX, 12, 'MINUTES', 

0099 


rSx, 'LdiNOltUDE : , IX. Al. IX, l4, ' 6fGrtSES^ IX, 12, 'HlNUtSS ',/> ) 
WRITE<6. 1900) <C(K), K-1, 6). RMSE 

53^'S 

i 500 

f6/?mat< u', 'COO-', ^6‘S; 3xr'"ci3i-'. f^S 2, 3x, 'Oib-7 fe. 2, Ox. 

1 'C02"'. F8. 2, 3X, 'C20-', F8. 2, 3X, 'Cll-', F8. 2, 3X, 

009. 


i 'MSg^i'.f^B, 2) 

WR17E<6. 1990)CN00, CNOl, CNIO, CN02, CN20, CNll. SUM 

5o9ff " 

O 

n 

F0RMAT<1X, 'CNOO'.Ffy, 4, 3X. 'CNOl'.FB 4, 3X. 'CNIO'. 
1F8. 4, 3X, 'CM02', FP. 4, 3X. '.CN20', F0. 4, 3X, 'CNll 

0099 


1 F8 4. 3X, 'SUMSR ', F8. 2) 

WRITE<6, 1979)SAS1, SASX, EASY, SASX2, SASY2, SASXY 

0100 

1979 

F0RMAT<IX, 'SASl ', FB 4, 3X, 'SAGX'.FH 4, 3X. 'SASY '. 

1 F8 4, 3X, 'SASX2', F8. 4, 3X. 'SASY2', FS 4, 3X, 'SASXY ' . 

0101 


2 F8. 4) 

IF(REP0RT(1 ). EQ 'S')00 TO 900 


C 

C 

WRITE LONO REPORT DATA 

0102 

C 

URITE<6. 1600) 

0103 

0104 

1600 

FORMAT(//, lOX, 'PATCH', lOX. 'VIN', lOX, 'VlNHAT'. 12X. 'R ' ) 
DO 400 I»l, 16 

0109 

0106 


DO 400 Jal, 16 
IP«(1-1)*I6+J 

0107 

0109 

1700 

WRITE <6. 1700) IP, IGINd, I . J) , V INH< I , J, 1),R<I. J) 
FORMATdlX. 13, llX. IS.BX.FB 2, 8X, F8, 2) 


400 

900 

CONTINUE 

CONTINUE 


C 

c 

USE THIS GO TO FOR 2 PASS RESULTS FOR CHECK UUT OF ALGORITHM 
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i FORTRAN 

IV-PLUS 

V02-91E 14;33;.33 3O-OCT-0O PAGE 4 

plSstpp. 

1 

FTN 

/TR. OLOCK^/WR 

out 

C 

6o Tb 3d 

GO TO 700 


C 

C 

URITE REPORT FOR PASS 2 

0112 

C 

600 

MRITE(6» 1800) 

0113 

0114 

1800 

FORMATdHl) 

URITE(6. 1900) (C<K).K«1. 6). RMSE 

0119 

0116 


URXTE(6. 1590)CN00« CNOl. CNIO. CN02. CN20* CNll. SUM 
1F(REP0RT<1). EG. 'S' >00 TO 700 

0117 

0118 


WRITE (6. 1600) 
DO 650 I-li 16 

0119 

0120 


bb 6^6 ~J-1< 16 
IP»<I-1)*16+J 

0121 

0122 

1790 

WRITE<6. 1750) IP. VINH( I. J. 1) . VINH( 1. J. 2) . R ( t. J) 
FORhAT< IIX. 13, 6X. FS. 2, 8X, F8. 2, BX, F8. 2) 

0123 

0124 

650 

700 

CONTINUE 

CONTINUE 

0125 


END 
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• o: a; 


IN 4 
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<o ^ 
o 
o 

S3 
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u . 

ICi < 
|0 Q 
U X 
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ju 
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lOi 
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rd 


n 
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POSTPP. FTN /TR-rOLOCK^/WR 


FUNCTIONS AND SUBROUTINES REFERENCED 
OPENt 


TOTAL SPACE ALLOCATED - 034214 7238 

, TT17: C320# 403A. D/-SP-SYO: C320. 403P03TPP 
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